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CHAPTER 1 
INTRODUCTION 



1 - 1 . General. 

cl This series of manuals, entitled Designing Fa- 
cilities to Resist Nuclear Weapon Effects, is or- 
ganized as follows: 

TM 5-858-1 Facilities System Engineering 
TM 5-858-2 Weapon Effects 
TM 5-858-3 Structures 
TM 5-858-4 Shock Isolation Systems 
TM 5-858-5 Air Entrainment, Fasteners, 
Penetration Protection, 
Hydraulic-Surge Protective 
Devices, EMP Protective 
Devices 

TM 5-858-6 Hardness Verification 
TM 5-858-7 Facility Support Systems 
TM 5-858-8 Illustrative Examples 

A list of references pertinent to each manual is 
placed in an appendix. Additional appendixes and 
bibliographies are used, as required, for documen- 
tation of supporting information. Pertinent biblio- 
graphic material is identified in the text with the 
author's name placed in parentheses. Such biblio- 
graphic material is not necessary for the use of 
this manual; the name and source of publications 
related to the subject of this manual is provided 
for information purposes. 

b. The purpose of this series of manuals is to 
provide guidance to engineers engaged in designing 
facilities that are required to resist nuclear weapon 
effects. It has been written for systems, structural, 
mechanical, electical, and test engineers possessing 
state-of-the-art expertise in their respective dis- 
ciplines, but having little knowledge of nuclear 
weapon effects on facilities. While it is applicable 
as general design guidelines to all Corps of En- 
gineers specialists who participate in designing 
permanent military facilities, it has been written 
and organized on the assumption a systems- 
engineering group will coordinate design of the 
facilities. 

c. Technical Manual 5-858 addresses only the de- 
signing of hardened facilities; other techniques to 
achieve survival capacity against nuclear weapon 
attacks are deception, duplication, dispersion, 
nomadization, reconstitution, and active defense. A 
facility is said to be hardened if it has been de- 
signed to directly resist and mitigate the weapon 
effects. Most of the hardening requirements are 
allocated to the subsidiary facilities, which house, 
support, and protect the prime mission 
materiel/personnel (PMMP). This manual is ap- 



plicable to permanent facilities, such as those asso- 
ciated with weapon systems, materiel stockpiles, 
command centers, manufacturing centers, and com- 
munications centers. 

<L The nuclear weapon threats considered are 
listed below. Biological, chemical, and conventional 
weapon attacks are not considered. 
—Weapons aimed at the facility itself or at 

nearby targets 
—A range from many, relatively small-yield 

weapons to a single super-yield weapon 
—Weapon yields from tens of kilotons to 

hundreds of megatons 
—Weapon delivery by aerial bombing, 

air-to-surface missile, surface-to-surface 

missile, or sattelite-launched vehicle 
—Detonation (burst) of a weapon in the air, at 

the ground surface, or beneath the ground 

surface 

—Direct-overhead bursts for a deep-buried 
facility 

—Near-miss bursts for a near-surface facility, 
producing peak over-pressures from tens 
to thousands of psi at the facility 

e. The designing of facilities resistant to nuclear 
weapon effects is an evolving speciality using a 
relatively narrow data base that incorporates both 
random and systematic uncertainties. The range of 
these uncertainties may vary from significant 
(order of 1 to 2 magnitudes) to normal (10 to 100 
percent variation from average values). The ap- 
plicable uncertainty value depends on the specific 
weapon effect or hardening objective under consid- 
eration. Loading uncertainty is generally more sig- 
nificant than resistance uncertainty. Awareness of 
the appropriate uncertainty (extent of ignorance) 
factor is essential not only for system engineering 
trade-offs, but in the utilization of available analy- 
sis or test procedures. Studies and experiments are 
being conducted to improve methodology, to better 
define random uncertainties, and to reduce system- 
atic uncertainties. This manual will be revised as 
significant improvements occur in either method- 
ology or data base. 

1-2. TM 5-858-6: Hardness 
verification. 

This volume presents methodology for verifying 
the hardness of the facilities. The methodology 
comprises both physical and mathematical simula- 
tions. A synopsis of available testing facilities is 
given. 
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CHAPTER 2 
GENERAL METHODOLOGY 



2-1. Introduction. 

cl A hardness verification program for protective 
facilities and their systems, subsystems, and ele- 
ments measures the ability to survive an attack 
that has been specified in terms of weapon size, 
range, and burst conditions and of particular site 
characteristics. When nuclear test ban treaties pre- 
cluded atmospheric nuclear testing, alternative 
testing techniques and analytical procedures were 
developed to assess the probability of a facility 
surviving a prescribed attack. Most protective fa- 
cilities are extremely complex, containing major 
systems that require methodical testing at the sub- 
system level and integration of the results; this 
analytical methodology is called hardness verifica- 
tion. 

b. Hardness verification must be made at regular 
intervals during the facility develop- 
ment/design/construction cycle. Simple 
verifications will be performed during the initial 
design phases. As designs become firm and fabrica- 
tion begins, more comprehensive verifications will 
be made. It is essential that hardness verification 
be used periodically to monitor the iterative pro- 
cess of definition, systhesis, design (redesign), ana- 
lysis, test, and evaluation that transforms mission 
requirements to demonstrable and acceptable fa- 
cility survivability. This will secure continuing in- 
tegrity between design and construction, will verify 
that each element, subsystem, and system has been 
combined in a manner that properly accounts for 
uncertainties in responses and the interdependence 
of various physical parts of the facility. 

c The analyst must know the failure modes and 
resistances for the relevant elements of the sub- 
system or higher level assemblages, and the local 
weapon-effect loads on those points. Whether from 
experimental tests or from analyses, the data pro- 
vided must include uncertainties associated with 
each mean value. This volume deals primarily with 
procedures to evaluate failure modes and resis- 
tances; techniques to deal with the random nature 
of the process are discussed in chapters 5 and 6. 

2-2. Identification and organization of 
system elements. 

a. The first task in system hardness verification 
is to define the physical system. The relationship 
of each element to mission critical functions must 
be accurately defined and the associated local loads 
and failure modes must be specified. For complete 
examination of a facility, network logic analysis 



including the use of Boolean algebra is strongly 
recommended. The construction of fault trees 
within the context of multilevel system organiza- 
tion is a procedure will suited to the solution of 
complex problems dealing with verification activi- 
ties. 

6. The fault-tree approach provides for the analy- 
sis of all elements, subsystems, and systems, and 
includes every factor that influences the hardness 
(or failure) of each element and consequently each 
subsystem and system. It graphically represents 
the logic that relates the failure mode(s) of system 
elements to a particular weapon effect(s). That is, 
the fault tree organizes failure-mode/weapon-effect 
combinations into a logic network, thereby allowing 
use of probabilistic information in descriping the 
hardness of each element, subsystem, and system. 

c The multilevel system organization views the 
protective facility as a hierarchy of systems, sub- 
systems, and elements: 

—Level 1: Complete Faciltiy 

—Level 2: Complete Sytems 

—Level 3: Subsystems within Complete Systems 

—Level 4: Elements within Subsystems 

(1) System Level 1 represents the finished fa- 
cility, complete with protective and protected sys- 
tems. 

(2) System 2 includes complete functional sys- 
tems that make up a particular facility. Included in 
this category are both the protective and the pro- 
tected systems. There are generally nine protective 
systems: 

—Structure 

—Shock isolation 

—Air entrainment 

—Anchors, mounts, and fasteners 

—Penetration 

— Hydraulic-surge protective devices 
— EMP protective devices 
—Radiation shielding 
—Thermal shielding and fire barrier 
The protected systems include: 
—Power supply 
—Power distribution 
—Cooling 

—Heating and ventilation 

—Water supply 

—Sewage disposal 

—Lighting 

— Communication 

—Prime-mission materiel 

—Personnel 
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If required by the facility mission, other protected 
systems may also exist in the facility. 

(3) System Level 3 includes the complete func- 
tional subsystems that make up a system. For ex- 
ample, the various subsystems of the air- 
entrainment system are: 

—Intake structure 
—Expansion chamber 
—Blast attenautor 
—Blast valve 
—Delay line 
—Filtration 

(4) System Level 4 includes the components of 
the various subsystems. For example, the blast 
valve subsystem of the air-entrainment system 
would have components such as: 

—Valve hardware elements 
—Attachment hardware parts 
—Actuation components 
—Sensor elements 

d Typical hardened facilities will include thou- 
sands of structure and equipment items. It is im- 
practical and unnecessary to verify the hardness of 
each item. The analyst must define some consistent 
method for the selection (or screening) of items. 
The specifics of the method may differ from sys- 
tem to system, but a fairly universal approach 
includes these four basic steps: 

(1) Identify those items of equipment that are 
critical mission success. (Note: the noncrit- 
ical items must not pose a threat to the 
critical items when damaged or destroyed 
during attack.) 

(2) Identify those critical items that are suscept- 
ible to the weapon effect. 

(3) Identify those items that have marginal 
hardness. Omit items having a large factor 
of safety. List items that obviously will fail 
and must be redesigned. 

(4) From each population of different items re- 
miaining, randomly select a statistically sig- 
nificant sample for test or perform probabil- 
istic computations to verify hardness. 

2-3 Weapon-effect loads. 

cu After the items have been selected for ver- 
ification, identify local loads induced by one or 
more of the free-field nuclear weapon-effect envi- 
ronments. Identify the total signal path and path 



segments leading from the free field to the ele- 
ment. Determine the transfer function for each 
path segment. In some instances, the inputs to be 
used for hardness verification will be the same as 
those used for design purposes; in other instances, 
the inputs will be determined by extensive testing 
and analytic programs. 

6. Identifying the total transmission path is best 
performed by applying network logic analysis 
(chap. 6). 

c. Determining transfer functions is the process 
of relating free-field nuclear environments to local 
loads in order to assess hardness levels. Transfer 
functions may be determined experimentally; more 
often, an analytical procedure will be employed. 

2-4. Quantifying resistances. 

a. The resistance of a system, subsystem, or ele- 
ment is defined as its ability to withstand nuclear 
weapon-effect loads. The load-resisting characteris- 
tics may be quantified by experimental data 
(laboratory or field) or by analysis of mathematical 
models. Scarcity of experimental data may require 
augmentation from computational models to accrue 
a statistically significant body of data. 

b. Data from either laboratory or field tests are 
preferred for determining resistance. Tests can be 
conducted on prototypes or scale models; but in 
order to be statistically meaningful, a significant 
number of tests must be conducted. (See chap 3 for 
the statistical implications of testing; see chap 4 
for the types of tests available.) If testing is not 
possible because of technical difficulties or lack of 
equipment, time, or funds, a computational analysis 
will be required. 

c. Analytical procedures, especially numeric tech- 
niques, are the backbone of many of the verifica- 
tion studies applied to hardened systems. Although 
deterministic analytical procedures have reached 
an advanced state of development, probabilistic 
analyses that allow for the computation of uncer- 
tainties and survival problbilities are markedly less 
familiar. Thus, the performance of verification 
analyses according to the guidelines in chapter 6 
will require considerable ingenuity. Some guide- 
lines to the development of probabilistic procedures 
that can be applied to real systems are presented 
in chapter 5. 
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CHAPTER 3 
STATISTICAL PROCEDURES 



3-1. Introduction. 

a. Statistical techniques are usually associated 
with laboratory testing. Application of these tech- 
niques to field testing, and especially to the statis- 
tical inferences of calculations, is still in an early 
state of development; consequently, the verificaiton 
anslyst will not have the benefit of extensive pre- 
vious work. 

6. Hardness verification is performed to answer 
either of two basic survivability questions: (1) How 
hard is the system? or (2) Is the system hard to a 
specified level? In answering question (1), the re- 
sistance to pertinent weapon effects is determined. 
This provides a more generally comprehensive and 
useful answer. Then, if the mean resistance is de- 
termined to be greater than the mean specified 
hardness goal, the answer to question (2) is auto- 
matically (and deterministically) "yes"; if it is de- 
termined to be lower, the answer is a no. w In a 
properly designed system, however, a deterministic 
"yes" is not inconsistent with a probabilistic "no" 
and vice versa; there is some probability that a 
system will fail (or survive) regardless of its most 
probable behavior. 

c. Hardness verification generally is not 
conducted to answer Question (1) because it is 
much more costly to determine what the resistance 
is than it is to determine whether the resistance is 
above or below some specified level. If undertaken, 
the result will be a mean resistance with a cor- 
responding variance or uncertainty for each 
weapon-effect failure-mode combination. By com- 
bining the resistance for all t modes of failure, the 
mean-resistance/uncertainty for a particular sys- 
tem or subsystem will be determined. 

d. Verification analysis to answer question (2) 
will result in a hardness statement that presents 
the probability of success (probability of the par- 
ticular failure-mode resistance being equal to or 
greater than the stated hardness goal) with a cor- 
responding level of confidence (e.g., there is a 90 
percent probability that the resistance is >800 psi, 
with a confidence of 70 percent.) 

e. Throughout design and development of the 
hardened system, it is important to acquire . re- 
sponse data to identify failure modes. If possible, 
avoid "go/no-go" tests restricted to a single level in 
order to minimize retesting downstream as ele- 
ments of the prototype design are changed. 



/ The hardness statement must have statistical 
validity regardless of the point of view taken in b 
above. A purely deterministic statement will, not 
be acceptable unless the system, subsystem, or ele- 
ment is superhard or supersoft. 

3-2. Probabilistic distribution. 

a. Ask the following questions: Is it necessary to 
define a specific type of distribution? If so, what 
type of distribution should be selected? 

6. For many (if not all) of the analyses to be 
conducted, the definition of the exact distribution 
will not critically affect the results. (See Aitchison- 
Brown, 1969; Benjamin-Cornell, 1970.) This is true 
for the following circumstances: 
—If the distribution is not grossly nonlinear 

near the mean 
—If the uncertainty parameter Q — is small 
(i.e., Q 2 «1.0) 

where 

a =The true standard deviation 
M —The true mean 

These assumptions are equivalent to ignoring the 
tails of the distribution (i.e., at those responses 
exceeding, say, ±3<r). These assumptions will al- 
most always be acceptable for verification analyses. 
Furthermore, the normal or log normal distribu- 
tion will provide an adequate approximation in 
most cases. Whether the normal or log normal 
model is adopted because of physical considerations 
or as an approximation to another distribution, it 
will be sufficiently accurate in practical applica- 
tions. 

c. If a normal distribution cannot be adopted, 
then either conduct enough tests (usually >30) 
that the assumption of a particular distribution 
can be defined (using, for instance, a chi square 
goodness-of-fit test), or transform the data so that 
the normal assumption may be made, or use a 
distribution-free analysis. Others may be used, but 
only if these are not suitable. 

<L The normal, or t statistic, distribution is the 
single most-used model in applied probability anal- 
yses (Benjamin-Cornell, 1970). The normal distribu- 
tion asymptotically represents the sum of a num- 
ber of random events. For example, it represents 
how a manufactured or constructed item deviates 
from a specified performance value because of 
flaws or errors in the many separate components 
of the item. Therefore, the normal distribution is a 
prime candidate to represent the random 
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deviations of resistance from a specified level, since 
this resistance nearly always results from a sum of 
contributing components. 

e. Adopt the normal distribution whenever as- 
sumptions in b above can be accepted. Use it also 
whenever adequate tests are available and the as- 
sumption of normality is confirmed through a chi- 
square (X 2 ) test. (The performance of the chi- 
square test is described in most of the references 
in statistics included in this volume.) 

/ If the variables are known to result from the 
multiplication of many effects (still assuming 6 
above), use the lognormal distribution instead of 
the normal. However, when Q 2 <<1, there will be 
no significant difference between the normal and 
lognormal results. In this volume only the normal 
(t) distribution is discussed; a complete discussion 
of the lognormal distribution is presented by 
Aitchison-Brown (1969). 

g. Select the nonparametric distribution instead 
of the normal when assumptions in b cannot be 
accepted and another particular distribution cannot 
be confirmed. 

h. Select the binomial distribution only when 
"fail or succeed 0 information is acquired through 
tests at a specified level. 

3-3. Analysis for resistance. 

a. Experiments can be designed to determine the 
resistance of a single weapon-effect/failure-mode 
combination for a simple element, or of combined 
resistance for a total system or complex subsystem. 
In either case, conduct enough independent trials 
or tests to statistically determine the resistance(s) 
of interest, usually >30 samples. 

6. In experimental verification it will not be pos- 
sible to "exactly" identify the resistance of interest. 
The best that can be done is to bracket the resis- 
tance. Ideally, bracketing is done by starting the 
testing at a level where design analysis has shown 
there will be a high probability of success, and 
then increasing the level in small increment* until 
failure occurs. The resistance then lies between the 
level where failure occurred and the level just be- 
fore failure. Obviously, the smaller the incremental 
change in load level, the more accurately the resis- 
tance can be bracketed for the item. Generally, 
specify the resistance as the last level where suc- 
cess occurred. If, however, large increments 
(greater than % of the resistance) must be used 
because of costs, a higher level might be taken. But 
never specify a level more than one-half of the 
increment above the last level of success; and sup- 
port such a resistance level by rigorous analysis. 



c Because testing total systems is expensive, the 
resistance of a total hardened system must usually 
be determined by combining the results of many 
element or subsystem tests. The resistance of an 
element, subsystem, or system is always related to 
one or more weapon-effect/ failure-mode combina- 
tion^). It is critically important to adequately sim- 
ulate the weapon effect of interest and accurately 
measure the resulting repsonses. 

<L The use of statistics in testing for resistance is 
different than statistics used in calculating resis- 
tance. For experimental verification, a number of 
tests are performed on the system, subsystem, or 
element to define its statistical performance. In 
calculated verification, the mathematical model for 
the system, subsystem, or element is deterministic 
but its constitutive parpameters such as stiffness, 
weight, and strength are allowed to have a statisti- 
cal character so that for each selection of param- 
eters chosen a unique response will be obtained. 
The calculations that result will have much the 
same character as test data that result from ex- 
perimental verification. 

3-4. Mean and variance by normal, t, 
distribution. 

cu Given a sample of n data (such as resistance 
measurements), the sample mean value and vari- 
ance are defined as 



n 

1-.L£ x, (3-1) 
i=l 

n 

s 2 = (Xi-5)2 (3 " 2) 

i=l 

or 

n 

i=l 

where 

n —Number of samples 

s 2 ,8* —Sample variances 

x —Sample mean 

Xi —Individual sample values 

Equation 3-3 is presented as a better estimate of 
the variance for n ^30. 
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(Table of t a . v — the 100 a percentage point of the t distribution 
for v degrees of freedom [v = n - 1]) 
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6. The optimum procedure for testing the hy- 
pothesis that the mean of a normal distribution 
has some specified value ji =p 0 ) is based on the t 
test statistic, where: 



t^JEz^iL^n = 2L^L^ (3 . 4) 



c: The acceptance region for a two-sided proce- 
dure (where the hypothesis is stated as — 0 ) i9 



a/2,n - 1 <t<l a/2,n - 1 3-5 



where t is calculated from equation 3-4 with n =m 0 
and a is the level of significance for the test; i.e. 
there is a probability of 1 — a of accepting the 
hypothesis ^ = /i 0 when it is true. The values of 
t&u-i are taken from tables of t statistics such as 
table 3-1. By combining equations 3-4 and 3-5, 



For resistance verification analysis, a should range 
between 0.05 and 0.10. The smaller the value of a, 
the larger the confidence interval. In the extreme 
case it would be possible to have a probability of 
0.9999 (1 — a ) that the true mean falls in an 
interval so large as to be meaningless. This is 
particularly true for small n. 

/ A second probability is called the Type II 
error, which is denoted by 0. The analyst selects P 
to represent the probability of accepting the hy- 
pothesis when it is considered important to detect 
the value of the ratio 0^— At</a, where a= true 
standard deviation of the population and is un- 
known. In other words, /S is the probability of 
accepting the hypothesis Hq: /x = /x 0 when actually 
fi=tii. The analyst must choose values of p x and <r 
such that the ratio (fi x -fi 0 ) is meaningful to the 
analysis. 

g. Once 0 and the ratio m x =At</a have been 
selected, operating characteristics curves such as 
those shown in figures 3-1 and 3-2 can be used to 
determine n. This procedure is illustrated by the 
following example, for which three problems and 
solutions are presented. 



o S 

x-to/2, n-1 <n< x+t^n-i -P= 



8 

The values x±t /2^-rVn" are the confidence 
limits and bound the confidence interval. Equation 
3-6 states that lies within the given confidence 
interval. This statement will be true 1 -a fraction 
of the time. 

d In many instances the analysis to determine 
the resistance will require answering this question: 
Is the resistance equal to or greater than a specific 
value? This allows use of the one-sided procedure 
based on the hypothesis: = <>. In this case, the 
acceptance region is 

t^-t^ 

The hypothesis is accepted if equation 3-7 is satis- 
fied, but is rejected otherwise. This procedure also 
has the probability of 1 - a of accepting the hypoth- 
esis when it is true. 

e. As stated above, 1 -a is the confidence level or 
the probability of accepting the hypothesis when it 
is true. Conversely, the probability of rejecting the 
hypothesis when it is true is given by a. This is 
called a Type I error. The value of a is selected by 
the analyst, using input from the systems manager. 



Example: 

Item Air-entrainment port 

closure 

Weapon Effect Airblast 

Failure Mode Stress 

Resistance Design Goal Pressure wave: 

Rise time (t,) 
=0.002 sec 
Max. pressure (P m ) 
=1000 psi ( o 
=1000 psi, e.g.) 
Decay time (U) 
=0.100 sec (time 
to decay to P m /2) 



Problem 1: The closure was designed to resist 
1000 psi. Can the design be considered 
successful at a level of significance of, 
say, a =0.05? Four resistance values 
have been determined. 

Calculate x (mean value of sample) 
x (variance of sample) 
C.I. (confidence interval) 
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Test Data x x =900 psi 
x 2 =1100 psi 
x 3 =950 psi 
X4 =1000 psi 



- 9UU+11CXH 
nL" Xi 4 



900+1100+950+1000 



i=l 



=987.5 psi 



(using equation 3-1) 



n-1 ^ 



xf (using eq. 3-3) 



i=l 



=1 |(900-987.5) 2 +(1100-987.5) 2 +(950 -987.5) 2 



8=85.4 psi 



x-M 



987.5-1000 
85.4 



X2 =-0.29 

(using eq.3-4) 



<*/2=0.05/2=0.025 
n-l=3 

Prom table 3-1, to.025, 8 =3.182 

Since -3.182 <-0.29 <3.182, the hypothesis that 
M = 1000 is accepted at a =0.05 



C.I.=987.5-3.182X^- to 987.5 +182x-^ 
to 1123 psi. 1 * 



=852 



Therefore, it can be stated that the mean resis- 
tance of the population lies between 852 psi and 
1123 psi, with a confidence of of 95 percent that 
the statement is true. (Had a =0.10 been selected 
for the analysis, then it could be stated that the 
mean resistance of the population lies between 887 
psi and 1088 psi, with a confidence of 90 percent 
that the statement is true; and for a =0.20, the 
CI. =918 to 1057 psi.) 



Problem 2: Can the mean resistance be considered 
equial to or greater than 1000 psi with 
a =0.05? 

One-Sided Procedure: 

The acceptance region is t >— t^-i 



From the example: 

t =-0.29 
a =0.05 
n-1 -3 

Prom table 3-1: toow =2.353 

Since -0.29 > 2.353, the hypothesis is accepted at a 
level of a =0.05. 



C.I. =x-t a< 



mi-i 



to 



Or it can be considered that the C.I. has a lower 
limit only: 



x-t M , -p- =987.5 - 2.353x^p 

Vn 2 

(3-6) Y 



=887 psi 

It can be stated that the mean resistance of the 
population is equal to or greater than 88 psi, with 
a confidence of 95 percent. 



Problem 3: What is the value of 0 for the above 
sample of n =4? 

Two-Sided Procedure: 

Let d= =2.0 (from / above) 

If it is assumed that s is a reasonable estimate of 

a, then 

<r=85.4 psi. 

Since Mo =1000 psi, then Mi =2 X 85.4 +1000 
=1171 psi. 

0 gives the probability of accepting the hypothesis 
that m =1000 when actually m =1171 (or 829). 0 is 
found by using figure 3-1. The intersection of the n 
=4 curve and the d =2 value shows that 0 would 
be approximately 0.244. Therefore, for the sample 
tested there is a probability of 0.244 of accepting 
M=1000 psi, when actually is could be as low as 
829 psi or as high as 1171 psi. 

One-Sided Procedure: 

Let d =^_L1 =2.0 

Again, a =85.4 psi 
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Mi =829 psi 

/S is found by using figure 3-2. 0 is seen to be 
equal to approximately 0.25. Therefore, for the 
sample tested there is a probability of 0.25 of ac- 
cepting /~>1000, when actually it could be as low 
as 829 psi. 



3-5. Nonparametric distributions. 

a. Distribution-free methods are utilized when 
the conditions specified in paragraph 3-26 do not 
apply or if the data set does not readily fit any 
other recognizable distribution. 

6. Use of nonparametric methods requires the 
same kinds of tests as the t statistic described 
above. These tests would result in a measurement o 
the resistance level. The nonparametric results as- 
sert that a particular proportion (p) of the popula- 
tion has failure thresholds falling within a certain 
range. The statement has an associated confidence 
level (P). Procedures for calculation of the mean 
and variance are not defined. 

a Using this method, it can be stated with the 
probability (P) of being correct that no less than a 
particular proportion (p) of the population has fail- 
ure thresholds between the maximum and mini- 
mum value of the failure levels observed in a sam- 
ple of n tests. As before, the probability (P) is the 
level of confidence. The range between the maxi- 
mum and minimum is analogous to the confidence 
interval.* 

<L Assuming that the distribution is continuous 
and that x denotes the proportion of the population 
values that falls within the maximum and mini- 
mum values of any random sample of size n, the 
distribution of x (see Fisz, 1963) is given by 

f(x) =n(n-l)x n - 2 (l-x) (3-8) 

Then 

P =1 -Hp 0 " 1 +(N-l)p n 
An approximation for n from equation 3-8 is 



n 0.24 X* 4 — - +0.5 (3-9) 
1-p 

where 

p =The proportion of the population that will 

fall within the sample range 
Xj4 =The chi-square distribution for 4 degrees 

of freedom 



1 — a =The confidence level (a is the level of 
significance) 

Equation 3-9 is plotted in figure 3-3 for selected 
values of P =1 —a. 



•The interval between any two observations may be used as 
a confidence interval, but here only the sample range is dis- 
cussed, of Sample Size 

e. Analyses can also be performed based on the 
distribution of exceedances. Some significant infer- 
ences can be made from limited sample sizes. Har- 
ris (1952) develops the equation 

P =l-(l-p) n 

where P is the probability that in a future very 
large sample, the proportion p or less of the ob- 
servations will fall below the minimum value * 
observed in a trial sample of size n. (It is also the 
proportion that will exceed the maximum value 
observed in the sample.) Equation 3-10 is plotted in 
figure 3-4 for selected values of n. 



•For the r* smallest value, P - 2 p M-p)"~**) n. 

8— r 



Example: 



Data: 



n 


-6 




= 850 psi 




= 900 psi 




= 1050 psi 


Xmjn 


= 850 psi 




= 1100 psi 


x 2 


= 1100 psi 




= 1000 psi 


*6 


= 950 psi 



Problem: What proportion of the population can 
be expected to have resistance levels 
in the range from 850 to 1100 psi with 
a confidence level of 80 percent? 

The solution can be calculated from equation 3-9 
with 

n =6 

a =0.2 

X 1 4=5.989 (from table 3-2) 

or can be found from the curves of figure 3-3. 
Prom either approach it is found that p =57 per- 
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Figure 3-2. Operating Characteristics Curves for Different 
Values of n 
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Figure 8-1. Operating Characteristics Curves for Different Values of n 
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Figure 3-3. Percentage of Population Within Sample Range as a 
Function 
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cent. Thus, it can be asserted with a confidence of 
80 percent that 57 percent at least of the popula- 
tion will have resistances falling in the range from 
850 psi to 1100 psi. No inferences can be made for 
the remaining 43 percent. (Harris, 1952) 

/. Inferences using distribution of exceedances 
can be made using the same data presented above. 



Problem: With a confidence of 80 percent, what 
is the maximum proportion of a fu- 
ture very large sample which could be 
expected to fall below the 850-psi 
minimum of the sample of 6 observa- 
tions? 

The answer can be calculated from equation 3-10 
above by solving for p with n =6 

P =0.80 

Or the curves in figure 3-4 can be used. From 
either approach, 

p ~ 23.5percent 

Thus, no more than 23.5 percent of a large sample 
from the population can be expected to have resis- 
tances below 850 psi with a confidence of 80 per- 
cent. 



3-6. Binomial distribution. 

a. As stated previously (para. 3-1), the verifica- 
tion of whether the resistance of a system is equal 
to or greater than a specified hardness level re- 
quires a somewhat different approach than deter- 
mining the level of the resistance. This approach 
uses the statistical inferences of a binomial dis- 
tribution. 

b. The data that will be used in the analysis will 
be "failure-success" type of information: Did the 
system, subsystem, or element fail or succeed under 
test or in the calculation when subjected to the 
specified environment? When performing verifica- 
tion testing this does not eliminate the need to 
fully instrument the test item to acquire all per- 
tinent response data. Comprehensive instrumenta- 
tion becomes more important as the specimen and 
test costs increase. 

c. Verification tests or analyses for this type of 
investigation are conducted at a specified hardness 
level, usually the design hardness goal. However, 
tests or calculations can also be conducted at levels 
exceeding the hardness goal. Such a decision is 
made by the verification analyst in conjunction 



with the designer and system manager. If tests or 
calculations are to be performed at levels other 
than the design goal, the level must be related to 
the design goal through analysis that is performed 
separately from the verification analysis presented 
in this volume (i.e., a full quantal response analy- 
sis). 

d. Acquisition of resistance data and the analysis 
for determining whether the resistance is at least 
the specified hardness level will follow basically 
the same sequence of steps as for the resistance 
level itself (para. 3-4). The hardness compliance 
data will almost always be composed of fewer than 
30 samples. However, as previously stated, the sta- 
tistical inferences will be made from a binomial 
distribution because of the "go/no-go" or 
"pass-fail" nature of the data. The basic informa- 
tion required for the verification analysis is pro- 
vided below. 

e. The confidence intervals can be determined for 
two conditions: 



P(s=r|n ' p)= 7!^jT prqa " r (3 " n) 



which gives the probability of realizing exactly s=r 
successes in n trails (where the probability of suc- 
cess for any one trial is p, [q =1— pt p can take 
any value from 0 to 1; 0<p<l. Alternatively, 



n n j 

P(s>r | n,p) = £ s!(n ^ 9)! P' q" (3-12) 
s=r 

which gives the probability of realizing r or more 
(s>r) successes in n trials. 

/. Both equations 3-11 and 3-12 can be evaluated 
directly or tables may be used. Extensive tables are 
included in Abramowitz-Stegun (1972), and for val- 
ues up to n =10 are included for convenience in 
appendix A. 

g. Generally, equation 3-12 will be used in the 
verification analysis, since the question usually 
posed to the verification analyst is: "Is the resis- 
tance equal to or greater than r, where r — r/n?", 
rather than, "Is the resistance exactly equal to r?" 

h. In the population being considered (the total 
number of like elements, subsystems, or systems in 
the total deployed system), an unknown proportion 
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Table S-S. The x* Distribution 
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Figure 3-i. Percentage as a Function of Sample Size and Probability 
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of items will withstand the imposed environment 
(design goal environment). This proportion will al- 
ways remain unknown unless the entire population 
is tested. However, this unknown proportion (p) 
can be bounded by confidence limits for a specified 
confidence coefficient (P) (e.g., 90 percent such that 
p^<p <p M or p >p^in which p^and p M are the lower 
and 90 upper bounds, respectively. 

i. The discussion in paragraphs 3-4 e and / re- 
garding a and 0 are also applicable to binomial 
distributions. The basic approach (Mace, 1964) is 
summarized below. 

j. The hypothesis to be tested is p >p 0 , where 
p is the true but unknown proportion of successes 
in the population and p 0 is the least favorable value 
that is acceptable. As before, a is the probability 
of rejecting H, when actually it is true, and 0 is 
the probability of accepting p>p 0 when actually p 
=Pi, where px <p 0 



Let: 



n = Number of tests in sample 
s = Number of successful tests 



Define: 



en 



1-P. 
1-Pi 



m = 



P. 1-Pi 



(3-13) 



0 



1- a 



Pl_ + 
P. + 



1-P. 
1-Pl 



(3-14) 



^ 



1-/8 



(3-15) 



1-Pi 



Then 

Accept H„: p >p 0 if s >h 0 +nm. 
Reject H 0 : p >;p 0 if s <hj + nm. 

Conduct an additional test if 
hj +nm =<s <ho =nm 



Example- 



Problem: 



Test the hypothesis: H„:p <(p 0 
for a =0.10. 



=0.90) 



Let p,= 0.70 and /3=0.l 



(Note. /8 =0.30 implies that there is not much 
concern about accepting p>0.90 when actually p 
=0.70. If there is concern, choose fi <0.05.) 

Data: 
n =5 
s =4 



0.1 
0.3 



m= 



' n 09 + '"61" 



m=- 



-1.0986 



-1.0986 



-0.2513+(-1.0986) -1.3497 



=0.814 



h„=- 



-1.3499 
0.7 



h,= PA. = it><tOT = -1.442 

-1.3499 -1.3499 

\ +nm =-1.442 +(5 x 0.814) =2.628 



-1.0986 
-1.3499 

1.9459 



=0.814 



ho +nm =0.814 +(5x 0.814) 
2.68 <s =4 <4.884 



=4.884 



The data show that h, =nm <s ho +nm. There- 
fore, an additional test must be conducted to verify 
the hypothesis at the a and p probabilities se- 
lected. 

For n =9 

s =8 

h„ +nm =0.814 +(9 x 0.814) =8.140 
s <ho +nnv 
s <hj +nm 
For n =10 

s =9 

h, +nm =0.814 +(10 x 0.814) =8.954 
s<ho +nm 

Accept H o :p>0.90 

Thus it would take 10 tests with 9 successes to be 

Google 
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able to assert that p >0.90, with a confidence of 90 
percent, that the hypothesis will be accepted when 
it is true, and with a probability of 0.3 of accepting 
the hypothesis when actually p could be as low as 
0.70. 



k. Figure 3-5 shows plots of the functions (h/n) 
+m and (h^n) +m for values of n for p 0 -0.90, 
0.10, and p —0.25. The curves are developed 
from equations 3-10 and through 3-14 and can be 
used as follows. 

2. Select a value of px that is desired for the 
analysis and calculate Pi/p* From the (h</n) +m 
set of curves, find the value of s/n required for 
acceptance of Ho:p>Po for the desired n. From the 
(h/n) +m set, find the value of s/n that would 
require rejection of Ho for the desired n. 



Example: 

Let Pi=0.70 

Px/po-0.78 (P 0 =0.90) 



Accept Ho:po^0.90 
when 



s/n> 



Reject Hoipo>0.90 
when 



s/n< 



10 0.91(10 10) 

7 0.94 (7/7) 

5 1.00 (5/5) 

4 1.05- 

3 1.12- 



10 0.66 (6/10) 

7 0.60 (4/7) 

5 0.52 (2/5) 

4 0.44(1/4) 

3 0.32 (0/3) 



The above data show the following: 

- For p 0 =0.90, px -0.70, <* =0.10 and 0 
—0.25, it is impossible to verify the 
hypothesis with fewer than 5 tests, i.e., 
it is impossible to have s>n. 
For a sample of 5, all tests must be successful 

to accept H^p>0.90 
As tests are conducted, the criteria for 

rejection can be constantly monitored; 
for example, if after 4 tests, 3 have 
failed (s =1, s/n — W), the hypothesis 
must be rejected for the parameters 
selected above. 

3-14 



Problem: What is the proportion (p) of systems 
that will survive the design-goal test 
environment with an associated con- 
fidence coefficient of 90 percent? Thus 
we must find p such that p >P with 
a confidence of 90 percent. 

Data: n =3 (tests conducted at the design- 

goal level) 

r =2 (minimum number of successful 
tests) 

From table A-2, column (n =3, r =2), find the 
value nearest to 0.10 (0.10 -1.0 - 0.90). The value is 
0.104. In the p column read the corresponding 
value, which is 0.20. (Note: the table values can be 
interpolated if desired to get the value of p for 
exactly 0.10.) This value (0.20) is p . The confidence 
interval is p >0.20 with an associated e confidence 
coefficient of 0.90. The statement can now be made 
(based on the test data) that the proportion of 
successes to be expected in the population is 
greater than or equal to 0.20 and there is a 90 
percent confidence that this is true. 

m. When the system manager assumes a priori, 
what constitutes an acceptable confidence interval 
and when the computed interval seems too low, two 
alternatives exist Conduct more tests with the ex- 
pectation that some of those tests will be success- 
ful, or decrease the expectation of the confidence 
limit and accept the greater risk. For example, if 
two or more tests are conducted and both are 
successful, then 

n =5 

r =4 

From table A-2, the 90 percent confidence limit is: 
p>0.42 

If no additional tests could be conducted, a lesser 
confidence must be accepted to decrease the inter- 
val or increase p. For 

n -3 

r -2 

and a 70 percent confidence, the confidence limit is: 
p>0.37 

For 50 percent confidence, the confidence limit is: 
p>0.50 

The results of this example support the observa- 
tions of the previous section, i.e., it is impossible to 
make high confidence statements with a meaning- 
fully narrow confidence interval without conducting 
a significant number of tests (even if all tests are 
successful). This is further demonstrated by the 
curves in figure 3-6. 
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Figure 8-5 Determination of Minimum Sample Size 
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CHAPTER 4 

VERIFICATION REQUIREMENTS AND EXPERIMENTAL METHODS 



4-1. Introduction. 

a. In verifying whether systems will survive a 
postulated nuclear attack, both experimental and 
analytic (mathematical) methods may be required. 
The selection of the methods to be implemented is 
the responsibility of the systems manager and the 
verification analyst. The decision is based on sys- 
tem design and mission requirements as well as on: 

—Accuracy of mathematical models 

—Cost of calculations 

—Accuracy of experimental simulation of 
environments 

—Cost of experiment 

—Number of items available for analysis 

6. Although protective and protected systems 
must be capable of withstanding the 13 weapon 
effects described in TM 5-858-2. each system will be 
more susceptible to some effects than to others. If 
the facility and its contents are resistant to par- 
ticular effects, the system is defined as super-hard 
with regard to those effects, and verification will 
not be required. Conversely, some systems may be 
very susceptible to some effects, and redesign 
rather than verification is required. It is primarily 
toward those systems whose survivability is mar- 
ginal that the verification program is directed. 

c. Verification studies must consider the free- 
Held weapon-effect load, the local load at the sys- 
tem, subsystem, or element level via the transmis- 
sion path between the free field and the local load 
point, and the resistance of the system in question. 
System managers must define the free-field 
weapon-effect load, the transmission path charac- 
teristics, and the system resistance probabilistical- 
ly, i.e., in terms of their mean values and their 
coefficients of variation (COV's). In general, a de- 
terministic assessment of hardness will not be ac- 
ceptable. 

4-2. Simulation requirements. 

a. Definition of the local environments that ex- 
cite systems, subsystems, and elements is a signifi- 
cant effort that must be accomplished to provide 
input to the verification analysis. Although it may 
not always be the case, the verification analysts 
likely will inherit responsibility for defining local 
loads, since they have access to the data from 
which input loads would be derived and they are in 
the best position to understand the true nature of 
those environments. 



6. The transmission of the primary weapon-effect 
loads from the free field to the input points of the 
systems is highly dependent on the threat (weapon 
size and number, height or depth of burst), the 
siting conditions (geology), and the configuration of 
the facility. The general characteristic of the free- 
field nuclear weapon effects are presented in TM 
5-858-2. Although it is not practical to present a 
comprehensive discussion of the methods for simu- 
lating the transmission of the particular weapon 
effect from the burst point to the hardened target, 
a recommended methodology is presented in chap- 
ter 6. 



c. The system/weapon-effect combinations of 
protective system for which hardness assessment 
must be accomplished are presented in table 4-1. It 
must be recognized that each of these systems will 
contain numerous subsystems and elements. In 
many instances it will be necessary to perform the 
verification analysis at the element or subsystem 
level. For each separate item to be analyzed, spe- 
cific requirements and methods will need to be 
defined. Simulation requirements and techniques 
must be fully developed, evolving finally into a 
comprehensive test and analysis plan to guide the 
conduct of the analysis. 



<L Protected systems must also be subjected to 
hardness verification. These systems are not sepa- 
rately addressed in this manual. However, in most 
cases they must be hardened only to the local 
environment within the protective structure or on 
the shock-isolation system. 



e. Experimental techniques are used to generate 
input that simulates the overpressure, ground 
shock, and EMP nuclear weapon effects for ver- 
ification analyses. Exclusion here of techniques for 
the other effects presented in TM 5-858-2 does not 
imply that these effects can be ignored. The sus- 
ceptibility of each subsystem to all effects must be 
determined and verification testing conducted if 
necessary. However, in general the subsystem de- 
sign^) will be governed by these three weapon 
effects. The GE-TEMPO (1972) and Bednar (1968) 
reports present detailed information for some of 
the techniques that were in use in 1972 and 1967, 
respectively. Excerpts from both are presented in 
appendix B. 
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4-3. Simulating and testing 
overpressure. 

Design. The airblast environments will have ex- 
treme variations. Surface/flush structures/ systems 
or the surface elements of air-entrainment systems 
may be designed to survive peak overpressures of 
1000 psi and greater, whereas aboveground facili- 
ties may only be required to survive tens of psi or 
less. No single technique or tool is adequate to 
providing the simulation-testing capabilities 
required for verification throughout this range. 

(1) Because large field tests usually are pro- 
vided by a central agency, involve numerous other 
agencies, and may be available only once during 
the total verification process for a particular fa- 
cility, the opportunity should be optimized. This 
means that the resistance of all critical elements of 
the system or subsystem to be tested must be 
estimated beforehand with corresponding 
uncertainties. Furthermore, the incorporation of 
these resistances into a complete system, as de- 
scribed in chapter 6, must show a high probability 
that the system will survive. The system is then 
fielded such that the design environments are sim- 
ulated as closely as possible. Replication of the 
tested system at different ranges from ground zero 
is very desirable in order to acquire data above and 
below the design environment. 

(2) It is most important that the actual local 
load inputs and the failure-mode responses be ac- 
curately measured to provide results of high statis- 
tical confidence that will credibly support a final 
post-test verification statement. Depending on the 
degree of simulation achieved for the environments 
of interest, it may be necessary (as discussed 
above) to base the final verification statement on a 
calculation from a verified mathematical model. 

6. High-exjkjsive simulation technique (HEST) 
tests. The basic configuration of a HEST test is 
illustrated in figure 4-1. A cavity is constructed 
over the test item. The amplitude of the peak over- 
pressure is controlled by the density of explosives, 
generally Primacord, in the cavity. The positive- 
phase duration is governed by the depth of the 
surcharge, and the propagation velocity is governed 
by the wrap angle of the Primacord. 

(1) Combustion products of the Primacord are 
released at a rate that is proportional to the 
Primacord wrap angle. These products act as a 
piston to drive a shock wave into the undisturbed 
cavity volume. As the generated pressure acts 
against the overburden, the cavity volume increases 
as the overburden rises, causing the prssure am- 



plitude in the cavity to decay. It is this process 
that produces the fast rise and exponential decay 
to simulate the nuclear overpressure. 

(2) A HEST test can be configured to produce 
a blast wave that provides, for short durations, 
good simulation of the overpressure from weapons 
that extend up to the megaton range. Because of 
the limited duration that can be simulated, stress 
attenuates in the ground more rapidly with depth 
than for the nuclear event and the displacements 
will be smaller. Hence, the early dynamic responses 
of the tested system will approach those that 
would be induced by the nuclear event, and good 
simulation of peak accelerations and velocities is 
achieved; but late-time rigid body motions are not 
well simulated. 

(3) The technique is applicable primarily for 
testing surface-flush and shallow-buried structures 
where the principal failure mode considered is di- 
rectly related to overpressure loading. Dynamic 
pressure effects are not well simulated by the 
HEST test, and failure modes related to this effect 
must be tested separately. Because the peak over- 
pressure is essentially the same over the entire test 
area, specimens of large dimensions can be realisti- 
cally loaded in a HEST test. Both full-scale and 
model structures can be tested in the HEST test 
bed. 

(4) Even though the positive-phase duration 
does not simulate the nuclear event at late times, 
the duration is usually sufficient to provide ade- 
quate loads for direct verification when the over- 
pressure impinges directly on the test item. Thus, 
good verification data can be obtained for most 
deformational failure modes for surface-flush sys- 
tems such as structures and closures. Further in- 
formation on HEST is available in D'Arcy et al. 
(1965), D'Arcy-Clark (1966), GE-TEMPO (1972), 
and ble Traindafilidis-Zwoyer (1968). 

c. High explosive (HE) field tests. Many HE tests 
have been conducted to provide overpressure and 
ground-motion environments in sizes varying from 
a few hundred pounds to hundreds of tons of equiv- 
alent TNT. For overpressure studies the explosives 
are placed in spherical or hemispherical arrange- 
ments at or above the -earth surface, typically as 
shown in figure 4-2. 

(1) The blast wave produced in an HE test 
provides good simulation of the peak amplitude for 
both the overpressure and dynamic pressure. The 
peak pressure amplitudes from an HE burst decay 
more rapidly with distance from the center of the 
charge than do nuclear bursts because usually less 
weapon yield is available. Because of this, only the 
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Figure U-l. Sketch of HEST Facility 
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larger tests (e.g., 500 tons) are adequate for testing 
full-scale elements. For large systems, the system- 
size/peak-pressure combination of interest will 
generally dictate that testing be performed on 
models of subsystems. Certain elements such as 
closures and blast valves may be subjected to good 
simulation of peak loads for direct verification. 

(2) HE tests provide an airblast environment 
that is comparable" to the nuclear environment es- 
pecially at peak overpressures below 300 psi. Good 
simulation of overpressure is provided for a wide 
range of yields but particular attention must be 
directed to proper scaling. The yield dependence is 
discussed in chapter 7 of TM-5-858-2. Further in- 
formation on HE tests is available in Kingery 
(1968) and in Reisler et al. (1975). 

d. Airblast simulators. In addition to the large- 
scale field tests described above, there are fixed 
facilities designed to produce blast waves that sim- 
ulate to some degree the blast wave from a nuclear 
event. Many of these facilities are listed in appen- 
dix B and are useful for component or element 
testing. Two of the larger facilities of this type are 
GRABS and LBLG, described below. 

e. Giant reusable airblast simulators (GRABS). 
The GRABS facility is located at the Kirkland AFB 
and is operated by the Air Force Weapons Labora- 
tory (AFWL). The basic configuration of the 
GRABS facility is shown in figure 4-3. The facility 
provides for a soil depth of up to 30 feet for 
emplacement of the specimen to be tested. The 
amount of explosives, cavity volume, and soil sur- 
charge depth are selected to produce the pulse de- 
sired within the range of the facility capabilities. 
The mechanism of producing the pulse is similar to 
that of the HEST test except that the pressure 
pulse does not horizontally traverse the specimen 
but loads are areas of the test bed simultaneously 
and essentially equally. Simulation of the peak 
overpressure history is similar to that achieved in 
a HEST test. The size of the facility allows for 
testing of full-scale elements; subsystems and sys- 
tems would generally be modeled for testing. The 
test s design must consider the reflecting bound- 
aries of the facility. The s response time(s) of the 
failure mode being studied must be such that the 
critical parameters have reached peak amplitudes 
before reflected waves reach the location of inter- 
est. Further information on GRABS is available in 
Bartel and Jackson (1973) and Jackson et al. (1973). 

/ Large blast load generator (LBLG). A number 
of blast load generators listed in appendix B and 
discussed by GE-TEMPO (1972) and Bednar (1968), 
of varying sizes and designs have been used to 



conduct verification tests on full-scale elements and 
models of elements and subsystems. One of the 
largest of such facilities is the Large Blast Load 
Generator (LBLG) located at the U.S. Army En- 
gineer Waterways Experiment Station (WES). The 
basic configuration of the LBLG is shown in figure 
4-4. Further information on LBLGs is available in 
Albritton (1965) and in Kennedy et al (1966). 

(1) A soil bed of up to 10 ft deep can be loaded 
for emplacement of specimens to simulate buried 
conditions. Peak pressure and positive phase dura- 
tions can be controlled within the range of the 
facility capabilities. Durations equal to or exceed- 
ing those from megaton nuclear events can be 
achieved. 

(2) Airblast loads for elements and subsystem 
models can be provided to accurately simulate nu- 
clear effects. As with the GRABS facility, the test 
must be designed such that reflections from the 
facility boundaries will not affect the failure mode 
response being studied until maximum amplitudes 
have been reached. 

g. Shock tribes. A number of shock-tube facilities 
are listed in appendix B and discussed by GE- 
TEMPO (1972) and Bednar (1968). In general, such 
facilities provide a good simulation of the airblast 
characteristics of a nuclear event, namely peak 
overpressure, decay rate, positive-phase duration, 
and dynamic pressure. Shock tubes are used for 
testing subsystem models and small full-scale sys- 
tem elements. However, depending on the response 
of the failure mode being studied, verification can 
be accomplished through scaling. Direct full-scale 
verification testing of the performance of air- 
breathing equipment can also be accomplished in 
shock tubes. Shock tubes have also been used to 
perform debris-impact tests. Further information 
on shock tubes is available in Lane (1971). 

4- 4. Simulating and testing ground 
shock. 

a. Overpressure-induced and crater-induced. The 
techniques and facilities that are used to simulate 
the ground-shock environment from a nuclear 
event are the same whether the shock is induced 
from the overpressure or the crater formation. The 
environments, motion and stress, created by the 
two sources are very different as presented in TM 

5- 858-2, and are dependent on combinations of 
weapon nt size, height or depth of burst, and sur- 
rounding media. Some of the tests t described 
above for simulation of overpressure are also used 
for ground shock. It is generally the objective, par- 
ticularly when testing full-scale subsystems or ele- 
ments, to simulate the airblast and ground shock 
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Figure 1>-S. The Giant Reusable Airblast Simulator (GRABS) 




U.S. Army Corps of Engineers 
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Figure U-2. Typical HE Field Test Configuration (Spherical 
Charge) 
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simultaneously. Only in rare circumstances is this 
actually accomplished. Usually one environment or 
the other must be compromised. Whether over- 
pressure or ground shock is most important is a 
decision that must be made by the verification 
analyst. The decision is based on the susceptibilit 
of the failure mode(s) being studied to the particu- 
lar design environments. Many times, separate 
tests will be required to achieve complete verifica- 
tion of the system. 

b. HEST test The HEST test described above for 
simulation of overpressure is one of the best tech- 
niques for studying overpressure-induced vertical 
ground motion over relatively large areas. How- 
ever, as with all simulation techniques, a com- 
pletely satisfactory simulation cannot be achieved. 

(1) Because of the complex geometry of the 
test facility and the relatively small surface on 
which the pressure load acts, the ground motions 
that occur in a HEST test probably are more com- 
plex than comparable responses in a nuclear envi- 
ronment. The HEST test does not simulate all as- 
pects of nuclear ground shock; because of the finite 
dimensions of the test cavity "edge effects" distort 
the measured response. 

(2) Analyses performed for HEST programs 
indicate that, in some regions beneath the test bed, 
peak vertical and horizontal-longitudinal (in the 
direction of the blast-wave propagation) accelera- 
tions, velocities, and stresses are reached before 
edge-induced (relief wave) responses can influence 
overpressure-induced responses. Figure 4-5 shows 
these regions. The analyses also show that the peak 
displacements, which occur late in time and include 
edge effects, do not simulate nuclear responses. 
Other responses, such as the horizontal-transverse 
(perpendicular to the direction of the blast-wave 
propagation) and beneath-berm motions, are also 
HEST-unique and have no direct relationship with 
the nuclear problem. 

c. High explosive (HE) field tests. HE field tests 
designed to simulate overpressure also produce 
ground shock. However, the limited yield of the HE 
device usually precludes testing of items other than 
elements or models. Moreover, ground motion is 
critically dependent on the geology, so it is difficult 
to locate an HE test site that scales geologically to 
an operational site. 

(1) The predominance of one motion over the 
other (overpressure-induced and crater-induced) de- 
pends on the different attenuation rates and con- 
sequently on the distance from ground zero where 
motion is measured, and on the particular motion 
parameter considered and the associated direction. 



Overpressure may dominate the acceleration at all 
ranges of interest to this volume except deep-based 
systems. However, the crater-induced accelerations 
may cause larger maximum displacements. In gen- 
eral, overpressure-induced motions will contain 
higher frequency components than the crater- 
induced motions, especially near the ground sur- 
face. 

(2) The familiar HE test configuration shown 
in figure 4-2 is not effective in simulating the 
crater-induced ground-motion weapon effects. To 
circumvent this, special configurations of charge 
burial are used as shown in figure 4-6. Using finite 
element computer codes, the configuration of the 
charge is calculated for the particular geology, 
yield, and HOB/DOB combination of the nuclear 
weapon to be simulated. This technique produces 
good simulation of the crater-induced motions but 
does not simulate the overpressure. At this time 
(1976), accurate verification of both effects cannot 
be achieved in a single test event. Further informa- 
tion on HE tests is available in GE-TEMPO (1967, 
1970, 1973), Stubbs et al. (1974), and Rooke et al. 
(1974 and 1976). 

<L HEST-BLEST. The HEST test is limited in 
effectiveness by the size of the test bed/cavity; the 
HEST-BLEST concept shown in figure 4-7 was de- 
veloped to overcome this deficiency. Here, the 
HEST cavity is constructed over a relatively small 
area where it is necessary to produce the overpres- 
sure directly on the surface of the structure sys- 
tem. Separate shallow-buried, high-explosive 
charges are then placed in a much larger surround- 
ing area. This technique of loading the area is 
called a Berm Loaded Explosive Simulation Tech- 
nique (BLEST). 

(1) The BLEST charge array is designed to 
simulate the stress environment in the upper sur- 
face of the test area that would result from nuclear 
overpressure loading. This loading combined with 
the direct air-bias from the HEST cavity simulated 
the overpressure from a nuclear s event over a 
much larger area and for much larger times than 
in a simple HEST test. The peak displacements are 
not influenced by the unloading waves from the 
edge of the HEST cavity. 

(2) This test does not simulate crater-induced 
ground shock. Further information on HEST- 
BLEST is available in Schrader et al. (1976). 

e. HEST-DIHEST. As evidenced by the above 
discussions on HEST and HE tests, there is no 
single technique that produces a totally satisfac- 
tory simulation of both the overpressure-induced 
ground shock and the crater-induced ground shock 
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Figure 4-6. HE Test Configuration to Simulate Crater-Induced 
Ground Shock 
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Figure i-7. HEST/BLEST T Configuration 
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effects. Therfore, tests have been designed that 
combine the HEST test with DIHEST 
(Direct-Induced High-explosive Simulation 
Technique) to produce both effects simultaneously. 

(1) The basic configuration of a typical HEST- 
DIHEST test is shown in figure 4-8. This test is 
designed to combine overpressure and limited 
ground shock from the HEST test with crater- 
induced ground shock from the DIHEST test. The 
DIHEST test uses a buried array, usually planar, 
of high explosives that are detonated 
simultaneously, to produce a predominately hori- 
zontal shock wave on a test system. Timing of the 
detonation of HEST and DIHEST explosives sys- 
tems and location of the DIHEST array require 
extensive pretest experiments and calculations. 

(2) Tests have been conducted in both rock and 
soil media. The size of the HEST cavity, as well as 
the design of the DIHEST array, influences the 
duration over which the ground motion is simu- 
lated. Generally acceptable simulation of accelera- 
tion, velocity, and stress peak amplitudes are ac- 
heived. 

(3) Further information on HEST-DIHEST is 
available in Schlater (1972) and Blouin (1969). 

/ Mechanical tests. The tests discussed for the 
simulation of overpressure and ground shock apply 
to scale model and full-scale testing of systems and 
subsystems. However, before field or large-scale 
laboratory tests are conducted, and before a system 
has been fabricated, comprehensive testing will be 
performed on elements and small subsystems. 
These tests, conducted to obtain data for basic 
development and verification, are mechanical tests, 
in which the excitation force is applied directly to 
the specimen through a known mechanical connec- 
tion. 

(1) Mechanical tests are generally the simplest 
and least expensive to conduct; the specimens are 
relatively small (i.e., measured in feet rather than 
tens or hundreds of feet), and the test equipment is 
readily available. Because of these factors, this 
type of testing is popular, and most of the statisti- 
cal data to determine the mean resistance and 
uncertainty associated with a particular failure 
mode will be acquired with these tests. 

(2) The greatest limitation of the mechanical 
test is that it is not directly comparable to the 
free-field weapon-effect ground shock. Therefore, it 
is necessary that supporting analyses and tests 
(often very complex and costly) be conducted to 
translate the free-field ground shock to the local 
input environment. Nevertheless, compared to the 



unwieldy field test and the large-scale laboratory 
test, the mechanical test will generally be the most 
used tool in the verification process. 

(3) There are numerous techniques and facili- 
ties used to conduct mechanical tests. A discussion 
of some of the most applicable facilities is pre- 
sented in appendix B and in GE-TEMPO (1972) and 
Bednar (1968). However, in many instances it will 
be necessary to apply basic force-producing equip- 
ment such as shakers or rams, in a manner that is 
unique to the element-input requirement combina- 
tion being addressed, and no existing fixed facility 
would be adequate. 

(4) In many instances it will be required to 
perform a verification test of a large subsystem in 
a full-scale operational mode (e.g., of critical con- 
trol and communication equipment supported on a 
shock-isolated floor or platform). In such cases, 
where fixed-facility tests are prohibitively costly, 
pulse-train simulation tests have been conducted 
for verification. This technique applies specifically 
designed force-pulse trains at the required number 
of locations, directly on the platform structure, 
such that the platform response simulates the re- 
sponse experienced by the total shock-isolated sys- 
tem as a result of the motion induced in the pro- 
tective structure by the ground shock. The tech- 
nique uses measurements of impedence and transfe 
functions on the as-built, in-place subsystem to 
transfer the motion from the attachment point of 
the protective structure to critical locations on the 
platform. With the equipment elements in place, 
these motions are then simulated on the platform 
by application of the force-pulse trains. 

(5) Further information on mechanical tests is 
available in Safford and Walker (1975a and b). 

4-5. Simulating and testing EMP. 

a. Uniqueness. The waveform and corresponding 
frequency spectrum for a nuclear EMP differ sig- 
nificantly from other man-made electromagnetic 
disturbances or from natural phenomena such as 
lightning. The pulse rises very rapidly 
(nanoseconds) to the maximum level (hundreds of 
kilovolts per meter for the electric field) and then 
decays exponentially. The frequency range is exten- 
sive, varying from UHF to VLF. Also, the EMP 
field is widely distributed, whereas g lightning is 
localized. Because of these characteristics and dif- 
ferences, designing to protect against EMP or to 
simulate for verifying hardness requires EMP- 
unique approaches. The verification of EMP hard- 
ness will us both analysis and testing techniques. 
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6. Verification analysis. The analysis techniques 
presently available are sophisticated enough to 
identify areas of weakness in the system design, to 
guide test selection/design, and to confirm the ac- 
curacy of test data. However, specific points of 
weakness and specific quantitative levels of hard- 
ness cannot be calculated. 

c Experimental analysis. The analysis of the 
system can be accomplished using computer codes 
that model the many elements of the system. An 
experimental analysis can also be performed that 
consists primarily of measuring the impendance 
and transfer functions between critical system 
junctions and then involving these functions in 
time, or multiplying them in the frequency domain 
with the postulated threat environment inputs. The 
predicted total response will be very accurate as 
long as the responses of all elements in the mea- 
sured path remain in the linear range. 

d. Verification testing. Both laboratory and field 
tests are conducted to provide the required simula- 
tion for hardness assessment. Hardness evaluation 
for EMP may also use scale-model and component 
testing; however, these tests are generally not used 
for final verification analysis. In general, all sys- 
tem components will be separately tested in the 
laboratory prior to being included in a system or 
subsystem test, the test data being integrated. 
Testing is used to verify calculations to confirm 
general hardness of the design and identify specific 
weaknesses, and to provide bounds for final ver- 
ification testing. 

c. Subthreat-level testing. Subthreat-level tests 
performed on elements or on systems to define the 
manner in which the EMP field couples to the test 
item and to expose weaknesses that were not dis- 
covered in the analysis. Techniques used for 
subthreat-level testing are the following: 

(1) Low-level transient: The threat time his- 
tory is simulated but the amplitudes are below the 
threat level. 

(2) Repetitive pulse: A train of pulses (10 to 
100 per second) is applied. The threat time history 
may not be simulated, but the frequency range 
must be the same as that of the threat pulse. 



(3) Continuous wave (CW): A continuous wave 
is applied to the system. The wave may be swept 
across a broad frequency range or applied at dis- 
crete frequency steps. 

(4) Direct injection: Any of the above signals 
can be directly applied at one or more parts of the 
systems. 

(5) Stationary field: Small items are subjected 
to field-coupling tests. 

/ Threat level. Threat-level testing requires that 
the simulation of the threat be applied over a 
volume that is large relative to the total system 
being analyzed. When this is accomplished, the 
hardness (or lack of hardness) is verified by ob- 
serving selected critical system responses. Since the 
verification is accomplished through observation of 
system responses, the threat-level test environment 
must accurately simulate at least the following: 
—Propagation direction of both the electric 

and the magnetic fields 
—The pulse shape and frequency spectra of 
both fields 

—Peak amplitudes greater than the threat 
magnitude 

—Relative magnitudes of electric and magnetic 
fields 

(1) The above requirements must be satisfied 
for each threat scenario for which hardness is to 
be verified, and the electric and magnetic fields 
must be developed over a volume larger than the 
volume of the system to be tested. The arrival 
direction and polarization must be addressed if 
total hardness is to be verified. It is also necessary 
that the on-line operation of the system be ac- 
curately defined such that the simulation can be 
imposed at the time(s) of maximum susceptibility. 

(2) High-level transient field tests are required 
for total threat-level testing. However, some degree 
of verification can be accomplished at threat level 
using direct injection. 

g. Simulation facilities. There are numerous fa- 
cilities that have been developed to provide EMP 
simulation for various systems and components. 
Many of these facilities are listed in appendix B. A 
comprehensive treatment of EMP design and ver- 
ification can be found in COE (1974) IITRI (1973), 
Whitson (1973), and Schlegel, et al. (1972). 
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CHAPTER 5 
ANALYTIC TECHNIQUES 



5-1. Introduction 

a. The analytical approach is an important and 
well-established tool to support hardness verifica- 
tion. In situations where a physical system does 
not yet exist or, if it does, cannot be tested, the 
analytical approach may be the only way to assess 
hardness, survivability, vulnerability, etc. Unlike 
experimental verification which naturally lends it- 
self to a statistical presentation, computational 
verification almost always depends on deterministic 
tools to provide system response to nuclear weapon 
effects. These deterministic tools generally produce 
sophisticated information of the most probable re- 
sponse of the system (the mean), but quantifying 
uncertainty will involve a laborious selection of 
values for the random nature of system param- 
eters. 

b. Although the design load for a system, sub- 
system, or element will sometimes be the input for 
verification analyses or tests, it may also occur 
that updated information will be available to better 
define the local load environment. Thus, the ver- 
ification analyst may not only be responsible for 
assessing the survivability of the system to its 
design loads, but he may also be required to re- 
assess the local design load to reflect more realistic 
conditions. Part of the reassessment will come 
about as a natural consequence of verification stud- 
ies because the response of one system may be the 
load to another. Verification of the first system 
may indicate that the load to subsidiary system is 
different than expected. If the verification analyst 
is also charged with the responsibility for deter- 
mining system survivability, he must be prepared 
to view the verification study from a system en- 
gineering point of view, since he will be the one 
most likely to intelligently assess the adequacy of 
previously established design criteria. 

5-2. Computing uncertainties and 
survivability. 

a. Uncertaintie&The total coefficient of variation 
of a random variable X, designated by ft x is ob- 
tained from 




Here, 5 X models the natural randomness of X, 
whereas A x represents the uncertainty arising from 
errors in estimation. If X itself is a function of 
several random variables, 

X = f ( Y lf Y* Y*..., Yn,..., Y N ) (5-2) 



then using first-order linear approximations 
(Ang-Cornell, 1974) 

x =f(yi,y2,y3,.»,yn,...,yN) (5-3) 

and 

n=l x 70 n^m v ' 

where 

ft f =Total coefficient of variation associated 
with the functional form of X 

fly = Total coefficient of variation associated 
' withY„ 

- y*».+Ai. 

Pnm= Correlation coefficient of Y n and 
Y^-l^p^l 

The subscript w o" denotes that 6f/6Y n is to be 
evaluated at the mean values of the variables. Im- 
plicit in equations 5-3 and 5-4 are the assumptions 
that ftY n 2<<1 and nonlinearity in f near the 
mean values is not large. 

6. Survival probability. Consult chapter 9 of TM 
5-858-1 for a description of survivability models 
and chapter 11 therein for a presentation on sur- 
vivability allocation. 

5-3. Types of analytic techniques. 

a. In performing verification calculations the out- 
standing problem stems not from finding a suitable 
deterministic model, but rather— once having se- 
lected a model— from interpreting results pro- 
babilistically. There have been attempts to produce 
compendia of the many analytical techniques that 
have been applied to calculating the response of 
systems undergoing nuclear attack (AE WES, 
1972). There is almost never an obviously "best" 
technique to solve a particular problem. The "best" 
solution depends on understanding all available 
techniques; the time frame and cost constraints; 
the accuracy desired, i.e., whether the assessment 
is preliminary of final; the relative elegance of 
concomitant analyses; and the state of the art. 
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6. To be meaningful, verification analyses must 
be performed probabilistically, i.e., local loads and 
system responses must be expressed in terms of 
uncertainties, and system survivability must be ex- 
pressed as a probability. Uncertainty is defined as 
the standard deviation of the load and response 
normalized to the mean. It usually occurs that the 
the best estimate of the calculated mean response 
is obtained with the most complete (and usually 
most expensive) mathematical models. However, 
the more complicated the models, the more dif- 
ficult it is to estimate variations in the response 
due to variations in model parameters. 

c. Virtually all analytical methods have been de- 
veloped to solve differential equations or to repre- 
sent solutions of those equations. For most real 
systems, the solution of differential equations is 
accomplished via numerical techniques, using dif- 
ferencing techniques. Alternatively, the computer is 
used to numerically evaluate integral solutions or 
closed-form solutions. Volumes TM 5-858-3 and -4, 
GE-TEMPO (1974b), and AE WES (1972) reference 
computer programs pertaining to various weapon 
effects as well as descriptions of the phenomenol- 
ogy computer codes that calculate the free-field 
loadings, which the verification analyst may find 
useful. 

d. It is important that the verification analyst 
recognize that using more sophisticated computa- 
tional tools is equivalent to the reduction of (or the 
attempt to reduce) systematic uncertainties. In 
other words, these techniques are designed to mini- 
mize bias in the results. This reduces the ignorance 
factor and the uncertainties. An obvious conse- 
quence is the saving of dollars spent avoiding fa- 
cility overdesign that otherwise would occur due to 
use of poor computational tools. The analyst should 
also recognize that bias and random uncertainty 
can be reduced by using better data, i.e., more 
accurate loads and specifications of material prop- 
erties. 

e. Another resource available to the analyst are 
functional relationships derived from regression 
equations applied to experimental or computed 
data in which the response of systems is defined in 
terms of the important parameters of the system. 
These solution techniques often are empirical or 
semi-empirical, and are found in various sources, 
such as Crawford et al. (1974) and TM 5-585-3. The 
verification analyst will find most of this informa- 
tion unsuitable for performing deterministic ver- 
ification analyses; i.e., for calculating the mean 
response. Nevertheless for the probabilistic analy- 
ses discussed subsequently, these techniques may 
be quite useful. 



/. The verification analyst may avail himself of 
the following techniques for determining the mean 
response of systems that are subject to nuclear 
attack 

—Explicit solutions of differential equations 
—Numerical solutions of differential or 

integral equations 
—Semi-empirical relationships developed from 
regression analyses applied to 
experimental and calculated data 

The last technique is most often applied to the 
design of systems and the first method is usually 
too idealized to be much pratical usefulness; the 
second approach offers the greatest potential to 
generally address the verification problem. 

g. Compared to those methods available for cal- 
culating the mean response of a system, estimating 
the expected randomness due to natural variations 
in material properties and in geometries is grossly 
more difficult, a task for which little experience 
exists and few tools have been developed. Although 
there is awakening of interest in the probabilistic 
aspects of system response, most of the current 
statistical work is too fundamental to be of much 
direct application to hardened-facilities 
verifications. Therefore, what is lacking in formal 
computational procedures must be compensated for 
with a resourceful application of those analytical 
tools that are available. 

h. The deterministic methods previously outlined 
genrally involve the solution of differential equa- 
tions with variable coefficients, whereas deriving 
the probabilistic response of systems involves the 
solution of those same equations with both variable 
and stochastic coefficients. The latter property 
sharply limits the methods for calculating the re- 
sponse as a random variable. In verification analy- 
sis there are, for all practical purposes, only three 
fundamental techniques for calculating the random 
response of systems: 

—Monte Carlo methods 

—Engineering judgment 

—Use of simplified uncertainty analyses 

generally involving sensitivity (changes 
in response due to changes in system 
parameters such as stiffness, strength, 
etc.) 

The most versatile and reliable of these techniques 
is the Monte Carlo method, which can be applied to 
almost any problem that can be solved determinis- 
tically. However, there are certain negative aspects 
to this technique that should be considered (para 
5-4). Due to cost and time constraints or simply 
because no better alternative seems to exist, the 
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probabilities of the system response may be simply 
estimated by applying the judgment of acknowl- 
edged experts (para 5-5). On the other hand, if it is 
possible to explicitly and functionally express the 
relationship between the input and the response of 
a system, then direct differentiation of the function 
may provide a direct means for calculating uncer- 
tainties, and subsequently the survival probabilities 
of the systems (para 5-6). 

5-4. The Monte Carlo method. 

a. The correspondence between the experimental 
techniques described elsewhere in this volume and 
the mathematica approaches that seek to simulate 
the experimental method is achieved by a proce- 
dure called "mathematical experimentation." The 
form most widely used is called "Monte Carlo" or 
the Monte Carlo simulation. 

b. A simple application of the Monte Carlo pro- 
cess consists of the selection of random variable 
from an appropriate probability-density distribu- 
tion that describes the loading function and the 
geometrical and constitutive parameters of a sys- 
tem, and subsequently calculating the deterministic 
response for each of the random selections. 

c. More complicated applications of the Monte 
Carlo technique involve multiple systems interre- 
lated to each other. By linking the results of var- 
ious Monte Carlo solutions together, verification 
analyses of large complex systems can be 
performed. 

<L Although the Monte Carlo method can be used 
in many applications, there are disadvantages that 
must be considered: Cost, loss of visibility, and 
requirement for statistical input. In order to con- 
trol costs, the verification analyst may select a 
sophisticated analytical procedure for computing 
the mean response of the system (to reduce bias 
errors), and select an "equivalent" but simpler 
mathematical model for producing Monte Carlo t 
solutions (to provide information on the random 
nature of the response). 

e. In linked analyses, in which the results of the 
Monte Carlo problem are used as input to another 
Monte Carlo problem, the uncertainties in the final 
calculation may be governed by errors in the as- 
sumed nature of the randomness of the system 
parameters. To minimize costs it is desirable to 
reduce the number of Monte Carlo solutions as 
much as possible, but the use of a crude Monte 
Carlo mesh introduces further uncertainties. These 
can be mitigated by increasing the number of solu- 



tions, but the improvement in accuracy of the ran- 
dom response is relatively insensitive to this tech- 
nique. 

/. When Monte Carlo calculations are performed, 
the phenomenological relationship between the ran- 
dom variations in the system parameters and the 
random nature of the response may be obscured. 
Auxiliary sensitivity studies using simpler models 
can help in understanding the parent Monte Carlo 
calculations. 

g. The Monte Carlo method requires that the 
parameters of the system be defined statistically. 
These statistical data are obtained experimentally 
or are estimated. Errors can be reduced by con- 
ducting more extensive tests to determine the true 
nature of the system variables. 

h. The verification analyst will want to study the 
specialized techniques that have been developed for 
the practical implementation of the Monte Carlo 
method. One such well-known technique is Rowan 
(1974). 

t. A modified Monte Carlo approach can be 
adopted in which the partial derivatives in equa- 
tion 5-4 can be numerically obtained by randomly 
varying the independent variables Y n in the vicinity 
of the mean value to obtain the corresponding vari- 
ation in the function f. The rpocedure is a for- 
malization of the method discussed in paragraphs 

4- 2 e and / of TM 5-858-2, and noted in paragraph 

5- 6 c below. 

5-5. Engineering judgment. 

a. In many cases, especially during the early 
phases of hardness verification, it will be judged 
not feasible to embark on a Monte Carlo approach 
or to use the approach discussed in paragraph 5-6. 
It may be necessary to estimate the probabilistic 
nature of the system response whose mean value 
was obtained by deterministic computations. Be- 
cause of his or her familiarity with the solution 
technique and a grasp of the range of system pa- 
rameter variations, the analyst will often estimate 
the distribution of responses around the mean 
value. This estimate will draw almost exclusively 
on experience, coupled perhaps with auxiliary cal- 
culations that provided insight into the sensitivity 
of the parametrs in the analysis. When using this 
method, the analyst should avoid selecting the 
most adverse condition of each parameter to assess 
the severest system response. The probability of 
uniformly encountering the most detrimental ex- 
tremes of each parameter is very remote and cor- 
responds to a very small probability of 
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occurrence— much less than the probabilities of 
each independent parameter selection. 

5-6. Direct differentiation of functions. 

a. When the solutions to differential equations 
can be expressed in a functional form, then the 
uncertainties can be calculated by evaluating the 
partial derivatives in equation 5-4. For example, 
the atmospheric transmission of thermal energy is 
often expressed in the simple equation 

T=e- <7R s 

where is a coefficient and R, is the slant range. 
The total uncertainty will then be 



where fi f is the total uncertainty of the functional 
form of equation 5-5, and fl a and fl r, are the total 
uncertainties of o* and R, respectively. 



6. As a pratical matter, the simplified relation- 
ships in the response equations may not yield the 
best prediction of the mean response, since other 
important variables may be overlooked or not com- 
pletely represented. Most such simple relationships 
are probably more useful for estimating 
uncertainty than for calculating the mean response. 
It would be more appropriate to calculate the mean 
response by using one of the more complete deter- 
ministic procedures referenced in paragraph 5-2. 



c. For systems that are not grossly nonlinear, 
equation 5-4 may be solved numberically by vary- 
ing the parameters of the system in small incre- 
ments to evaluate the various partial derivatives, 
df/dXj. Usually, each independent parameter 
would be varied at least once above and once below 
its mean value so that the derivatives could be 
computed from polynomials fitted to the response 
calculations. 
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CHAPTER 6 
VERIFICATION METHODOLOGY 



6-1. Introduction. 

a. In order to perform a verification analysis, the 
analyst must define the extent of the system to 
which the verification procedures will be applied. 
The "system" to be verified may be a complete 
facility or a functional unit that forms part of the 
facility, or a submit, subsystem, or element. 

6. Most facilities are so complex that it is neither 
practical nor desirable to consider independently 
all elements or even all subsystems or systems in a 
verification analysis. Usually, many parts of a fa- 
cility will not be susceptible to every nuclear 
weapon effect or, occasionally, some parts will not 
be susceptible to any nuclear weapon effect. Thus, 
it is necessary to select (based on estimate of de- 
sign loads and resistances) which parts of the com- 
plete facility are susceptible and, of these, which 
are most prone to failure under a specified attack. 
Those parts that are neither susceptible nor threat- 
ened will not need further analysis. The remaining 
parts will then be subjected to a verification analy- 
sis as subsequently described. 

6-2. Diagrams, network logic, fault 
trees, and Boolean algebra. 

cu Procedure. Applying the various techniques 
described in previous sections of the manual is 
facilitated by following this procedure: 
—Construct functional block diagrams 
—Construct network logic or fault-tree diagrams 
—Write and simplify Boolean equations 
—Determine resistances and compute element 

uncertainties 
—Define loads and uncertainties 
—Compute subsystem and system uncertainties 

and compute survival probabilities where 

required. 

6. Functional block diagram. The functional block 
diagram is prepared to orient the analyst to the 
interrelationships between the physical or 
functional parts of a system. It is most commonly 
exemplified by schematic diagrams, which are then 
used to construct logic diagrams. 

c. Network logic and fault-tree diagrams. Net- 
work and tree diagrams, both of which can be 
useful in hardness-verification programs, are de- 
picted in figure 6-1. However, for complete facili- 
ties, systems, subsystems, and even elements, the 
tree diagram is most often used, since there is 
usually no complicated logic that would require 
network diagrams. Only tree analyses are discussed 
in this volume. 



d. Fault tree. An example of a generalized fault 
tree is shown in figure 6-2. It is diagrammatic 
representation of the interrelationships between 
the various failure modes of a system and the 
environments that potentially could induce failure. 
Generally, it can be said that a fault tree for a 
complete facility consists of N systems, each of 
which has M subsystems, each of which is in turn 
composed of Q elements. An element will have z 
functional failure modes, each of which is excited 
by S loads. In reality, each functional failure mode 
may have contributory components either in par- 
allel to each other or in series. Also, the fault tree 
is often expanded to show the resistances that 
react against the loads. The analyst may construct 
a tree for any point in the system, ranging from 
the entire facility to any element within a sub- 
system. Figure 6-2 is a fault tree for an entire 
facility. 

(1) When all of the elements that are to be 
included in the fault tree have been identified and 
diagrammed to show their relationship to each 
other (this information is obtained from the func- 
tional block diagram discussed in b above), then 
the failure modes for each element should be iden- 
tified if they bear some positive relationship to the 
mechanical, radiative, thermal, and other thermo- 
nuclear elements. The local environment, i.e., the 
loads exciting each response leading to a failure 
mode, should also be identified. The failure modes 
and loads are then added to the fault tree as shown 
in figure 6-2. More than one element, subsystem, or 
system can be considered in the analysis. The lar- 
ger the physical system involved, the more com- 
plete will be the fault tree. However, the tree will 
also become more complex. Initially, the fault trees 
are constructed without regard to simplification 
because their primary purpose is to provide a com- 
plete (and perhaps redundant) description of the 
system, its failure modes, the loads acting on it, 
and the resistances. 

(2) When constructing the fault trees, the over- 
all complexity may be minimized by removing 
superhard or supersoft systems, subsystems, or ele- 
ments from consideration; and by beginning at the 
lowest level of the tree that is practicable, i.e., at 
the element level. Nevertheless, for a system, sub- 
system, or element of any real complexity, it is 
inevitable that the first construction will contain 
redundancies and interdependencies that would not 
be immediately obvious. A simplifying procedure is 
required. 
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Figure 6-1. Types qf logic diagrams. 
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Figure 6-1. Typical generalized fault tree. 
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e. Boolean equations. Fortunately, the fault tree 
can be used to develop Boolean equations that will 
reveal interdependencies and eliminate redundan- 
cies. This is accomplished by writing the Boolean 
equations for the original fault tree and then re- 
ducing the equations to their simplest form. A new 
fault tree can then be constructed from the new 
Boolean equations, and this revision will be void of 
redundant information. 

(1) In order to implement Boolean algebra, the 
fault trees are usually presented in the particular 
form typified in figure 6-3. The symbols are called 
"gates;" the Boolean logic depends on the type of 
gate involved. Hammer (1972) describes the special- 
ized fault tree and Boolean algebra. An application 
to hardened systems is presented by Collins (1975). 

(2) An extremely important feature of Boolean 
algebra is that the equations can be used to des- 
ignate the survival or failure probability of the 
system under study. (See Collins, 1975.) 

6-3. Elements of verification analysis. 

a. Resistance. One of the principal purposes of 
this volume has been to establish that resistance is 
a random variable. The character of the load and 
the resistance are expressed as uncertainty factors 
that are defined as a measure of scatter (variance) 
normalized to the mean-squared response. More- 
over, uncertainty is generally composed of a ran- 
dom component (a reflection of natural variation) 
and a systematic component (an expression of ig- 
norance). How these uncertainties can be deter- 
mined and how survival probabilities can be com- 
puted were discussed in paragraph 5-2. 

6. Local load environment. The nuclear weapon 
effects exciting a system usually can be considered 
random variables, since free-field weapon effect 
loads are transmitted through the parts of a fa- 
cility whose transmission characteristics, 
determined by its geometrical and material prop- 
erties, are at best only statistically known. Hence, 
the local environments are probabilistic, and are 
modifications (amplification or attenuation) of one 
or more of the free-field, nuclear weapon-effect 
loads. If only the primary environment is known, 
than either the local environment must be arbi- 
trarily specified or it must be determined. If it is 
to be determined, the values must be derived by 
the use of transfer functions, which relate free- 
field environments to local-load environments. Be- 
cause of the nature of the real systems, the trans- 
fer functions are also random variables. The rela- 
tionship among these variables is shown in figure 
6-4. 



c. Transfer functions. Transfer functions link the 
primary nuclear weapon-effect loads to the local 
loads acting on the system, subsystem, or element 
under study. When performing an actual verifica- 
tion analysis of a system, the general procedure is 
first to complete verification of critical elements, 
then progress to subsystems, systems, and finally 
the facility, in building-block fashion. In so doing, 
however, defining the loads at the element level 
becomes more difficult, since the analyst will not 
yet have determined how the primary nuclear 
weapon-effect loads have been modified by the ac- 
tual intervening systems and subsystems. Thus, the 
verification analyst either must utilize the existing 
design loads or must develop transfer functions. 

(1) Whichever method is used depends on the 
particulars of the problem under investigation. For 
example, if the design-load specification represents 
the most advanced (best) estimate of the load, the 
analyst may continue to use this as the load speci- 
fication. However, if the outgoing verification anal- 
yses indicate that the original estimate was in 
error, a better estimate of the local loads would be 
required. Later, as progress is made by solving the 
fault tree equations, the transfer functions can be 
corroborated by analyses; also new verification at 
less comprehensive levels can be conducted when- 
ever the ongoing analyses indicate that the trans- 
fer functions were erroneously or incompletely 
specified. 

(2) It is again emphasized that transfer func- 
tions have their own uncertainties, consisting of 
random and systematic components. These uncer- 
tainties add another variation to the local load 
environment, beyond the variations (if any are as- 
sumed) existing in the free-field nuclear-weapon 
effect loads themselves. 

6-4. Correlations. 

a. An additional factor to be considered is that 
local loads acting on a system, subsystem, or ele- 
ment may be transmitted via one or more transfer 
functions, in series or in parallel, that may be 
correlated with each other. The simplest configura- 
tion comes from Boolean equations derived from 
the fault tree. Note that the Boolean algebra auto- 
matically accounts for serial and parallel connec- 
tions. However, it will not account for the degree 
of correlation between transfer functions; these 
must be calculated separately. 

b. A given system may be suceptible to a number 
of nuclear-weapon-effect loads. At the particular 
point where verification is being performed, there 
may be some degree of correlation between local 
environments, transfer functions, or resistances. 
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Figure 6-S. Typical Fault Tree with Network Logic Symbols 
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Figure 6-4. Loads, Transfer Functions, and Resistance. 
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The degree of correlation is specified by the cor- 
relation coefficient p u . Correlation coefficients vary 
in the range of — l<p<l. If strong positive cor- 
relation occurs, the system is more likely to survive 
than if the variables were not correlated at all. On 
the other hand, if strong negative correlation oc- 
curs, the system is less likely to survive than if the 
variables were not correlated at all. Therefore, in 
evaluating equation 5-9 the correlation coefficients 
between various elements of the local loads and the 
resistance must be determined. 

c. The covariance (which is a more general quan- 
tity to be used in equation 5-4) is defined as 



where Y n and Y m are any two of N variables cor- 
responding to the resistance function or the load 
function. The notation y indicates selected values of 
the model representing the random variable Y, and 
y is the mean of those selected values. The quan- 
tity Oy is the standard deviation of the variable Y. 
Equation 5-7 is evaluated by randomly selecting 
values for the variables y n and y m . For example, y n 
may represent random values for the stiffness of 
concrete, and y m may represen random values for 
the damping properties of the same concrete. If the 
stiffness and damping are totally uncorrected (i.e., 
totally independent of each other), the quantity 
within the summation of equation 5-7 will tend to 
zero as N becomes large and Pnm ~0 ; on the other 
hand, if the stiffness and damping are negatively 
correlated (i.e., if selections of y„>y m are more 
often than not produce the condition that y m <y m , 
then the quantity within the summation will tend 
toward a negative number and P y will be a nega- 
tive number for large values of N. Coversely, posi- 
tive correlation (p nm >0) would indicate that as the 
concrete stiffness increases, the damping also in- 
creases. The limiting conditions for the correlation 
coefficient are -1 <p nin <1. 
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APPENDIX A 
BINOMIAL DISTRIBUTION 



This appendix contains two tables in support of 
statistical analysis, Chapter 3. Use Table A-l in 
determining the probability of s=r successes in n 
trials, table A-2 for s>r successes in n trials. 
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TABLE A-l. INDIVIDUAL TERMS 



PROBABILITY OF REALIZING s * r SUCCESSES IN n TRIALS 
P(s = r/n ' P> * r!(n - r)i P q 



p 


n = 2 


n = 2 


n = 3 


n = 3 


n = 3 


n = 4 






r = l 


r - u 


r = z 


r = l 


r = u 


r = 3 


P 


• 01 


• 0196 


• 9601 


• 0003 


• 0294 


• 9703 


• 0000 


• 01 


• 02 


• 0392 


• 9604 


• 0012 


• 0576 


• 9412 


• 0000 


• 02 


• 03 


• 0562 


.9409 


• 0026 


• 0647 


• 9127 


• 0001 


• 03 


• 04 


• 0766 


• 9216 


• 0046 


• 1106 


• 6647 


• 0002 


• 04 


• 05 


• 0950 


• 9025 


• 0071 


• 1354 


• 6574 


• 0005 


• 05 


• 06 


• 1126 


• 6636 


• 0102 


• 1590 


• 6306 


.0006 


• 06 


• 07 


• 1302 


• 6649 


• 0137 


• 1616 


• 6044 


• 0013 


.07 


• 06 


• 1472 


• 6464 


• 0177 


• 2031 


• 7767 


• 0019 


• 06 


• 09 


• 1636 


• 6261 


• 0221 


• 2236 


• 7536 


• 0027 


• 09 


• to 


• 1600 


• 6100 


• 0270 


• 2430 


• 7290 


• 0036 


• 10 


• 1 1 


• 1956 


• 7921 


• 0323 


• 2614 


• 7050 


.0047 


• 1 1 


• 12 


• 2112 


• 7744 


• 0360 


• 2766 


• 6615 


• 0061 


• 12 


• 13 


• 2262 


.7569 


• 0441 


• 2952 


• 6565 


• 0076 


• 13 


• 14 


•2406 


.7396 


• 0506 


• 3106 


• 6361 


• 0094 


• 14 


• 15 


• 2550 


• 7225 


.0574 


• 3251 


• 6141 


• 0115 


• 15 


• 16 


• 2666 


• 7056 


• 0645 


• 3367 


• 5927 


• 0136 


• 16 


• 17 


• 2622 


• 6669 


• 0720 


• 3513 


• 5716 


• 0163 


• 17 


• 16 


• 2952 


• 6724 


• 0797 


• 3631 


• 5514 


• 0191 


• 16 


• 19 


• 3076 


• 6561 


• 0677 


• 3740 


• 5314 


• 0222 


• 19 


• 20 


• 3200 


• 6400 


• 0960 


• 3640 


• 5120 


• 0256 


• 20 


• 21 


• 3316 


• 6241 


• 1045 


• 3932 


• 4930 


• 0293 


• 21 


• 22 


• 3432 


• 6064 


• 1133 


• 4015 


• 4746 


• 0332 


• 22 


• 23 


• 3542 


.5929 


• 1222 


.4091 


• 4565 


.0375 


.23 


• 24 


• 3646 


.5776 


• 1313 


.4159 


• 4390 


• 0420 


• 24 


• 25 


• 3750 


• 5625 


• 1406 


.4219 


• 4219 


.0469 


• 25 


• 26 


• 3646 


.5476 


• 1501 


• 4271 


• 4052 


• 0520 


• 26 


• 27 


• 3942 


• 5329 


.1597 


• 4316 


• 3690 


.0575 


• 27 


• 26 


• 4032 


• 5164 


• 1693 


• 4355 


.3732 


• 0632 


• 26 


• 29 


• 4116 


• 5041 


.1791 


• 4366 


• 3579 


.0693 


• 29 


• 30 


• 4200 


.4900 


• 1690 


• 4410 


• 3430 


.0756 


• 30 




• 4276 


• 4761 


.1969 


• 4426 


• 3265 


.0622 






• 4352 


.4624 


• 2069 


• 4439 


• 3144 


• 0691 


'.11 




• 4422 


.4469 


.2169 


• 4444 


• 3006 


• 0963 






• 4466 


• 4356 


• 2269 


• 4443 


• 2675 


• 1036 






• 4550 


.4225 


• 2369 


• 4436 


• 2746 


• 1115 






• 4606 


• 4096 


• 2466 


• 4424 


• 2621 


• 1194 






• 4662 


.3969 


.2567 


• 4406 


• 2500 


• 1276 






• 4712 


• 3644 


• 2666 


• 4362 


• 2363 


• 1361 






• 4756 


• 3721 


• 2763 


• 4354 


• 2270 


• 1447 






• 4600 


.3600 


• 2660 


• 4320 


• 2160 


• 1536 






• 4636 


• 3461 


.2975 


• 4262 


• 2054 


.1627 






• 4672 


• 3364 


• 3069 


• 4239 


.1951 


.1719 






• 4902 


• 3249 


.3162 


• 4191 


• 1652 


• 1613 






• 4926 


• 3136 


• 3252 


• 4140 


.1756 


• 1906 






• 4950 


• 3025 


• 3341 


• 4064 


• 1664 


• 2005 






• 4966 


• 2916 


• 3426 


• 4024 


.1575 


.2102 






• 4962 


• 2609 


• 3512 


• 3961 


• 1469 


• 2201 






• 4992 


• 2704 


• 3594 


• 3694 


• 1406 


• 2300 






• 4996 


• 2601 


• 3674 


• 3623 


• 1327 


• 2400 






• 5000 


• 2500 


.3750 


• 3750 


• 1250 


• 2500 


1 50 
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n = 4 


n = 4 


n = 4 


n = 5 


n = 5 


n = 5 




P 


r = 2 


r = 1 


r = 0 


r = 4 


r = 3 


r = 2 


P 


• 01 


• 0006 


• 0366 


• 9606 


• 0000 


• 0000 


• 0010 


• 01 


• 02 


• 0023 


• 0753 


• 9224 


• 0000 


• 0001 


• 0036 


• 02 


• 03 


• 0051 


• 1095 




• 0000 


• 0003 


• 0062 


• 03 


• 04 


• 0066 


• 1416 


*6493 


• 0000 


• 0006 


• 0142 


• 04 


• 05 


• 0135 


• 1715 


*6145 


• 0000 


• 0011 


• 0214 


• 05 


• 06 


• 0191 


• 1993 


*7607 


• 0001 


• 0019 


• 0299 


• 06 


• 07 


• 0254 


• 2252 


*7461 


• 0001 


• 0030 


• 0394 


• 07 


• 06 


• 0325 


• 2492 


*7 164 


• 0002 


• 0043 


• 0496 


• 06 


• 09 


• 0402 


• 2713 


*6657 


• 0003 


• 0060 


• 0610 • 


• 09 


• 10 


• 0466 


• 2916 


*6561 


• 0004 


• 0061 


• 0729 


• 10 


• 1 1 


• 0575 


• 3102 


*6274 


• 0007 


• 0105 


• 0653 


• 11 


• 12 


• 0669 


• 3271 


*5997 


• 0009 


• 0134 


• 0961 


• 12 


• 13 


• 0767 


• 3424 


*5729 


• 0012 


• 0166 


• 1113 


• 13 


.14 


• 0670 


• 3562 


*5470 


• 0017 


• 0203 


• 1247 


• 14 


• 15 


• 0975 


• 3665 


°5220 


• 0022 


• 0244 


• 1362 


• 15 


• 16 


• 1064 


• 3793 


*4979 


• 0026 


• 0269 


• 1517 


• 16 


• 17 


• 1195 


• 3666 


# 4746 


• 0035 


• 0336 


• 1652 


• 17 


• 16 


• 1307 


• 3970 


*4521 


• 0043 


• 0392 


• 1766 


• 16 


• 19 


• 1421 


• 4039 


*4305 


• 0053 


• 0450 


• 1919 


• 19 


• 20 


• 1536 


• 4096 


4096 


• 0064 


• 0512 


• 2046 


• 20 


• 21 


• 1651 


• 4142 


*3695 


• 0077 


• 0576 


• 2174 


• 21 


• 22 


• 1767 


• 4176 


3702 


• 0091 


• 0646 


• 2297 


• 22 


• 23 


• 1662 


• 4200 


3515 


• 0106 


• 0721 


• 2415 


• 23 


• 24 


• 1996 


• 4214 


*3336 


• 0126 


• 0796 


• 2529 


• 24 


• 25 


• 2109 


• 4219 


*3164 


• 0146 


• 0679 


• 2637 


• 25 


• 26 


• 2221 


• 4214 


*2999 


• 0169 


• 0962 


• 2739 


• 26 


• 27 


• 2331 


• 4201 


*2640 


• 0194 


• 1049 


• 2636 


• 27 


• 26 


• 2439 


• 4160 


# 2667 


• 0221 


• 1136 


• 2926 


• 26 


• 29 


• 2544 


• 4152 


2541 


• 0251 


• 1229 


• 3010 


• 29 


• 30 


• 2646 


• 4116 


*2401 


• 0263 


• 1323 


• 3067 


.30 


• 31 


• 2745 


• 4074 


*2267 


• 0319 


• 1416 


• 3157 


• 31 


• 32 


• 2641 


.4025 




• 0357 


• 1515 


• 3220 


• 32 


• 33 


• 2933 


• 3970 


•|jj5 


• 0397 


• 1613 


• 3275 


• 33 


• 34 


• 3021 


.3910 




• 0441 


• 1712 


• 3323 


• 34 


.35 


• 3105 


• 3645 




• 0466 


• 1611 


• 3364 


• 35 


• 36 


• 3165 


• 3775 




.0537 


• 1911 


• 3397 


• 36 


• 37 


• 3260 


• 3701 




.0590 


• 2010 


• 3423 


• 37 


• 36 


• 3330 


• 3623 




.0646 


.2109 


• 3441 


• 36 


.39 


• 3396 


• 3541 




• 0706 


• 2207 


• 3452 


.39 


• 40 


• 3456 


• 3456 




• 0766 


• 2304 


• 3456 


• 40 


• 41 


• 3511 


• 3366 




• 0634 


.2399 


.3452 


• 41 


• 42 


• 3560 


.3276 




.0902 


• 2492 


• 3442 


• 42 


• 43 


• 3604 


.3165 




• 0974 


.2563 


• 3424 


• 43 


.44 


• 3643 


• 3091 




• 1049 


• 2671 


• 3400 


• 44 


.45 


• 3675 


• 2995 




• 1126 


.2757 


• 3369 


• 45 


.46 


• 3702 


• 2697 




• 1209 


• 2636 


• 3332 


• 46 


.47 


.3723 


.2799 




• 1293 


• 2916 


• 3269 


• 47 


.46 


• 3736 


.2700 




• 1360 


.2990 


• 3240 


• 46 


.49 


• 3747 


.2600 




• 1470 


.3060 


• 3165 


• 49 


• 50 


• 3750 


• 2500 


I0625 


• 1562 


.3125 


• 3125 


.50 
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• vaav 


1 31 A 
• 1 ? 1 o 


?96Jt 


ac 
• c? 




1A9A 


PPI 9 


0053 
• V V J J 


0379 


1 fldl 
• 1 aca 


1 A A 1 
• ?W A 1 


3^ 
• CO 


?7 


3A3A 


P073 
• C V f 3 


0063 
• WO J 


0A?5 
• V*C J 


1531 
• 1331 


1 1 A5 
.31 V3 


• Cf 


PA 


376? 


1 935 


A07 A 
• V V » H 


047 A 
• V*t * o 


• loa" 


1 1 &A 
. 3 IOW 


3 A 
• CO 


*P9 


36A5 


1 AAA 
• lOVA 


AAA7 
• VVOr 


• Vaaa 


i 7 
• 1 1 AO 


13 A4 

. 3c0o 


30 

• c? 




360? 


1 6A 1 
• I OO I 


A t A3 
• VI Vc 


ACQ! 
• W373 


• 1 03C 


.3CA1 


1A 
• 30 




• 33 1 3 


1 56A 


. V 1 I ▼ 


A66A 
• VOO V 


1 OA7 

• 1 "ar 


13A7 


• 31 




34P1 


t A5A 
• i *3* 


A 1 37 
• VI Jf 


A737 
• V f C f 


• CVOl 


13 A A 
. ?COA 


ia 
• ac 




33?5 


i ten 
. & asu 




A700 


9 1 ka 
• cloc 


1303 
. 3C?C 


ti 

• 33 




- 3PP6 


1 P5P 


0 1 AO 
. V A O V 


AA7 3 
• VO f J 


aa^n 
• ccov 


139A 
. 3C^V 


1A 
• 3A 


1 35 


• 3124 


• 1160 


.0205 


• 0951 


.2355 


.3260 


• 35 




• 30cO 


1 A 7 A 

• 1 V r A 


. VC JC 


t a la 
• IV ac 


9 It AX. 

.c44o 


• 3261 


• 36 




• 2914 


• 0992 


• 0262 


• 1116 


• 2533 


• 3235 


.37 




• 2606 


.0916 


• 0295 


• 1202 


• 2616 


• 3201 


• 36 




• 2700 


• 0645 


• 0330 


• 1291 


• 2693 


• 3159 


.39 




• 2592 


• 0776 


• 0369 


• 1362 


• 2765 


• 3110 


• 40 




• 2464 


• 0715 


• 0410 


• 1475 


• 2631 


• 3055 


• 41 




• 2376 


• 0656 


• 0455 


.1570 


• 2691 


• 299A 


• 42 




• 2270 


• 0602 


.0503 


• 1666 


.2945 


• 2926 


• 43 




• 2164 


• 0551 


.0554 


• 1763 


.2992 


• 2656 


.44 




.2059 


• 0503 


• 0609 


• 1661 


.3032 


• 2780 


.45 




• 1956 


• 0459 


.0667 


• 1956 


.3065 


.2699 


• 46 




• 1054 


• 0416 


.0729 


.2056 


.3091 


• 2615 


.47 




.1755 


• 0360 


.0795 


• 2153 


• 3110 


• 2527 


.46 




.1657 


.0345 


• 0664 


• 2249 


• 3121 


.2436 


.49 




• 1563 


.0313 


.0937 


• 2344 


• 3125 


• 234A 


.50 
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n — A 
n — O 


n - A 
11 — O 


n - 7 


n = 7 
n — / 


n - 7 


n - 7 


P 


p 


r = 1 


r = 0 


r = 6 


r = 5 


r = 4 


r = 3 




• 01 


• 0571 


• 9415 


• 0000 


• 0000 


• 0000 


• 0000 


.01 


• 02 


• 1085 


• 6656 


• 0000 


• 0000 


• 0000 


.0003 


• 02 


• 03 


• 1546 


• 6330 


• 0000 


• 0000 


• 0000 


• 0006 


• 03 


• 0* 


• 1957 


• 7626 


• 0000 


• 0000 


• 0001 


.0019 


.04 


• 05 


• 2321 


.7351 


• 0000 


• 0000 


• 0002 


• 0036 


.05 


• 06 


• 2642 


• 6699 


• 0000 


• 0000 


• 0004 


.0059 


• 06 


• 07 


• 2922 


.6470 


• 0000 


• 0000 


• 0007 


• 0090 


.07 


• 08 


• 3164 


• 6064 


• 0000 


• 0001 


• 0011 


• 0126 


• 06 


.09 


• 3370 


• 5679 


• 0000 


• 0001 


• 0017 


.0175 


• 09 


• 10 


• 3543 


• 5314 


• 0000 


• 0002 


• 0026 


• 0230 


• 10 


• 1 1 


• 3665 


• 4970 


• 0000 


• 0003 


• 0036 


• 0292 


• 11 


• 12 


• 3600 


• 4644 


• 0000 


• 0004 


• 0049 


• 0363 


• 12 


• 13 


• 3666 


.4336 


• 0000 


• 0006 


• 0066 


• 0441 


• 13 


• 14 


• 3952 


• 4046 


• 0000 


• 0006 


• 0066 


• 0525 


• 14 


• 15 


• 3993 


.3771 


• 0001 


• 0012 


• 0109 


• 0617 


• IS 


• 16 


• 4015 


• 3513 


• 0001 


• 0016 


• 0136 


• 0714 


• 16 


• 17 


• 4016 


• 3269 


• 0001 


• 0021 


• 0167 


• 0616 


• 17 


• IS 


• 4004 


• 3040 


• 0002 


• 0027 


• 0203 


• 0923 


• 16 


• 1* 


• 3975 


• 2624 


• 0003 


• 0034 


• 0242 


• 1033 


• 19 


• 20 


• 3932 


• 2621 


• 0004 


• 0043 


• 0267 


• 1147 


• 20 


•21 


• 3677 


• 2431 


• 0005 


• 0054 


• 0336 


• 1263 


• 21 


• 22 


• 361 1 


• 2252 


• 0006 


• 0066 


• 0369 


• 1379 


• 22 


• 23 


• 3735 


• 2064 


• 0006 


• 0060 


• 0447 


• 1497 


• 23 


• 24 


• 3651 


• 1927 


• 0010 


• 0097 


• 0510 


• 1614 


• 24 


• 25 


• 3560 


• 1760 


• 0013 


• 0115 


• 0577 


• 1730 


• 25 


• 26 


• 3462 


• 1642 


• 0016 


• 0137 


• 0646 


• 1645 


• 26 


• 27 


• 3356 


• 1513 


• 0020 


• 0161 


• 0724 


• 1956 


• 27 


• 28 


• 3251 


• 1393 


• 0024 


.0187 


• 0603 


• 2065 


• 26 


• 29 


• 3139 


• 1261 


• 0030 


• 0217 


• 0666 


• 2169 


.29 


• 30 


• 3025 


• 1176 


• 0036 


• 0250 


.0972 


• 2269 


• 30 


• 31 


• 2909 


• 1079 


• 0043 


• 0266 


• 1062 


• 2363 


• 31 


• 32 


• 2792 


• 0969 


• 0051 


• 0326 


• 1154 


• 2452 


• 32 


• 33 


• 2673 


• 0905 


• 0061 


• 0369 


• 1246 


• 2535 


• 33 




• 2555 


• 0627 


• 0071 


• 0416 


• 1345 


• 2610 


• 34 


1 33 


• 2437 


• 0754 


.0064 


• 0466 


.1442 


.2679 


.35 




• 2319 


• 0667 


• 0096 


• 0520 


• 1541 


• 2740 


• 36 




• 2203 


.0625 


• 0113 


• 0576 


• 1640 


.2793 


.37 




• 2069 


.0566 


.0131 


• 0640 


.1739 


• 2636 


• 36 




• 1976 


.0515 


.0150 


• 0705 


.1636 


.2875 


.39 




• 1666 


• 0467 


.0172 


.0774 


.1935 


• 2903 


.40 




.1759 


• 0422 


.0196 


.0647 


• 2031 


• 2923 


.41 




• 1654 


.0361 


.0223 


.0923 


• 2125 


• 2934 


.42 




.1552 


• 0343 


• 0252 


.1003 


• 2216 


• 2937 


• 43 




• 1454 


• 0306 


• 0264 


.1066 


• 2304 


• 2932 


• 44 




• 1359 


.0277 


• 0320 


.1172 


.2366 


.2918 


.45 




• 1267 


• 0246 


• 0356 


• 1261 


• 2466 


.2897 


.46 




.1179 


• 0222 


• 0400 


• 1353 


.2543 


.2667 


.47 




• 1095 


• 0196 


• 0445 


.1447 


.2612 


630 


.48 




• 1014 


• 0176 


.0494 


• 1543 


• 2676 


.2766 


.49 




• 0936 


• 0156 


.0547 


• 1641 


• 2734 


.2734 


.50 
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• 0000 
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• 0000 
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• 0000 
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• 33tO 
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A A A « 
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AA 


AO 
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A A A A 

• 0000 


AAAA 

.0000 


AA A3 
.0002 


AO 
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773A 
• 3f CO 
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AAAA 
.0000 


AAAA 
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# 1 A 


• J * 
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• 302/ 
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A A A A 
.0000 


AAAA 
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• 3701 
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. COO f 
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.0000 
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A A AO 
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• 3TCO 


. 3 f f D 


AAAA 
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. VO 1 c 




• J J 


t 07A 


YOAC 
.3V0D 


7ATO 
. 3C f 7 


AAAA 
. 0000 


A A A3 

• 000c 


A A i O 
.001™ 


•}! 


•J' 


3AOT 

• 20t7 


• 3tO0 


73 A A 

• 3200 


AAAA 
. 0000 


A A A3 

• 000c 


A A 31 
• 0020 




"It 


• CCCO 


tott 
• 3 TOD 


3QC t 
.CVDI 


AAAA 
. 0000 


A A A 7 

• 0003 


A At C 

.003D 


•}: 




.23T1 


7 AO t 
.3071 


97 fl A 

• cf ic 


AAAA 
.0000 


A A AC 

• 000D 


A A AC 

.00CD 




•[J 


ia in 
• 3030 


DAOI 
• COTD 


AAAA 
. 0000 


A AAl 

• OOOO 


A ACA 

. 00DO 


" 1 A 




• coo j 


«3'D0 


33AA 

• ccoo 


AAAA 
.0001 


A A AO 

• 000t 


A A 7 A 
. 00 fC 
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AAAA 
.0001 
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A AOS 


DO 
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• 1 TcO 


AA A t 

.0001 


AAl 1 

.001 D 
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• 0113 


Dl 
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307C 


1A&A 
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• 1 1 Df 


A A A3 
• 0002 


A A f O 
«00lT 


n t 77 
• 0137 
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Pi 
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• 300 f 


.33D0 


1 AAC 
• 1 OOD 


AA A3 
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• V VCD 
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• 010D 




da 
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7377 

.3237 


« AAC 
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A A A 7 

• 0003 


A A7 A 
.0031 


A t Ol 
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DA 
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so 


• 3113 


li it 
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A 11C 

• 133D 
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• OOOC 


A A 7 A 

.0030 


A37 % 

• 0231 
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• 31 90 
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• IclD 
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AAA7 

.00*7 


A37 A 
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DA 




• 31 7 C 


•cooo 


f 1 AC 

• 1 1 0D 


A AAl 

• 0000 


A AC A 

. 00DO 


A 7 t 7 

• 0313 




DA 


• 3lOD 


37 1 1 
.2731 


t AAl 

• 1003 


AAAA 

• 0000 


A A 7 A 
.00 '0 


A 71 A 

• 0300 


DA 


3 o 


1 i AA 

. 3100 


31 A A 

.2000 


AO A A 

.0*910 


A A * A 
• 0010 


AAAA 

• OOOC 


AAA! 

• ocii 


DO 


«A 

• 


• 31 f t 


• 247 1 


A A 3 A 

• 0024 


A A f 3 

• 0012 


.0 1 00 


A A 1 T 

• 0CO7 


*a 




7 t oa 
• 31 DO 


3 7 A 3 

• 2302 


A7 AC 

. 07 CD 


AAll 

• 0013 


At 1 A 

.0110 


AC37 

• 0D27 


*7f 


73 
•™ 


• 3127 


3 3 t K 


.0072 


• 001t 


At tS 

.0139 


A OA * 

• ODTl 


*73 




t AAA 

• 30OO 


3 AO A 

• 20t0 


A1A1 

. oooo 


A A 3 7 

• 0023 


A t 1 3 

• OlOc 


aico 
• 0ODT 




7 


1AAA 

.3000 


1 047 
• 1tO# 


. 0DCO 


A A 3 A 

• 0020 


A 4 A A 

.0100 


AT 73 

• 0 f 32 


7A 


1 35 


• 2965 


• 1646 


• 0490 


.0033 


.0217 


.0606 


1 33 




• 2922 


• 1732 


• 0440 


• 0040 


.0250 


.0666 






• 2653 


• 1619 


• 0394 


• 0046 


• 0265 


.0971 






.2776 


• 1511 


.0352 


• 0057 


• 0324 


.1056 






• 2696 


.1407 


• 0314 


• 0067 


• 0367 


.1147 






• 2613 


• 1306 


• 0260 


.0079 


.0413 


.1239 






• 2524 


• 1211 


• 0249 


• 0092 


• 0463 


.1332 






• 2431 


• 1119 


• 0221 


• 0107 


.0517 


• 1426 






• 2336 


• 1032 


• 0195 


.0124 


• 0575 


• 1525 






• 2239 


• 0950 


• 0173 


• 0143 


• 0637 


• 1622 






•2140 


• 0672 


• 0152 


• 0164 


• 0703 


.1719 






• 2040 


.0796 


• 0134 


• 0166 


.0774 


• 1616 






• 1940 


.0729 


• 0117 


.0215 


.0646 


• 1912 






• 1640 


.0666 


• 0103 


.0244 


• 0926 


• 2006 






.1740 


• 0604 


• 0090 


.0277 


• 1006 


• 2096 






• 1641 


• 0547 


• 0076 


.0312 


.1094 


• 2167 
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p 


n = ft 
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11 — o 
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n - ft 
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p 


r = 4 


r = 3 


r = 2 


r = 1 


r « 0 


r = 8 


• 01 


• 0000 


• 0001 


• 0026 


• 0746 


• 9227 


• 0000 


• 01 


• 02 


• 0000 


• 0004 


• 0099 


• 1369 


• 6506 


• 0000 


• 02 


• 03 


• 0001 


• 0013 


• 0210 


• 1939 


• 7637 


• 0000 


• 03 


• 04 


• 0002 


• 0029 


• 0351 


• 2605 


• 7214 


• 0000 


• 04 


• 05 


• 0000 


• 0054 


• 0515 


• 2793 


• 6634 


• 0000 


• 05 


• 06 


• 0007 


• 0069 


• 0695 


• 31 13 


• 6096 


.0000 


• 06 


• 07 


• 0013 


• 0134 


• 0666 


• 3370 


• 5596 


• 0000 


• 07 


• 00 


• 0021 


• 0169 


• 1067 


• 3570 


• 5132 


• 0000 


• 06 


• 09 


• 0031 


• 0255 


• 1266 


• 3721 


• 4703 


.0000 


• 09 


• 10 


• 0046 


• 0331 


• 1466 


• 3626 


• 4305 


• 0000 


• 10 


• 11 


• 0064 


• 0416 


• 1664 


• 3692 


• 3937 


• 0000 


• 11 


• 12 


• 0067 


• 0511 


• 1672 


• 3923 


• 3596 


• 0000 


• 12 


• 13 


• 0115 


• 0613 


• 2052 


• 3923 


• 3262 


• 0000 


• 13 


• 10 


• 0147 


• 0723 


• 2220 


• 3697 


• 2992 


• 0000 


• 14 


• 15 


• 0165 


• 0639 


• 2376 


• 3647 


• 2725 


• 0000 


• 15 


• 16 


• 0220 


• 0959 


• 2516 


• 3777 


.2479 


• 0000 


• 16 


• 17 


• 0277 


• 1064 


• 2646 


• 3691 


• 2252 


• 0000 


• 17 


• 10 


• 0332 


• 1211 


• 2756 


• 3590 


• 2044 


• 0000 


• 16 


• 19 


.0393 


• 1339 


• 2655 


• 3477 


• 1653 


• 0000 


• 19 


• 20 


• 0459 


• 1466 


•2936 


• 3355 


• 1676 


• 0000 


• 20 


•21 


• 0530 


• 1596 


• 3002 


• 3226 


• 1517 


• 0000 


• 21 


•22 


• 0607 


• 1722 


• 3052 


• 3092 


• 1370 


• 0000 


• 22 


•23 


• 0669 


• 1644 


• 3067 


• 2953 


• 1236 


• 0001 


• 23 


•20 


.0775 


• 1963 


• 3106 


• 2612 


• 1113 


• 0001 


• 24 


•25 


• 0665 


•2076 


• 3115 


• 2670 


• 1001 


• 0001 


•25 


• 26 


.0959 


•2164 


• 3106 


• 2527 


• 0699 


• 0001 


• 26 


•27 


• 1056 


• 2265 


• 3069 


• 2366 


• 0606 


• 0002 


•27 


• 20 


• 1156 


• 2379 


• 3056 


• 2247 


• 0722 


• 0002 


• 26 


• 29 


• 1256 


• 2464 


• 3017 


• 2110 


• 0646 


• 0003 


• 29 


• 30 


• 1361 


• 2541 


• 2965 


• 1977 


• 0576 


• 0004 


• 30 


• 31 


• 1465 


• 2609 


• 2904 


• 1647 


• 0514 


• 0005 


• 31 


• 32 


• 1569 


• 2666 


• 2635 


• 1721 


• 0457 


• 0007 


• 32 


• 33 


• 1673 


• 2717 


• 2756 


• 1600 


• 0406 


• 0006 


• 33 




• 1775 


• 2756 


• 2675 




• 0360 


• 001 1 


• 34 


l55 


.1675 


• 2766 


• 2567 


Il373 


• 0319 


• 0013 


.35 




.1973 


• 2605 


• 2494 




• 0261 


• 0016 


• 36 




.2067 


• 2615 


• 2397 




• 0246 


• 0020 


.37 




.2157 


• 2615 


• 2297 




• 0216 


• 0024 


.36 




.2242 
.2322 


•2606 


• 2194 




• 0192 


• 0029 


.39 




• 2767 


• 2090 




• 0166 


• 0035 


• 40 




.2397 


.2759 


• 1965 




• 0147 


• 0042 


.41 




• 2465 


• 2723 


• 1660 




• 0126 


• 0051 


• 42 




.2526 


• 2679 


.1776 




• 0111 


• 0060 


• 43 




• 2560 


• 2627 


• 1672 




.0097 


.0071 


.44 




• 2627 


• 2566 


.1569 




• 0064 


.0063 


.45 




•2665 


• 2503 


.1469 




• 0072 


.0097 


.46 




• 2695 


.2431 


.1371 




• 0062 


.0114 


.47 




•2717 


• 2355 


• 1275 




• 0053 


.0132 


.46 




• 2730 


• 2273 


• 1163 




• 0046 


.0153 


.49 




• 2734 


• 2166 


• 1094 




.0039 


.0176 
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p 
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r — / 
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A A A A 


A A A A 
• 0000 


AAAA 
.0000 


AAAA 

.0000 


A A A 1 
.000 1 


A A 11 

• 0034 


A 1 
.01 


A S 


A A A A 
• OOOO 


A A A A 
• 0000 


AAAA 
.0000 


AAAA 
. 0000 


AAAt 
. 0000 


A 4 SC 

• 0 1 23 


A S 

.02 




A A A A 

• OOOO 


AAAA 

• 0000 


AAAA 

• 0000 


A A A • 

.0001 


AA t A 

• 001* 


A 3 1 9 

• 0262 


A 1 

.03 




A A A A 


AAAA 

• 0000 


AAAA 

.0000 


AAAt 

. 0003 


A A A 3 

• 0042 


A A 11 
. 0433 


A A 

.04 


* AC 


A AAA 

• OOOO 


AAAA 

• 0000 


AAAA 
.0000 


AAAt 

.0006 


A ATT 

• 0077 


A t SB 

• 002* 


AC 

.03 


A 


A A A A 


AAAA 

• 0000 


A A A 4 
.0001 


A A 1 3 

• 0012 


A « 3C 

• 0 123 


AAAA 

• 0040 


At 

. 00 


" A? 


AAA A 


AAAA 
• 0000 


A A A3 
• 00 02 


A A3 1 
• OOCl 


A 1 AJ» 
.01 OO 


1 A At 4 

. 1 00 1 


AT 
.Or 


AC 


A A A A 
• 0000 


AAAA 

• 0000 


A A A 1 
. 0003 


A A 1 A 
• 0034 


A St 1 

• 02O1 


4 SAC 
. 1203 


A A 
.00 


* AO 


A AAA 
• 0000 


AAAA 

• 0000 


A A AC 
. 0003 


A ACS 
« 003c 


A 1 A A 
• 0340 


4 CAT 
• 1 30 7 


AO 

. o* 


" 1 A 


AAA A 
• OOO O 


A AAI 
• 0001 


AAAA 
• O OOO 


AAT A 
• OO f 4 


AAA J. 

• 044© 


1 T3S 
• 1 7 22 


1 A 
.10 


# 11 

* J * 


A A A A 
• 0000 


A A A 1 
• 0001 


A A 1 1 
• 00 1 3 


A 1 A 1 
. 0 103 


ASIA 

• 033O 


4 OST 
. 1*27 


4 4 
.11 




A A A A 

• 0000 


A A A S 

.0002 


A A 4 A 

.001* 


A • 11 

.0130 


At T A 

• 0074 


S4 4 A 

.211* 


• s 

.12 


* if 


A AAA 
• 0000 


A A A V 
• 0003 


A A3T 
.002 r 


A 1 TO 
• 0 1 7* 


AAAA 
• OOOO 


ssoc 
.22*3 


« i 
.13 


• J J 


A A A A 
• 0000 


AAAA 

• 000* 


A A 17 
. 003 f 


A33A 

«0ccO 


• 0133 


S ACC 
.2433 


4 A 

.14 


• J J 


A A A A 
• 0000 


A A At 

• 0006 


A ACA 
.0030 


A3A1 
. 0cO3 


1 AAkO 

. 1 0O* 


SCOT 
.23* 7 


.13 


•f: 


AAA 4 
• 000 1 


AAAA 
• 0000 


. OOOO 


At AC 

. 03C3 


4 SAO 

. 120* 


ST S A 
.2 7 20 


. io 




A A A 1 

• 000 1 


AAI ) 

• 001c 


A A AC 
.0003 


A A 1 C 
• 0C13 


i no 
. 1 341 


SAS1 

.2023 


.17 


* 1 A 

• Jr 


A A A 4 
• 000 1 


• 0010 


A 4 A A 

.01 00 


A AAA 

• 0410 


• AAA 

. 140* 


SOA A 

.2*00 


• i 
. 10 




A A A3 
• 0002 


A A3 1 

• 0021 


A 1 11 
• 013C 


ACT1 

.03 7 3 


4 fckST 

. 1027 


SOT 1 

• 2* 7 3 


4 O 
.1* 


sa 


AA A V 

• 0003 


A A3! 
. 0020 


A« i( 

«01O3 


• oooi 


4 7fc>S 

. 1 7 02 


1A S A 

• 3020 


S A 

.20 


si 


A A A A 
• 0004 


A A 14» 

• 0030 


A3 A A 
• 0200 


ATCA 
• 0 7 34 


4 AO 4 
. 10*1 


1 A AO 
• 304* 


s « 
.21 




A A AC 

• 0009 


A A AC 
• 00*3 


A3 A A 
. 0240 


AAC3 
• 003c 


3A 4 A 

• 20 1 4 


1A A* 4 
• 3001 


S3 

.22 


si 
••J 


A A AT 
• 000 f 


A ACT 

• 0037 


A3AC 
.0203 


AOCA 

• 0*34 


S 4 1 A 

• 2130 


lACt 

. 3030 


91 
• 23 


•J? 


A A 1 A 

• 0010 


A A 7 A 
• 00 70 


• 0333 


1 A A 

• 1 000 


SS1A 

• 2230 


1A IT 

.3037 


si 
• 24 




.0012 


A A AT 

• 0007 


aha 

• 030* 


• 1 160 


silt 
• 2336 


IAA1 

. 3003 


• 23 


st 

• CO 


A A 1 Jk 

• ooi o 


A t AC 
• 0103 


AAAA 
• 044* 


1 37A 

• 1270 


• 2424 


SQCT 

.2*37 


3t 

• 20 


• 2 7 


AA3A 

• 0020 


A 1 3T 

.0127 


AC 1 V 

• 0313 


4 111 

• 1 300 


• 2302 


SAOO 

• 20** 


• 27 


• CO 


• 0023 


A 1 C 1 
• 0131 


ACA T 

• 0303 


1 AOO 
• 14** 


SCtO 

• 2301 


SA 1 4 

• 2031 


• 20 


• 2* 


A A 1 1 

• 0031 


• 0 1 7* 


• 0637 


• 1600 


st s a 
.2624 


S* 1 M 

• 2734 


.2* 






A S 4 A 

• 0210 


AT 1C 

• 0733 


• 1713 


.2660 


St t A 

• 2660 


1 A 

.30 


14 


A A A 7 
• 004 7 


A S AC 

.0243 


A A 1 A 

• 0010 


4 A3A 
• 1020 


ST A 4 

.2701 


.2370 


.31 




AACf 

• 0057 


• 0204 


AAA m 

.0*04 


• 1921 


.2721 


s a V A 

• 2476 


.32 




.006* 


.0326 


AAA M 

• 0**4 


SA 4 V 

.20 17 


S W • 

.2731 


.2376 


.33 


1 


A A A 3 

• 000c 


A1T1 

• 0373 


4 A AC 

• 100D 


41 AA 

• 210* 


• 272* 


.2270 


1 A 

• 34 


1 35 


• 0098 


• 0424 


• 1161 


• 2194 


.2716 


.2162 


.35 




• 0116 


• 0479 


• 1276 


• 2272 


• 2693 


• 2052 


• 36 




• 0136 


.0539 


• 1376 


• 2344 


• 2660 


.1941 


• 37 




• 0158 


• 0603 


.1475 


• 2407 


• 2616 


• 1631 


• 36 




• 0184 


• 0671 


.1574 


•2462 


• 2567 


.1721 


• 39 




• 0212 


• 0743 


• 1672 


• 2506 


• 2506 


.1612 


• 40 




• 0244 


• 0619 


• 1769 


• 2545 


• 2442 


• 1506 


• 41 




.0279 


.0900 


• 1663 


• 2573 


.2369 


.1402 


• 42 




• 0316 


• 0963 


• 1955 


• 2592 


• 2291 


• 1301 


• 43 




• 0360 


• 1070 


• 2044 


• 2601 


• 2207 


• 1204 


.44 




• 0407 


• 1160 


• 2126 


• 2600 


• 2119 


.1110 


.45 




• 0456 


• 1253 


• 2207 


• 2590 


• 2027 


.1020 


• 46 




• 0512 


• 1346 


• 2260 


.2571 


.1933 


.0934 


.47 




.0571 


.1445 


• 2347 


.2543 


.1837 


.0853 


.46 




• 0635 


.1542 


• 2406 


• 2506 


.1739 


• 0776 


.49 




• 0703 


.1641 


• 2461 


.2461 


.1641 


.0703 


• 50 



Digitized by 



TABLE A-l. (CONTINUED) 





n - Q 


n - Q 
n — it 




n = 10 


n = in 


n = 10 




p 


r = 1 


r = 0 


r = 9 


r = 8 


r = 7 


r = 6 


P 


• 01 


• 0630 


• 9135 


• 0000 


• 0000 


• 0000 


• 0000 


• 01 


• 02 


• 1531 


.6337 


• 0000 


• 0000 


• 0000 


• 0000 


• 02 


• 03 


• 2116 


• 7602 


• 0000 


• 0000 


• 0000 


• 0000 


• 03 


• 04 


.2597 


• 6925 


• 0000 


• 0000 


• 0000 


• 0000 


• 04 


• 05 


.2965 


• 6302 


• 0000 


• 0000 


• 0000 


• 0000 


• 05 


• 06 


• 3292 


.5730 


• 0000 


• 0000 


• 0000 


.0000 


• 06 


• 07 


.3525 


• 5204 


• 0000 


• 0000 


• 0000 


• 0000 


• 07 


• 00 


• 3695 


• 4722 


• 0000 


• 0000 


• 0000 


.0000 


• 06 


• 09 


• 3609 


• 4279 


• 0000 


• 0000 


• 0000 


• 0001 


.09 


• 10 


• 3674 


• 3674 


• 0000 


• 0000 


• 0000 


• 0001 


• 10 


• 1 1 


.3697 


• 3504 


• 0000 


• 0000 


• 0000 


• 0002 


• 11 


• 12 


• 3664 


• 3165 


• 0000 


• 0000 


• 0000 


• 0004 


• 12 


• 13 


• 3640 


• 2655 


• 0000 


• 0000 


• 0000 


.0006 


• 13 


• 14 


.3770 


.2573 


• 0000 


• 0000 


• 0001 


.0009 


• 14 


• 15 


• 3679 


• 2316 


• 0000 


• 0000 


• 0001 


• 0012 


• 15 


• 16 


• 3569 


• 2062 


• 0000 


• 0000 


• 0002 


.0016 


.16 


• 17 


• 3446 


• 1669 


• 0000 


• 0000 


• 0003 


• 0024 


.17 


• 16 


• 3312 


• 1676 


• 0000 


• 0000 


• 0004 


.0032 


• 16 


• 19 


• 3169 


• 1501 


.0000 


• 0001 


• 0006 


• 0043 


.19 


•20 


.3020 


• 1342 


• 0000 


• 0001 


• 0006 


• 0055 


• 20 


•21 


• 2667 


.1199 


• 0000 


• 0001 


• 0011 


.0070 


• 21 


• 22 


• 2713 


• 1069 


• 0000 


• 0002 


• 0014 


• 0066 


• 22 


• 23 


• 2556 


• 0952 


• 0000 


• 0002 


• 0019 


.0109 


• 23 


• 24 


• 2404 


• 0646 


• 0000 


• 0003 


• 0024 


• 0134 


• 24 


• 25 


.2253 


• 0751 


• 0000 


• 0004 


• 0031 


• 0162 


•25 


• 26 


• 2104 


• 0665 


• 0000 


• 0005 


• 0039 


• 0195 


• 26 


• 27 


• 1960 


• 0569 


• 0001 


• 0007 


• 0049 


• 0231 


.27 


• 26 


• 1620 


• 0520 


• 0001 


• 0009 


• 0060 


• 0272 


• 26 


• 29 


• 1665 


• 0456 


• 0001 


• 0011 


• 0074 


• 0317 


.29 


• 30 


• 1556 


• 0404 


• 0001 


• 0014 


• 0090 


.0366 


• 30 


• 31 


• 1433 


• 0355 


• 0002 


• 0016 


• 0106 


• 0422 


• 31 


• 32 


.1317 


• 0311 


• 0002 


• 0023 


• 0130 


• 0462 


• 32 


• 33 


• 1206 


• 0272 


• 0003 


• 0026 


• 0154 


.0547 


• 33 


• 34 


• 1102 


• 0236 


• 0004 


• 0035 


• 0161 


• 0616 


• 34 




• 1004 


• 0207 


• 0005 


• 0043 


• 0212 


.0669 


• 35 


'.11 


• 0912 


• 0160 


• 0006 


• 0052 


• 0247 


.0767 


• 36 




• 0626 


• 0156 


• 0006 


• 0063 


• 0265 


.0649 


• 37 




• 0747 


.0135 


• 0010 


.0075 


• 0327 


• 0934 


• 36 




• 0673 


.0117 


• 0013 


• 0090 


.0374 


• 1023 


• 39 




• 0605 


• 0101 


• 0016 


• 0106 


• 0425 


• 1115 


• 40 




• 0542 


• 0067 


• 0019 


• 0125 


• 0460 


.1209 


• 41 




• 0464 


• 0074 


• 0024 


.0147 


• 0540 


• 1304 


• 42 




• 0431 


• 0064 


• 0029 


• 0171 


• 0604 


• 1401 


• 43 




• 0363 


• 0054 


• 0035 


• 0196 


• 0673 


• 1499 


• 44 




,0339 


• 0046 


• 0042 


• 0229 


.0746 


.1596 


• 45 




• 0299 


• 0039 


• 0050 


• 0263 


• 0624 


• 1692 


• 46 




• 0263 


• 0033 


• 0059 


• 0301 


.0905 


• 1766 


.47 




• 0231 


• 0026 


• 0070 


• 0343 


.0991 


.1676 


• 46 




• 0202 


• 0023 


• 0063 


• 0369 


• 1060 


• 1966 


• 49 




• 0176 


• 0020 


• 0096 


.0439 


• 1172 


• 2051 


• 50 



Digitized by 



TABLE A-l. (CONCLUDED) 



p 


n = 10 


n » 10 


n «= 10 


n b 10 


n = 10 


n * 10 


D 

r 


* 


r = 5 


r « 4 


r * 3 


r = 2 


r = 1 


r = 0 


• 01 


• 0000 


• 0000 


• 0001 


• 0042 


• 0914 


• 9044 


• 01 


• 02 


• 0000 


• 0000 


• 0006 


• 0153 


• 1667 


• 6171 


• 02 


• 03 


• 0000 


• 0001 


• 0026 


• 0317 


• 2261 


• 7374 


• 03 


• 04 


• 0000 


• 0004 


• 0056 


• 0519 


• 2770 


• 6646 


• 04 


• OS 


• 0001 


• 0010 


• 0105 


• 0746 


• 3131 


• 5967 


• 05 


• 06 


• 0001 


• 0019 


• 0166 


• 0966 


• 3436 


• 5366 


• 06 


• 07 


• 0003 


• 0033 


• 0246 


• 1234 


• 3643 


• 4640 


• 07 


• 09 


• 0005 


• 0052 


.0343 


• 1476 


• 3777 


• 4344 


• 06 


• 09 


• 0009 


• 0076 


• 0452 


• 1714 


• 3651 


• 3694 


• 09 


• 10 


• 0015 


• 0112 


.0574 


• 1937 


• 3674 


• 3467 


• 10 


• 11 


• 0023 


• 0153 


• 0706 


• 2143 


• 3654 
.3796 


.3116 


.11 


• 12 


• 0033 


• 0202 


• 0647 


• 2330 


.2765 


. 12 


• 13 


• 0047 


• 0260 


.0995 


.2496 


• 3712 


• 2464 


.13 


• M 


• 0064 


• 0326 


• 1146 


• 2639 


.3603 


.2213 


.14 


• 15 


• 0065 


• 0401 


• 1296 


.2759 


.3474 


.1969 


.15 


• 16 


• 0111 


.0463 


• 1450 


• 2656 


• 3331 


.1749 


. 16 


• 17 


• 0141 


.0573 


• 1600 


.2929 


.3176 


.1552 


• 17 


• 16 


.0177 


.0670 


.1745 


.2960 


• 3017 

• 2652 


.1374 


• 16 


• 19 


• 0216 


.0773 


• 1663 


• 3010 


• 1216 


• 19 


• 20 


• 0264 


• 0661 


• 2013 


• 3020 


.2664 


.1074 


• 20 


• 21 


• 0317 


.0993 


• 2134 


• 3011 


• 2517 


.0947 


• 21 


• 22 


• 0375 


• 1106 


•2244 


• 2964 


• 2351 


• 0634 


• 22 


• 23 


• 0439 


• 1225 


•2343 


• 2942 


• 2166 


.0733 


• 23 


• 24 


• 0509 


• 1343 


• 2429 


• 2665 


• 2030 


• 0643 


• 24 


•25 


• 0564 


• 1460 


• 2503 


• 2616 


.1677 


• 0563 


• 25 


• 26 


• 0664 


• 1576 


• 2563 


• 2735 


• 1730 


• 0492 


• 26 


•27 


• 0750 


• 1669 


• 2609 


• 2646 


• 1590 


• 0430 


• 27 


• 26 


• 0639 


.1796 


• 2642 


• 2546 


• 1456 


.0374 


• 26 


•29 


• 0933 


• 1903 


• 2662 


• 2444 


• 1330 


.0326 


• 29 


• 30 


• 1029 


•2001 


• 2666 


• 2335 


• 1211 


.0262 


• 30 


• 31 


• 1126 


• 2093 


• 2662 


• 2222 


.1099 


• 0245 


• 31 


• 32 


• 1229 


• 2177 


• 2644 


• 2107 


.0995 


• 0211 


• 32 


• 33 


• 1332 


• 2253 


• 2614 


.1990 


• 0696 


• 0162 


• 33 


• 34 


• 1434 


.2320 


• 2573 


.1673 


• 0606 


• 0157 


• 36 


•JJ 


• 1536 


.2377 


• 2522 


.1757 


• 0725 


.0135 






• 1636 


• 2424 


• 2462 


• 1642 


• 0649 


• 0115 






• 1734 


• 2461 


• 2394 


• 1529 


• 0576 


.0096 






• 1629 


• 2467 


• 2319 


• 1419 


• 0514 


• 0064 






• 1920 


.2503 


• 2237 


• 1312 


• 0456 


.0071 






• 2007 


• 2506 


• 2150 


• 1209 


• 0403 


• 0060 






•2067 


•2503 


• 2056 


• 1111 


• 0355 


• 0051 






•2162 


•2466 


• 1963 


• 1017 


.0312 


• 0043 






•2229 


• 2462 


• 1665 


• 0927 


.0273 


• 0036 






•2269 


• 2427 


,1765 


• 0643 


• 0236 


• 0030 






• 2340 


•2364 


• 1665 


• 0763 


• 0207 


• 0025 






• 2363 


• 2331 


• 1564 


• 0666 


• 0160 


• 0021 






• 2417 


• 2271 


• 1464 


• 0619 


• 0155 


• 0017 






•2441 


• 2204 


• 1364 


• 0554 


• 0133 


• 0014 






• 2456 


.2130 


.1267 


.0494 


• 0114 


• 0012 






• 2461 


.2051 


• 1172 


• 0439 


• 0096 


• 0010 





Digitized by 



TABLE A- 2. PARTIAL SUMS 



PROBABILITY OF REALIZING s > r SUCCESSES IN n TRIALS 





PCs > 


r/n, p) 

Mr J 


n 


n 1 
11 • 


s n-s 




LmJ 

O S T* 

s = r 


s!(n - 


S) . 




p 


n = £ 


_ _ o 
n = I 


n - o 


n - o 


n = 3 


n = 4 


p 


r = 2 


r = 1 


r = 3 


r = 2 


r = 1 


r = 4 


r 


• 01 


• 0001 


• 0199 


• 0000 


• 0003 


• Oct? 


• 0000 


• 01 


• 02 


• 0004 


• 0396 


• 0000 


• 0012 


• 0506 


• 0000 


• 02 


• 03 


• 0009 


• 0591 


• 0000 


• 0026 


• 00 7 3 


• 0000 


• 03 


• 04 


• 0016 


• 0764 


• 0001 


• 0047 


• 1 153 


• 0000 


.04 


• 03 


• 0025 


• 0975 


• 0001 


• 0073 


• 1426 


• 0000 


• OS 


• 06 


• 0036 


• 1 164 


• 0002 


• 0104 


• 16t4 


• 0000 


• 06 


• 07 


• 0049 


• 1351 


• 0003 


• 0140 


« All 

• It56 


• 0000 


• 07 


• 06 


• 0064 


• 1536 


• 0005 


• 0162 


• 2213 


• 0000 


• 06 


• 09 


• 0061 


• 1719 


• 0007 


• 0226 


• 2464 


• 0001 


• 09 


• 10 


• 0100 


• 1900 


• 0010 


• 0260 


• 27 10 


• 0001 


• 10 


• 1 1 


• 0121 


• 2079 


• 0013 


• 0336 


• 2"50 


.0001 


• 11 


• 12 


• 0144 


• 2256 


• 0017 


• 0397 


• 3105 


• 0002 


• 12 


• 13 


• 0169 


• 2431 


• 0022 


• 0463 


• 3415 


.0003 


• 13 


• 14 


• 0196 


• 2604 


• 0027 


• 0533 


• 363" 


• 0004 


• 14 


• 13 


• 0225 


• 2775 


• 0034 


• 0607 


V ACQ 

• 305V 


• 0005 


• 15 


• 16 


• 0256 


•2944 


• 0041 


• 0666 


• A 9 Y 

• 407 3 


• 0007 


• 16 


• 17 


• 0269 


• 3111 


• 0049 


• 0769 


A ^ A ^ 

• 4262 


• 0006 


• 17 


• 16 


• 0324 


• 3276 


• 0056 


• 0655 


.4406 


• 0010 


• 16 


• 19 


• 0361 


• 3439 


• 0069 


• 0946 


• 4606 


• 0013 


• 19 


• 20 


• 0400 


• 3600 


• 0060 


• 1040 


« A A A 

• 4000 


• 0016 


• 20 


•21 


• 0441 


• 3759 


• 0093 


• 1136 


• 5070 


• 0019 


• 21 


•22 


• 0464 


• 3916 


• 0106 


• 1239 


• 5254 


• 0023 


• 22 


• 23 


• 0529 


• 4071 


• 0122 


.1344 


Bate 

• 5435 


.0026 


• 23 


•24 


• 0576 


• 4224 


• 0136 


• 1452 


• 5610 


.0033 


• 24 


•25 


• 0625 


• 4375 


• 0156 


• 1562 


A "T A « 

• 5701 


• 0039 


• 25 


• 26 


• 0676 


• 4524 


• 0176 


• 1676 


e A • a 

• 5940 


• 0046 


• 26 


•27 


• 0729 


• 4671 


.0197 


• 1793 


• 61 10 


• 0053 


• 27 


•26 


• 0764 


• 4616 


• 0220 


• 1913 


• 6260 


• 0061 


• 26 


•29 


• 0641 


• 4959 


• 0244 


• 2035 


• 6421 


• 0071 


• 29 


• 30 


• 0900 


• 5100 


• 0270 


• 2160 


• 6570 


• 0061 


ill 


• 31 


• 0961 


• 5239 


• 0296 


• 2267 


• 6715 


.0092 




• 32 


• 1024 


• 5376 


• 0326 


• 2417 


• 6056 


• 0105 




• 33 


• 1069 


• 551 1 


• 0359 


• 2546 


• 6??2 


• 0119 




• 34 


• 1156 


• 5644 


• 0393 


• 2662 


.7125 


• 0134 






• 1225 




• 0429 


• 2617 
.2955 


v act 
• 725* 


.0150 




ill 


• 1296 


'Mil 


• 0467 


.7379 


• 0166 






• 1369 




• 0507 


• 3094 


• 7300 


.0167 






• 1444 
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• 1143 


• 26 


• 27 
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APPENDIX B 
TECHNIQUES FR SIMULATION TESTNG 



B- 1 • Simulation techniques. 

a. The catalog of techniques that follows was 
compiled primarily from two sources. "Nuclear 
Blast and Shock Simulators/ a report 
(GE-TEMPO, 1972) by Panel N-2 of the Tripartite 
Technical Cooperation Programme (whose member 
nations are the United States, Canada, and the 
United Kingdom), contains descriptions of tests 
from which summaries have been tabulated herein. 
The second source was Rowan (1974) which tab- 
ulates other items besides airblast and ground mo- 
tion that are common to the TTCP report. Specifi- 
cally, there are sections relating to simulation of 
nuclear radiation, spallation impulse, thermal radi- 
ation, EMP, and debris. While a bit dated in some 
cases, it is still a worthwhile compilation from 
which future work may proceed. Additional re- 
ferences are listed in the bibliography, appendix D, 
grouped into categories pertinent to testing for nu- 
clear weapon effects. 

b. An extensive catalog of techniques is presented 
in tables B-1 through B-10. Recommended test con- 
cepts range from laboratory and simple in-plant 
tests to numberous field tests that can stimulate 
all or portions of the nuclear environment exper- 
ienced by major components and subsystems of a 
facility. Laboratory and in plant test concepts. 
However, experience with field tests is limited; the 
cost of field tests will generally be high. 

c. The important test comcepts considered are as 
follows: 

— High-Explosive Simulation Technique 

(HEST) 

— Direct-Induced High-Explosive Simulation 

Technique (DIEHEST) 

— High-Explosive Contact Surface Burst 

— Underground Nuclear Tamped Burst 

— Underground Nuclear Tunnel Test 

— Giant Reusable Air-Blast Simulator 

(GRABS) 

— EMP Simulation Testing 

— Blast Simulation Technique for Testing Air 

Entrainment Systems and Blast 
Closures. 

The following subsections discuss these critical test 
concepts in greater detail than is possible in the 
catalog presented as tables. More detailed evalu- 
ation and development of each of these test con- 
cepts is necessary in order to select the appropriate 
hardness verification test procedures for a pos- 



tulated nuclear burst. 

B-2. High-explosive simulation 
technique 

a. The HEST concept has been used four times 
for tests on operational Minuteman sites. The tech- 
nique had been expanded for high-overpressure 
tests in a rock medium. Periodically, variations of 
the HEST test are utilized for special applications. 
The reader is urged to review the current literature 
to determine whether recent tests have relevance 
to his application. A summary of HEST experi- 
ments from 1964 to 1974 is listed in table B-ll. 

b. The objective of the technique is to simulate 
the overpressure and superseismic air-induced 
ground shock from a nuclear detonation 
Operational and small scale tests have 
demonstrated the feasibility of simulating over- 
pressures (for about the first 200 msec) from yields 
up to 10 Mt and for overpressures up to 3000 psi. 

c. HEST uses a confined detonation of explosive 
fuse (Primacord) to produce a pressure pulse de- 
signed to travel over the ground surface at the 
same velocity as an air shock wave of equal inten- 
sity, and to have a timedecay shape similar to the 
early part of the pressure pulse produced by a 
nuclear detonation. A HEST facility consists of a 
platform structure constructed above the surface of 
the ground over the installation to be tested. (See 
figure 4-1 for facility configuaration.) The platform 
supports an overburden of earth and forms a cavity 
between the bottom of the overburden and the 
ground. An earthen embankment is built around 
the perimeter of the platform. Primacord is 
wrapped on wooden racks that are suspended in 
the cavity. Since Primacord detonates at a velocity 
greater than the shock front velocity to be simu- 
lated, the cord is wrapped at an angle to the direc- 
tion of propagation. The intensity of the pressure 
pulse depends primarily on the loading density 
(amount of explosives per unit cavity volume). The 
overburden placed over the HEST platform is com- 
pressed and accelerated upward when the overpres- 
sure acts on its base, and this upward motion of 
the overburden causes the volume of the cavity to 
expand, with a corresponding decrease in pressure. 

d. Operational HEST tests have used a cavity 
that is normally 300 ft. sq. For HEST III the 
cavity was 5.5 ft deep and the overburden was 
about 10 ft thick. Primacord charge density was 
0.070 lb of explosive per cubic foot of cavity, and 
the Primacord weave angle was 8 deg 27 min. 
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Table B-l. Dynamic Pressure Tests Using Shock Tubes to Simulate AirbUut Effect 



TABLE B-l. DYNAMIC PRESSURE TESTS USING SHOCK TUBES TO SIMULATE 
AIRBLAST EFFECT 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


DASACON 

U.S. Naval Wea- 
pons Laboratory 
Dahlgren, VA 


Conical high- 
explosive- driven 
shock tube. 
Probably the 
world's largest 
(2455 ft long). 


Test section 7.5. ft 
x 10 ft dia., 5.5 bar 
(80 psi), 300 ms. 

Test section 9 ft 

x 15 ft dia., 2.7 bar 

(40 psi), 500 ms. 

Test section 12 ft 
x 22 ft dia., 1.4 bar 
(20 psi), 500 ms. 


Air-blast loading 
on full-scale and 
small models. For 
model sections of 
10-ft dia. and 
below, altitude 
simulation up to 
10 5 ft is available. 


Time dependent 
measurements of: 
Overpressure 
Dynamic Pressure 
Acceleration 
Velocity/ 

Displacement 
Strain 
Video camera and 
recorder. 

The above results 
are available as 
analog or 
magnetic tape. 
High-speed 
cine camera film 
(up to 10,000 
rps). 


URS Shock Tunnel 
San Francisco, 
CA 


Rectangular ex- 
plosively-driven 
reinforced con- 
crete shock tube. 
Expansion 
chamber 8.5 ft 
x 12 ft x 92 ft 
long. 

Driver chamber 
8.5 ft x 8 ft x 
63 ft 


Test section 8.5 ft 
by 12 ft dia., 40 ft 
long. Overpressure 
range 0.1 to 0.75 bar 
(1.5 to 11 psi). 
Duration range: 
80 to 50 ms. 


Frontal air-blast 
loading on struc- 
tural elements. 


Time dependent 
measurements of: 
Overpressure 
Dynamic Pressure 

Results on 
magnetic tape 
and high-speed 
cine cameras 1000 
fps color, 2000 
fps black and 
white. 


BRL, Aberdeen 
Proving Ground, 
MD 


84 ft x 500 ft 
total length 
Compressed air- 
driven dual shock 
tube. Can be usee 
either (a) as a 
separate shock 
tube or (b) with 
the expansion 
chamber joined to 
the 5.5 ft dia. 
shock tube (see 
below), joined 
through a large 
air breathing en- 
gine as a dual 
shock tube to 
test the behavior 
of such engines 
under shock 
loading. 


Test section 8 ft dia. 
No target fixing 
facilities— each 
target to be secured 
to test section walls. 
Overpressure range: 
0.3 to 1.8 bars 
(3 to 26 psi). 
1 sec approx. 


Air-blast loading 
on full-scale and 
small models. 
Simultaneous 
internal blast 
loading from in- 
let to outlet of 
air breathing en- 
gines 


Time dependent 
measurements of: 
Overpressure 
Dynamic Pressure 

No high-speed 
cine facilities. 


BRL, Aberdeen 
Proving Ground, 
MD 


5.5-ft-dia.x 610ft 
air-driven shock 
tube. Driven gas 
can be heated to 
eliminate density 
change on expan- 
sion of driven 
air. 


Test section 5.5 ft 
dia. Remaining 
details as for the 
8 ft dia. tubes. 


As above. 


As above. 
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TABLE B-l. (CONTINUED) 



FACILITY 


EQUIPMENT 


pr.RroRMANcr. 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


AFWL 

Kirtland AFB, 
NM 


6 ft dia. x 245 ft 
cylindrical high- 
explosive-driven 
shock tube. 


Test section flat 
platform 5 ft x 8 ft 
with 4.5 ft headroom. 
Overpressure range 
up to 7 bars (100 psi) 
Duration: 100 to 200 
ms. 


Ground shock and 
related studies. 
Air-blast loading. 


As above. 


Sandia Labs 

Albuquerque, 

NM 


19-ft-dia. high- 
explosive-driven 
blast tunnel. 
6 ft dia. x 50 ft 
driver. 
6 ft x 150 ft 
expansion. 40 ft 
divergent section. 
19 ft dia. x 66 ft 
expansion section. 


Test section 66 ft x 
19 ft dia. Initial 
pressure in the test 
section can be varied 
from 1 psia to 12 
psia. At 12 psia, 
static overpressure 
varies between 1 to 
6 bars (14 to 87 psi). 
Duration: 25 to 
50 ms. 


Large military 
units at large 
angles of attack. 


Time -dependent 
measurements of: 

Overpressure 

Dynamic Pressure 

Stagnation Press, 

Acceleration 

Displacement 

Strain 

Velocity 
Results on 
magnetic tape. 

High-speed cine 
camera film up to 
8,000 fps. 


januia LtaDS 

Albuquerque, 

m 


A_f t-.fii a y ?nn. 
ft high explosive- 
driven shock 
tube. 


target suspended by 
straps, released 
just before shock wave 
arrival, and loaded in 
free fall. Test sec- 
tion ambient pres- 
sures from 0.013 bar 
(0.2 psia). At ambient 
of 0.8 bars, overpres- 
sure range is 6 to 14 
bars (87 to 200 psi). 
Duration: ^ 6 ms. 


Ciind.il nil 11 ceil y 

structures 


i\s ior UilC 

19 ft dia. 




6-ft-dia. x 50- 
ft-high-explosive- 
driven shock 
tube. Expendable 
driven section. 


Test section 6-ft- 
dia. Target ac- 
commodation as for 
200-ft tube above. 
Test section am- 
bient pressure, 0.07 
bar, overpressure 
10 bars (145 psi). 
Duration: 2 to 8 ms. 


As above. 


As for 19 ft dia. 


Lovelace DNA 
Kirtland AFB, 
Albuquerque, 
NM 


Cylindrical com- 
pres sed- air-driv- 
en shock tube. 
3.5 ft x 15 ft 
compression 
chamber, 3.5 ft 
x 125 ft exp. 
chamber. 6 ft 
x 30 ft exp. 
chamber. 9 ft 
long coupling 
sect. 


Test section (s) 65 ft 
by 3.5 ft dia. Over- 
pressure: 0 to 2.3 bars 
(0 to 33 psi) . Duration 
100 ms on end plate 
terminating the tube. 


Effect of long- 
duration airblast 
on a variety of 
large animal 
species. 


Time dependent 
measurements of: 
Overpressure 
Dynamic Pressure 
Photography 
Oscillograph and 
Chart Recorders 
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TABLE B-l. (CONCLUDED) 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


DRES 
Ralston, 
Alberta, Canada 


6-ft-dia. high- 
explosive-driven 
shock tube, 
i/riveri irom t\UA/ 
TNT in a gun bar- 
rel, or Primacord 
line or grid 
charge at entrance 
to expansion 
chamber. 

16 in. naval gun 
with 165 in. long x 
6 ft dia. expans. 
chamber. 57 ft x 
16 in. permits re- 
coil less operation, 


Test sections all 
6 ft dia. 

Overpressure: 0.1 to 

"7 A hare M A t>r» 

psi). Duration: 50 
to 70 ms. 


Blast loading. 

A. High pressure 
structural response 
st ud ies . 

B. Soil medium 
studies. 

C. Tests on mili- 
tary equipment and 
material . 

D. Tests on targets 
in normal reflect i or 
mode. 


Time dependent 

measurements of: 
Overpressure 
uynainic rressure 
Acceleration 
Velocity 
Displacement 
Strain 

High-Speed Cine 

(20,000 fps) 

17-in. 

Shadowgraph 

Video camera 

and Recorder 


DRES 
Ralston, 
Alberta, Canada 


3-ft-dia. com- 
pressed-air or 
high- explosive- 
driven shock tube. 
Driver section 
similar to 6 ft 
dia. model, but 
uses a 14-in 
naval canon barrel 
Approx 175 ft. 


Test section 6 f t x 
3- ft dia. 

Compressed Air Driven. 
Overpressure 0.5 to 
1.0 bars 
(2 to 14 psi). 
Duration: 200 ms. 
RDX/TNT Driven. Over- 
. pressure 1.1 to 7.0 
bars (15 to 100 psi). 
Duration: 50 to 100 ms. 


Tests on military 
equipment and 
material and 
dynamic loading 
effects. 


Time dependent 

measurements of: 
Overpressure 
Dynamic Pressure 
Acceleration 
Velocity 
Displacement 
Strain 

High-Speed Cine 

(20,000 fps) 

17 in. Shadowgraph 

Video camera 

and Recorder 


AFWL 

Kirtland AFB, 
Albuquerque, NM 


13-in. dia. driven 
by explosive hy- 
drogen and oxygen 
gas mixture; 
Approx 220 ft over 
all length. 


Test section approx 2 ft 
x 13-in. dia. 
Overpressure range 
3.5 to 45 bars (50 to 
650 psi). Duration: 
1 to 4 ms. 


High- pressure 
tests of Air Force 
targets. 


Pressure-Tine 
oscillographs 


BRL 

Aberdeen Proving 
Ground, MD 


8- in. dia. detona- 
tion-driven shock 
tube. Hydrogen - 
oxygen driver 
mixture. 


Test section 55 ft x 
22 -in. dia. Over- 
pressure 14 bars 
(200 psi) . Duration: 
50 ms. 


High-pressure 
loading. 


Oscillographs, 
chart, and tape 
recorders . 


Boeing/SAMSO, 
Boeing Tulalip 
Test Site, 
Tulalip, Wash. 


24 in. , 36 in. , and 
30 in. x 52 in. 
tubes of undeter- 
mined length 


No data at hand. 






NOL 

Silver Springs, 
Maryland 


Conical shock 
tube. 


16-in. dia. test 
section at 75 ft; 
24 -in. dia. test 
section at 135 ft. 
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Table B-2. Ground Shock Using Blast Load Generators to Simulate Airblast Effect 



TABLE B-2. GROUND SHOCK USING BLAST LOAD GENERATORS TO SIMULATE 
AIRBLAST EFFECT 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


WES 

Vicksburg, MS 


Large blast 
load generator 
(LBLG) . 


Soil loading test 
section surface area 
410 ft 2 . Maximum 
available depth of 
burial in ft. Pres- 
tuin* to 70 hat*^ 

(1000 psi) . Dura- 
tion: to 2 sec. 


Tests for studies 
in the design and 
analysis of under- 
ground structures. 


Time dependent 
measurements of: 
Overpressure 
Soil Stress 
Acceleration 
Vel oci ty 
Displacement 
St r. tin 
Magnetic Tape 
Records. 




4-ft dia. 
blast load 
generator. 


Static pressures 

up to 135 bars (2000 

psi) of compressed 

air or water. 

Dynamic pressures 

to 17 bars (250 psi) . 

Durations: to 2 sec. 


Testing small 
buried structures, 
stress wave prop- 
agation, and soil 
structure 
interaction. 


As above. 




Static loading 
device. 


Static pressure to 
410 bars (6000 psi) . 


Static loading of 
buried model 
structures, 
structural ele- 
ments, and sub- 
merged objects. 


Static measure- 
ments of: 
Fluid Pressure 
Soil Stress 
Displacement 
and Strain 


Waterways Exper- 
imental Station 
Vicksburg, MS 


Vertical 
detonatable gas 
shock tube. 


Peak pressures 55 to 
105 bars (800 to 1500 
psi) . Durations 3 to 
6 ms. Test Cham- 
bers: 46.7 in. dia. 
Depth: (a) 1.75 ft. 

(b) 4 ft. 

(c) 5.75 ft. 


Shot duration 
pressure pulses on 
buried model struc- 
tures and structural 
elements. 


Time- dependent 
measurements of: 

Pressure j 

Soil Stress 

Acceleration 

Strain 
Magnetic Tape 
Records. 


NCEL 

Port Hueneme, 
CA 


Blast simu- 
lator, Primacord 
driven . 


Test pit: 9 ft wide x 
10 ft long x 12 ft 
deep. Peak Pressures: 
2 to 14 bars (30 to 
200 psi) . 

Duration: 0.4 to 7.0 
sec. 


Static and dynamic 
loads on structural 
elements such as 
beams, slabs, and 
model elements. 


Time-dependent 
measurements of: 
Pressure and 
Strain 
Oscillograph 
Records. 


AFWL 

Kirtland AFB 
Albuquerque, NM 


2-ft-dia. 
(0.6 m) high- 
explosive-driven 
vertical shock 
tubes . 


Test pit: 4 ft deep x 
2 ft dia. Surface 
loading. Peak Over- 
pressure 0.35 to 35 
bars (5 to 500 psi) . 
Duration: 30 ms. 


Instrumentation 
development, proof, 
and testing 


Time-dependent 
measurements of: 
Soil 

Pressure 


GRABS 

Air Force Weap- 
ons Laboratory 
Kirtland AFB, 
NM 


Explosive cavity 
above soil sample 
in large 
reinforced 
concrete- lined 
silo. 


18-ft-dia. x 48-ft- 
deep silo, soil test 
bed up to 30 ft of 
depth possible. 
Up to 1800 psi over- 
pressure possible. 


Buried systems or 
models of systems, 
and soil/structure 
interactions. 


See 

for discussion 
of GRABS 
technique. 
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Table B-S. Dynamic Loading of Material Using Special Test Machines to Simulate Airblast Effect 



TABLE B-3. DYNAMIC LOADING OF MATERIAL USING SPECIAL TEST MACHINES 
TO SIMULATE AIRBLAST EFFECT 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


Wyle Laboratories 
Norco, Calif. 


Parallel Pendulum 
Impact . 


Velocity peak 
400 in. /sec. 


Shock testing 
pieces of equipment 
or of shock- 
mounted systems 






Hammer 
shock machine 


Velocity peak 
200 in. /sec. 


As above . 




Kirtland AFB 
Albuquerque, NM 


HYGE shock 
tester 


200 g max 

Duration: 10-50 ms. 
Peak force: 47,000 lb 
Wave shape: 1/2 sine 
square, triangle 


As above. - 


-- 


Anywhere 


Quick -re lease 
twang test, from 
small lab tests 
to full-scale 
simulation 


Simulates motion 
by displacement of 
platform rather than 
supports . 


Shock-isolated 
platform evaluation 
tests. 

Very effective when 
properly conducted. 




WES 

Vicksburg, MS 


500-kip ram 
loader. 


Hydraulic actuator, 
pump, and control 
system can apply 
loads in excess of 
350,000 lb with 
approx 1/4- in move- 
ment and 80 ms rise- 
time. 


Applies loads to 
construction mate- 
rials and struc- 
tural elements 
where it is neces- 
sary to achieve 
slowly applied 
cyclic or random 
loads. 


Time -dependent 
measurements of: 

Load 

Strain 

Acceleration. 

Velocity 

Displacement 
d.c. to 40 kHz FM 
magnetic recorder; 




200-kip dynamic 
ram loader. 


Hydraulic ram 
applies 10,000 to 
20,000 lb in either 
tension or compres- 
sion, either slowly 
or within 2 ms. 


Tests structural 
elements and de- 
termines strength 
of materials under 
slow static or 
dynamic loads. 


Time- dependent 
measurements of: 

Load 

Strain 

Acceleration 
Velocity 
Displacement 
d.c. to 40 kHz FN 
magnetic record- 
ing 


NCEL, 

Port Huenene, CA 


SO-kip ram 
loader. 


Max force: 50,000 lb. 
Rise time: 2 to 
200 ms. Duration: 
0 to 2 sec. Head 
velocity: 3 to 1800 
in./min. 


Dynamic testing of 
metal, concrete, 
and other mate- 
rial specimens. 


Time-dependent 
measurements of: 
Head displace- 
ment 

Strain 

Head resistance 
d.c. to kHz FM 
magnetic recording 


AFWL 

Albuquerque, NM 


50- kip ram 
loader. 


Pneuma t i c - hydr au lie 
driver, 100 to 50,000 
lb. Nominal 59,000- 
lb maximum. Rise 
time: 1.5 to 50 ms. 
Duration: 20 ms to 
30 sec. Velocity: 4 
in./min to 450 in./min. 
Test section: 14-to 
74 -in. wide by 48 in. 
Used either force- 
time or constant 
velocity loader. 


Dynamic testing 
of properties of 
materials. 


Refer to AFWL . 



U.S. Army Corps of Engineers 
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TABLE B-3. (CONCLUDED) 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


WES 

Vicksburg, MS 


Dynamic 
triaxial test 
apparatus . 


Compressed -gas 

vert ical -load 

generator 10- to 

5000-lb load, 4-in. 

stroke, 10 ms to 

1 sec time to failure. 


Dynamic testing to 
determine shear 
strength of soils. 


Time -dependent 
measurements of: 
Top load 
Bottom load 
Top 
displacement 


WES 

Vicksburg, MS 


Dynamic load- 
test systems. 


Compressed-gas 
loaders apply con- 
trolled impulse or 
static loads 100 to 
115,000 lb on test 
specimens up to 
8-ft wide by 50-in 
high. 


To determine re- 
sponse of earth 
materials and soil 
structure systems 
to high amplitude 
impulse loading 
(nuclear blast 
simulation) . 


Time -dependent 
measurements of: 
Force genera- 
tion 
Uniaxial strain 
Vertical fluid 

pressures 
Displacements 
Triaxial chambei 

pressures 
Vertical load 
Vertical and 
horizontal 
displacements 
on triaxial 
specimens. 

Oscillograph and 
ana log- tape 1 
recording. 




Explosive- 
test facility. 


Explosive plane- 
wave lens system 
applies high pres- 
sures (600 kb] 
on up to 6-ft 2 sur- 
faces. Flying plate 
techniques also 
employed. 


j Determination of 
Hugoniot equations- 
of-state on: 
Cements 
Epoxies 
Geological 

material 
Concrete 

(Hardened 
specimens) . 


Time- dependent 
measurements of: 
Pressure 
Strain 
Time of 
arrival 
High-speed 
oscillographs 
and magnetic 
tape recording 
High-speed 
cameras. 


iment Station 
Vicksburg, f6 


Compre s s ed - 
gas and powder 
gun facility. 


Gas gun 20 ft x 2 -in. 
dia. operated by 
high-pressure helium 
or nitrogen. Speci- 
men launched into 
evacuated tube. Pow- 
der gun: 20 ft x 2.5- 
in. dia. Explosive 
driving of specimens 
up to 8000 ft/sec. 
Launch tube evacuated 
and detonation pro- 
ducts contained. 


Shock response or 
Hugoniot data on 
geological mate- 
rials, grouts, and 
concretes . 


Particle velocity 
time and 
pressure time. 

Oscillograph 
recording. 


NCEL 

Port Huenene, CA 


Dynamic 
slab loader. 


Primacord driver. 
Maximum pressure: 
300 psi. Maximum 
force: 1,560,000 lb. 
Rise time: 1 to 3 
ms. Decay (several 
sec) . 


Tests of reinforced 
concrete slabs. 


Oscillograph 
recording . 

Oscilloscope 
photography. 
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Table B-i. Full-Scale Soil Loading Tests ofAirblast Effect During Field Tests 



TABLE B-4. FULL-SCALE SOIL LOADING TESTS OF AIRBLAST EFFECT DURING 
FIELD TESTS 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


HEST 

Air Force Weap- 
ons Laboratory 
Kirtland AFB, 
NM 


Primacord laced 
around wooden 
racks and deton- 
ated in an air 
cavity benearth an 
earth overburden. 


Ground surface air 
overpressures in 
excess of 1000 psi 
possible. Accom- 
modates very large 
test structures. 


Actual prototype, 
buried systems, 
models of systems, 
and soil/structure 
interaction. 


DIHEST 

Air Force Weap- 
ons Laboratory 
Kirtland AFB, 
NM 


Arrays of explo- 
sives embedded in 
the ground. 


Horizontal longi- 
tudinal peak parti- 
cal velocities of 
nearly 100 fps at a 
range of 10 ft possi- 
ble. Accommodates 
full-size test struc- 
tures . 


Buried systems in 
a variety of geo- 
logical formations. 
DIHEST coupled 
with HEST provides 
simulations of 
air-blast- induced 
and direct- 
induced ground 
motion . 


DELTA 

Civil Engineer- 
ing Research 
Facility 
Kirtland AFB, 
NM 


Explosive cavity 
above or below 
test slab inside 
reinforced con- 
crete cylinder. 


30,000-psi maximum, 
13-ft-dia. test sec- 
tion, 2-day assembly 
for each test. 


Models of systems, 
test to failure. 


DRES Portable 
Ralston, Alberta 
Canada 


Horizontal layer 
of explosive with 
water overburden 
can be placed over 
any existing 
buried structure 
to be tested. 


30 to 100 psi over- 
pressure range. 
Positive duration, 
40 ms. Ground 
surface cover is 
17.5 x 16 ft. 


Portable ground- 
blast simulator 
places time- 
dependent over- 
pressure loads on 
existing buried 
targets. 


Blast-Di recting 
Techniques 

DRES 

Ralston, 

Alberta Canada, 

and 

AWRE 

Foulness, 
England 


HE arranged in 
vertical flat 
plate. 


Cost well below 
similar pressure 
resulting from 
hemispherical 
trials. 30 to 450 
psi range. Target 
dimension should 
be smaller than 
the HE reactangle. 


Test of full-scale 
equipment. 



U.S. Army Corps of Engineers 
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Table B-5. Dynamic Water Wave and Shock Tests, Simulating Airblast o Direct Shock 
Table B-6. Blow-Off Response Using Sheet Explosives to Simulate Nuclear Radiation Effect 



TABLE B-5. DYNAMIC WATER WAVE AND SHOCK TESTS, SIMULATING AIRBLAST 

ON DIRECT SHOCK U.S. Army Corps of Engineers 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM 
OF RESULTS 


WES 

Vicksburc MS 


Big black 
river test 
facility. 


Basin 1: 150 ft x 250 
ft x 22 ft deep. 
Detonations up to 150 lb 
HE. 


Underwater explo- 
sion effects . 


Time-dependent 

me asiiTeinents of • 

Pressure 

Surface Displace- 






Basin 2: 160 ft x 260 
ft x 12 ft deep. Deto- 
nation up to several 
hundred pound HE. 

Basin 3: Trapezoidal 
with sloping bottom 
550 ft long; 300 ft wide 
at 0.2 ft depth, 150 ft 
wide at 12 ft depth. 


Explosively gen- 
erated surface 
effects. 


ment 
Bottom Ground 

motion 
High-Speed 
c inephotographs 
up to thousands 
of frames/sec. 


Lovelace DNA 
Kirtland AFB 

A 1 Khoho rniiA KIM 


Water 

shock facility. 


Basin 220 ft x 150 ft 
at top, 30 ft deep over 

a Tfl ft v 1 DO ft nr»r- 
H Ou II a lull II p\Jl~ 

tion. Upper limit of 
detonations 27 lb. 


Biomedical inves- 
tigations of under- 
water Diasi errects . 


Time-dependent 
measurements of: 

r rcasurc 
Oscillograph 

Records. 


NCEL 

Port Hueneme, 
CA 


Impulse 
water wave 
facility. 


Basin 94 ft x 92 ft 
x 3 ft deep. Com- 
pressed- air- dri ven 
plunger. Single 
abrupt and omni- 
directional motions 
of the plunger are 
used. 


Effects of water 
waves on model 
waterfront 
structures . 


Cinephotography . 



TABLE B-6. BLOW-OFF RESPONSE USING SHEET EXPLOSIVES TO SIMULATE 
NUCLEAR RADIATION EFFECT UtS . ^ Covpe of Engineers 



FACILITY 


EQUIPMENT 


PERFORMANCE 


TYPE OF TESTS 


TYPE AND FORM j 
OF RESULTS 


Any facility 
familiar with 
HE techniques 


Explosive sheet 
in contact 
with specimen 


10** taps -0.2 ms 
(time proportional to 
impulse), 100 kb peak 
pressure . 

Maximum achievable is 
probably less than 10 
taps - 10 7 taps gives 
a 200-ms pulse. 


Radiation impulse 
loads of closures, 
exposed structural 
elements. 






Explosive sheet 
over neoprene 
foam 


3000 taps on up, 
2000-2 ms 
0.1 to 10 kb 








Light-initiated 
explosives (lead 
or silver azide) 


1/4 sec pulse 


Under development- 
dangerous to work 
with 






Explosive mesh 


5000 taps minimum 

1 ms rise time and 

2 x 10 8 simultaneity 
claimed 
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Table B-7. Testa of Nuclear Radiation Effect* 
Table B-8. Test* of Thermal Radiation Effects 



TABLE B-7. TESTS OF NUCLEAR RADIATION EFFECTS 



EFFECT 
TO BE 
SIMULATED 


SIMULATION 
TECHNIQUE 


ENVIRONMENT OR 
CHARACTERISTIC 


COMMENTS — REMARKS 


AGENCY 
TEST SITE 


Penetration 
of 

Neutrons 
and 

Neutron- 
Induced 
Radiation 


Plasma 

Focus 

Device 

(Zipper) 

Fissionable 
Plate ♦ N 

Source 


Total 10 9 -10 n w/shot 14 mev n 
1/2 sec duration * 

Fission spectrum probably 
largest flux 


Can be repeated at about 1 min 
intervals. Needs development to 
be portable. Kaman Nuclear has 
best capability. 


Kaman 

Nuclear 




Pulsed 
Reactors 


Fission spectrum up to 10 13 n/cm 2 
and up to 10' rad (c) 
SO- 100 sec pulse 


Prop 1963 survey Super Kuka 
Minimum of 10 sec irradiation 


iettelle Memorial 

Institute 




Steady- 
State 
Reactors 


Maximum of 6 x 10 15 n/cm2 sec 


To match 




Penetration 
of Gamma 
Facilities 


Underground 
Facilities 


Maximum of 10 s rad(c)/hr more 
common is 10° rad(c)/hr 


Minimum of 1/2 hr, question of 
spectra must be resolved. 


AEC, NTS 




Irradiation 
Facilities 


500 -Curie source irradiation 
20 in. spec at 10* rads/hr 
500-Curie source irradiation 
20 in. spec at 1.3 x 10 s rads/hr 


Numerous facilities and places 


Martin Co., 
Baltimore 


Permanent 

Neutron 

Daaage 


Pulsed 
Reactors 


lO^n/cm 2 (fission) 600 sec 


Use internal chambers for 
irradiation— largest Super Kuka 
to be operated in 1968. 


LRL, Super 
Kuka 


Transient 

Radiation 

Effects 

Electronic 

System 

(TREES) 


Prompt Gam— 
Simulator 


Optimum specimen size 
Length: SO cm 
Diameter: SO cm 
Target Distance: 75 cm 
Pulse Length: 70 x 10 9 sec 
Max 5000-7000 rads/pulse 
10 in. rad s/ sec 


(Circa 1967) 


USAF/AFSNL, 
Kirtland AFB 



TABLE B-8. TESTS OF THERMAL RADIATION EFFECTS 



EFFECT 
TO BE 
SIMULATED 


SIMULATION 
TECHNIQUE 


ENVIRONMENT OR 
CHARACTERISTIC 


COMMENTS — REMARKS 


AGENCY 
TEST SITE 


Thermal 
Radiation 
on Exposed 
Elements 


Rocket 
Engine 
Exhaust 
Chemical 




Can verify subsystem calculations. 
Simulation not really representa- 
tive in heat/time relations. 


Rocket Test 
Stands 


Thermal 
Radiation 
on Exposed 
Elements 


Thermal 
Simulator 


Thermal Jet 3-5 cal/cm 2 /sec for 
10- cm specimen. 


Not fully developed; requires 
3-5 year development program 
small samples; could provide 
repeatable tests. 


DASA 


Thermal 
Radiation 
on Exposed 
Elements 


Nuclear 
Explosion 
Shock Tube 
(NEST) 


Proposal -stage attempt to 
simulate effects of a surface 
nuclear burst including thermal 
radiation, 1/4 subscale model. 


Thermal radiation is also produced 
in this concept. See detailed 
write-up on 

(Above appraisals circa 1967) 


AEC/ DAS A, NTS 
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Table B-9. Tests of EMP Effects 



TABLE B-9. TESTS OF EMP EFFECTS 



EFFECT 
TO BE 
SIMULATED 


SIMULATION 
TECHNIQUE 


ENVIRONMENT OR 
CHARACTERISTIC 


COMMENTS— REMARKS 


AGENCY OR 
CONTRACTOR 


EMP 

E-Field 
B-Field 


Miniature 
Freae 


Max E Field: 2.5 x 10 s V/a 
Max B Field: 120 Gauss 
Rise Tiae: 0.1 x 10~° sec 
Duration: 5 x 10"° sec 
Spec Size: 5 it dia. x o it 
long 


Saall -scale coaponent or 
subsystem testing 


US Army, ERDL 
Ft. Bel voir, VA 


EMP 

E-Field 
B-Field 


Freae 


Max E Field: 10 s V/a 
Max B Field: 60 Gauss 

(unifora to 10%) 
Rise Tiae: 0.5 x 10 6 sec 
Max Spec Size: 50 ft dia. x 

60 ft long 


Full-scale saaples of subsystea 


US Any, ERDL 
rt. Be l voir, v/\ 


EMP 

E-Field 


ALECS 


Max E Field: 10 4 V/a 
Rise Tiae: 3-10 x 10" 9 sec 
Duration: 100 x 10° sec 


Subsystea or scale saaples 


AFWL, Kirtland 
AFB 


EMP 

E-Field 


ALECS 


Max E Field: 7.5 x 10* V/a 
Riae Tiae: 5-10 x 10' 9 sec 
Pulse Width: 150 x 10' 9 sec 


Subsystea or scale saaples 


AFWL/ L AS L 


EMP 

E-Field 


ARES 

(Advanced 
Research 
EMP 

Simulator) 


Pulse Characteristics- 
Peak Output Voltage: 45 MV 
Energy Storage: 50 kJ 
Field Strength in working 

Voluae: 110 kV/a 
Pulse Rise Tiae: 6 x 10~ 9 sec 
Duration: 100-500 x 10' 9 sec 


Facility has a working voluae 
40 a high, 30 a long, and 
40 a wide. 

Shielded instruaentation rooa 
housed beneath the facility. 


DNA, Kirtland 
AFB 


EMP 


Orange 
Bank 

Generator 


Pulse Length character! sties- - 
Longest: 500- V sec rise 

750- usee decay 
Shortest: 6- usee rise 

100- |i sec decay 
Variable frequency control 
Capacitor discharge 
7 kV - 45 kV 
20 kA aaxiaua 


Saall portable unit 
Coaponent Testing 


AFSWC/SWTVE, 
Kirtland AFB 


EMP 

E-Field 


Marx 

Generator 
(Saall) 


Pulse Length Characteristics- 
Longest: 500- usee rise 

750- usee decay 
Shortest: 6- usee rise 

100- p sec decay 
Variable frequency control 
Gap discharge 1.4 aV - 15 aV 

SO kA aaxiaua 
E Field 10 3 - 10 4 V/a 


Saall portable unit for 
subsystems and coaponent s 


AFSWC/SWTVE, 
Kirtland AFB 

Lightning 
Transient 
Research 
Institute 


EMP 

E-Field 


Marx 

Generator 
(Urge) 


Pujse Length Characteristics- 
Longest: 500- sec rise 

750- sec decay 
Shortest: 6-sec rise 

100- sec decay 
Variable frequency control 
Ignition discharge 
80 - 320 kV 
160 kA aaxiaua 


Large portable unit can be taken 
to sites for full-scale testing. 


AFSWC/SWTVE 
Kirtland AFB 
and Field Sites 


EMP 

E-Field 
H-Field 


Test 

Generator 


H Field - 1 kHz daaped 
since Nave Peaking at 100 
Gauss 

E Free Field - 10 kV/a at 2 a 


In-plant saall saaples 


Boeing Co., 

Seattle, 

Washington 
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TABLE B-9. (CONTINUED) 



EFFECT 
TO BE 
SIMULATED 


SIMULATION 
TECHNIQUE 


ENVIRONMENT OR 
CHARACTERISTIC 


COMMENTS- -REMARKS 


AGENCY OR 
CONTRACTOR 


EMP 


SCREEN 
Rooms ♦ 
Generators 
♦ Elect. 
Test 

Equipment 


Variable, deposing on location 
and contractor 


In-plant testing, subsystems 
and scaled components to full- 
scale components. 


Numerous 


EMP 


Synthetic 
Pulse 
Diagnosis 
(SPUD) 


Measurement of CW transfer 
functions inside an enclosure 
when transmitted at the surface 


Proposal by TRW full-scale tests 
at actual sites 


TRW Systems, 
Actual Test Site 


EMP 

E-Field 


HEMP 


Pulse Character! sties- - 
Max voltage: 400 kV 
Max field: 26 kV/m 
Rise time: 3 x 10" 9 sec 
Decay time constant: 
300 x 10-* sec 


No working volume section; just 
2 transition sections each 
6B m long. Junction of sections 
is IS m high and 24 m wide. 


SAFEGUARD 
Communication 

Agency, 

Fort Huachuca, 
AZ 


EMP 

E-Field 


SIEGE 

(Simulator- 
Induced EMP 
Ground 
Environment) 




Norking volume is 3m high and 
uses the earth as the lower 
plate. Multiple drives, using 
B ft transition. 

Designed to test a buried system. 


SAFEGUARD 

Communication 

Agency, 

Fort Huachuca, 
AZ 


EMP 

E-Field 


TEFS 

(Transport- 
able E lectio 
magnetic 
Field 
Simulator) 


Pulse Characteristics- 
Field strength: SO kV/m 
Rise time: 4 x 10" 9 sec 
Decay time constant: 
350 x 10" 9 sec 


Propagate transient in the 
vertical downward direction. 
Multiple feeds (144) each 
driving 4 transition sections 
(S76 total) 

Can illuminate a 40 m x 40 m area 


SAFEGUARD 
Communication 
Agency, 
Portable, can 
be erected on 
site 


EMP 

E-Field 


Sandia 
Long Mire 


Pulse Characteristics- 
Rise time: 10 x 10' 9 sec 
uuration ito luij . i x iu sec 
Peak field strength: 

* 1000 V/m 
at a point 100 ft from wire. 


Facility is 1000 ft long mounted 
on 40- ft poles. Powered by two 
20 kV power supplies. 


Sandia Corp., 
Kirtland AFB 


EMP 

E-Field 


Martin- 
Marietta 
Long Wire 


Pulse Characteristics- 
Rise time: variable 

5 to 30 x 10' 9 sec 
Pulse width: variable 

100 to 700 x 10" 9 sec 
Max field strength: 

1100 V/m at 100 ft 
Puis* repetition frequency: 

10 pps 


Facility is 1000 ft long and 
46 ft aboveground. 
Powered by two 12S kV power 
supplies. 


Martin-Marietta 
Corp., 
Orlando, FL 


EMP 

E-Field 


AFSMC 
Dipole 
(Long Wire) 


Pulse Characteristics— 
Rise time: 7 x 10~* sec 
Duration: 200 x 10" 9 sec 
Field strength: 

300 V/m at 90 m 
Repetition rate: 10 pps 


Wire is 320 ft long and 45 ft 
aboveground. Located adjacent 
to an aircraft parking pad at 
Kirtland AFB. 


AFSWC, 

Kirtland AFB 


EMP 

E-Field 


HDL 

Biconic 


Pulse Character! st ics— 
Rise time: 4 x 10" 9 sec 
Zero crossing at decay: 

900 x 10" 9 sec 
Field strength (at 90 m): 

4.S kV/m 


Biconic sections 9 ft dia. x 9 ft 
long. Overall antenna length is 
1000 ft mounted 100 ft 
aboveground. 


Harry Diamond Lab a 
Woodbridge, VA 


EMP 

E-Field 


RES 




A flyable (helicopter) Biconic 
radiator. Horizontal version is 
200 ft long, vertical version is 
600 ft long. Design is similar 
to HDL simulator. 


AFWL, Portable 
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TABLE B-9. (CONCLUDED) 



EFFECT 
IT) BE 
SIMULATED 


SIMULATION 
TECHNIQUE 


ENVIRONMENT OR 
CHARACTERISTIC 


COMMENTS- -REMARKS 


AGENCY OR 
CONTRACTOR 


EMP 

E-Field 


NOL/ITTRI 

Hybrid 

Antenna 


Uses a vertical conic antenna 

(50 ft high) and a horizontal 

fringing line (300 ft) 

Pulse Characteristics-- 
Rise time: 10 x 10" 9 sec 
Decay time: 1-100 x 10~ 6 sec 
Field strength: 100 V/m single 
shot, 10 V/m at 60 pps. 


Sub- threat -level facility uniform 
fields developed in the 20 x 50 m 
interaction area. 


NOL/ITTRI 
Crystal Lake, 
IL 


EMI' 

E- Field 


Hybrid 
Antenna 


Similar to NOL/ITTRI, above but 
larger and more powerful. Desigi 
criteria include: 

Length: 1300 ft (fringing line) 
Height: 100 ft (conic antenna) 
Field strength: 

1000 V/m at 300 m 


To be built at Solomons. Maryland 
for Navy use. (comment Circa 
1973) 


NOL, Solomons, MD 


EMP 

E-Field 


DELTA 

Function 

Simulator 


Pulse Characteri stics-- 
Vertically polarized 
Pulse width (at 50% points) 

2 x 10-9 sec 
Field strength: 5 V/m at 5 is 
Beam width: «20 deg. 


Installed on a mountain top at 
White Sands Missile Range to 
evaluate aircraft in flight. 


White Qanrl 
•nine dino 

Missile Range, 

New Mexico 


Oebri s 

Testing 

Components 


CER Filter 
and Dust 
Separator 
Test 


Critical component debris test 


Used to verify design and 
calculations. 


Contractor, 
In-Plant 


Maximun 

Debris 

Load-- 

Debris 

Removal 

System 


Operating 

Test 

Fixture 


Operation of debris removal 
systems under maximum load. 


Laboratory test to confirm design 
of debris removal system 


Contractor » 
In-Plant 


Debris 
Removal 


Operational 

Test to 

Demonstrate 

Debris 

Removal 

System 


Selected debris to simulate 
maximum load. 


Full-scale operational prototype 
system may be conducted in con- 
junction with other closure 
testing. 


Contractor plus 
OCE, Field 
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Table B-10. Chronology o/HEST and DIHEST Tests 



TABLE B-10. CHRONOLOGY OF HEST AND DIHEST TESTS 



DATE 


TEST 


LOCATION 


PURPOSE 


PIT SIZE. 
FT 


Feb to Aug 1964 


Phase I 
(Gaa iag) 


Kirtland API 


Experinont with (1) gaa mixture/water overburden and 
(2) detonating cord/send overburden. Selected latter 
method. 


20 x 40 ! 


Doc 1964 


Phase II 


Kirtland APB 


Determine pressure area and instrument requirements 
for a full-scale Minuteman facility, using a quarter- 
scale model. 


96 x 150 


Feb 1965 


CHBST-2) 


Kirtland API 


Study parameters controlling the HEST air-pressure 
time histories. 


52 x 36 


Mar 1965 


(HEST-3) 


Kirtland APB 


Study parameters controlling the HEST air-pressure 
time histories. 


40 x 48 


Nay 1965 


Phase IIA 


Kirtland AFB 


Double overpressure, change surcharge containment, and 
improve instruments, using same teethed and structures 
as for Phase II. 


88 x 100 


July 1965 


Collins 
Antenna Test 


Kirtland APB 


Using facility similar to Phase IIA, include a cavity 
for investigating response of LEB access and air 
ontra lament systems and antenna. 


40 x 96 


1965 


Parana tor 

Studies 


Kirtland APB 


Study parameters controlling the HEST air-pressure 
time histories. 


Various 


Oct 196S 


(HEST-1) 


Kirtland AFB 


Study parametera controlling the HEST air-pressure 
time histories. 


52 x 36 


Otc 1965 


HEST Test I 
(Quick HEST) 


Warren AFB 
Wing V 


OPERATIONAL TEST: Test sa operational Ninuteman site 
with launch facilities and a ground test missile on 
simulated alert. 


502 x 304 


Mar 1966 


(HEST-6) 


McCo neck's 
Ranch, 


Study free field ground motion. 




Nay 1966 


HIP 1 


Albuquerque 
Kirtland AFB 


Improve HEST environment. 


40 x 60 


June 1966 


HIP la 


Kirtland APB 


Improve HEST environment. 


40 X 60 


July 1966 


HEST Test II 


Warren AFB 
Wing V 


OPERATIONAL TEST: Demonstrate the degree of structural 
survivability of facilities and equipment; assess the 
hardness. 


304 x 352 


Sept 1966 


HEST Test III 


Grand Forks APB 


OPERATIONAL TEST: Substantiate the hardness of tost 
site to meet design attack threat; obtain data for 
extrapolation of Ninuteman force hardness; demonstrate 
missile launch capability after teat. 


302 x 304 


Dec 1966 


Drillhole 


McComack*s 

Ranch, 

Albuquerque 


Study free field ground motion. 


64 x 148 


July 1967 


Backfill 
(HEST-4) 


McComack's 

Ranch, 

Albuquerque 


Study free field ground motion. 


56 x 72 


Oct 1967 


HEST V 

Deaonstration 


Grand Forks APB 
Wing VI 


Demonstrate maximum SOR environment; evaluate surcharge 
disposal; evaluate gage placement techniques; provide 
planning basis for HEST Test V. Used smaller pit. 


64 x 83 


Apr 1968 


Mini-Can X 


Kirtland AFB 


Verify design for Prairie Flat overpressure facility 
for Project LN320; however, shock wave outran detona- 
front and broke detonating cord. 


31 x 75 


Nay 1968 


Mini -Can II 


Kirtland AFB 


Using same site as for Mini -Can I, try 1/16- in. aircraft 
cable support for detonating cord, but also broken by 
shock wave. 


31 x 75 
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TABLE B-10. (CONCLUDED) 



DATE 


TEST 


LOCATION 


PURPOSE 


PIT SIZE, 
FT 


May 1968 


Mini-Can III 


Kirtland APB 


Using sane site as Mini-Can I and II, try two methods 
for protecting detonating cord; one was successful for 
protection, but failed to aeet requirements for shock 
front velocity and impulse. 


31 x 75 


Sept 1968 


HEST Teat V 


Grand Forks API 
Wing VI 


• OPERATIONAL TEST: Determine structural survivability 
and functional capability of launch-essential systems; 
obtain data useful for force hardness assessment. 


300 x 300 


Sept 1968 


Prairia Flat 
HEST Taat 
(Proj. LK320) 


auzzieio mango, 
Canada 


ncal lest wiui mooei structure, oocain rree xieia 
motion and structure response data at specific range 
and overpressure levels to compare with those produced 
by 500-ton Prairie Plat Trial. 




Nov 1968 


Rockteat I 


Albuquerque 


Research in nuclear weapons effects and systems develop- 
ment, as applicable to structures in hard rock. 


160 x 208 


1969 


DATEX I 


Cedar City 


DIHEST Development: Provide information on direct- 
induced rock stresses for HANDEC I experiment. 




Nay 1969 


HANDEC I 


Cedar City 


Demonstrate Improved simulation of a nuclear blast by 
combining a HEST with DIHEST input to ground motion. 


40 x 60 


Aug 1969 


HAND6C II 


Cedar City 


Follow-on to HANDEC I, with reduced overpressure while 
DIHEST was increased nearly 20 times. 


60 x 90 


Oct 1969 


DATEX II 


Cedar City 


DIHEST Development: First employment of a slurry 
explosive as a loading source. 


185 x 200 


Mar 1970 


Rockteat II 


Cedar City 


Test realistic configurations of Minuteman launch 
structures to' survive the HEST-DIHEST environment in 
rock, using ten experimental structures. 


2S0 x 400 



U.S. Army Corps of Engineers 



e. In the testing community the validity of the 
equipment to produce simulated airblast-induced 
ground motions is much debated. The HEST load- 
ing is a moving pressure over a bed with finite 
width and length. The true nuclear-induced loading 
more correctly resembles a line load with a step 
pulse behind the front. The fact that the HEST is 
loading only a small area creates an outward and 
upward motion adjacent to the bed; the true nu- 
clear is essentially a one-dimensional motion. This 
outward and upward motion in HEST may signifi- 
cantly alter the vertical component of motion near 
the center of the bed and consequently there may 
be a low validity for HEST to simulate airblast- 
induced ground motions. 



/. BLEST (Berm Loaded Explosive Simulation 
Technique) is a technique to supplement the HEST 
method by using an array of shallow-buried 
charges in the area adjacent to HEST structure, 
initiated in proper sequence with the HEST Event, 
it generates an extended downward loaded area, 
minimizing, for a while, the undesired outward and 
upward motion in the vicinty of the HEST area 
itself. The design permits simulation detonations 
that allow the test structure to react fully to the 
primary stimulus before boundary effects distort 
the loading. A relatively new addition to the sim- 
ulation techniques, it holds promise of reducing the 
size of the HEST structure itself while improving 
simulation. For further information on the theory 
and on test results, see Schrader et al. (1976) and 
test reports on the HARDPAN I test series that 
used the method. 
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Table B-U. Chronology of High-Explosive Tests 



TABLE B-ll. CHRONOLOGY OF HIGH- EXPLOSIVE TESTS 



DATE 


TEST 


LOCATION 


SIZE OF CHARGE; PURPOSE OF TEST 


Feb 1964 


Flat Top II 


Nevada 


20t TNT, spherical, half buried, desert playa, dry. To develop 
theoretical methods for predicting ground notion and cratering. 


Mar 1964 


Flat Top III 


Nevada 


20tTNT, spherical, half buried, desert playa, wet. 


June 1964 


Flat Top I 


Nevada 


20t TNT, spherical, half buried, limestone, dry. 


July 1964 


Snow Ball 


Canada 


SOOt TNT, spherical, surface. To simulate a Ikt surface nuclear 
burst. 


1966-1967 


Distant Plain: 






1966-1967 


Event 1 


Canada 


20t TNT, spherical, on 85 ft tower. 


1966-1967 


Event 2a 


Canada 


Gas filled balloon, hemisphere on surface. Equivalent to 20t TNT. 


1966-1967 


Event 3 


Canada 


20t TNT, spherical, surface. To study close-in ground shock 
phenomena . 


1966-1967 


Event 4 


Canada 


SOt TNT, Hemisphere, surface. To study blow-down effect in a 
forest . 


1966-1967 


Event 5 


Canada 


20t TNT, spherical, surface, frozen ground. To study cratering 
and ground shock in frozen ground. 


1966-1967 


Event 6 


Canada 


lOOt TNT, spherical, surface. To study scaling factors by com- 
paring results for Events 6 and la, using identical instrument 
patterns. 


1966-1967 


Event la 




20t TNT, spherical, center 29. S ft above ground. To compare 
scaling factors with Event 6 above. 


Aug 1968 


Prairie Flat 


Canada 


SOOt TNT, spherical, surface. To further study scaling factors 
by comparison with Distant Plain Events 6 and la. 


Oct 1968 


Nine Under 


Utah 


lOOt INT, spherical, 2 x charge radius above surface. To study 
cratering and shock effects over rock surface. 


Nov 1968 


Mine Ore 


Utah 


lOOt TNT, spherical, 0.9 x charge radius below surface. To study 
cratering and shock effects over rock surface. 


Sept 1969 


mMiU WW 


Utah 


I4t slurry, sealed cavity, 100 ft below surface. To study ground 
motion with fully contained detonation in rock. 


Oct 1969 


Mineral B/vk 


Utah 


lOOt TNT, spherical, surface. Further study of ejects and motions 
in rock surface. 


Aug 1970 


Dial rSCK 




e/u\. TVT tnkaHnl curfaf* PtilTtiAr c t'll/iv e%t tri 1 i n ff t ACtOTS . 

for comparison with Distant Plain Events 6 and la, and Prairie 
Flat, using same configuration over same geologic medium. 


JUXy 1972 


Mixed 
Company 


Colorado 


)UV\ Bill, Vpnv* , WTIKV* 1 (J • 1 ■ O Kl ,iuiw>u ihm a wm ■ 

burst over a layered geologic medium. 


Sept 1971 


Middle 
Gust I 


Colorado 


20t TNT, spherical, half buried, clay over shale with water table 
10 ft below surface. Middle Gust series to study effect of water 
near surface. 


Dec 1971 


Middle 
Gust II 


Colorado 


lOOt TNT, spherical, 2 x charge radius above surface, clay over 
shale with water table 10 ft below surface. 


Apr 1972 


Middle Gust 
Gust III 


Colorado 


lOOt TNT, spherical, surface, wet, clay over shale with water 
table 10 ft below surface. 


June 1972 


Middle Gust 
Gust IV 


Colorado 


lOOt TNT, spherical, surface, clay over shale with no near-surface 
water. 


Aug 1972 


Middle Gust 
Gust V 


Colorado 


20t TNT, spherical, half buried, clay over shale with no near- 
surface water. 


Spring 1976 


Pre- Dice Throw 
I 1 II 


New Mexico 


120t ANFO (lOOt TNT equivalent) Cylinder plus hemisphere, surface; 
desert alluvium ground- stock measurement and exposure of structures 
and equipment. 


Oct. 1976 


Dice Throw 


New Mexico 


600t ANFO (SOOt TNT equivalent) Cylinder plus hemisphere, surface; 
desert alluvium air-blast exposure of equipment and structures. 



U.S. Army Corps of Engineers 
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B-3. Direct-induced high-explosive 
simulation technique (DIHEST). 

a. The DIHEST development began in 1967 and 
uses a buried array of high explosives to produce a 
specified particle velocity history in the soil at a 
given range from the array. Because of time con- 
straints, the developmental DIHEST study was re- 
stricted largely to simultaneous detonation of rec- 
tangular, planar, and vertical explosive arrays 
(table B-4). 

b. In a buried retangular array of N sperical 
charges, one can argue that the peak horizontal 
particle velocity from the array can never be 
smaller than that obtained from one of the N 
sources (for sufficiently close source spacing). Fur- 
thermore, this lower bound should be approached 
as one moves closer and closer to the array. On the 
other hand, the peak particle velocity away from 
the explosive array should be a function only of the 
array geometry and total yield, and should be in- 
dependent of the number of sources involved. In 
fact, it is expected that a reasonable upper bound 
(approached asymptotically with increasing range) 
is given by the peak particle velocity from a single 
spherical source of the total array yield. 

c. A summary of particle velocity histories mea- 
sured in the vertical symmetry plane at ranges of 
30 and 60 ft from the DIHEST array suggests that: 

—Horizontal particle velocity signatures were 
reasonably consistent within the lower 
two-thirds of the DIHEST depth. 

—Horizontal particle velocity signatures were 
essentially the same form to be expected 
from contained spherical charges (in 
particular, no appreciable shear effects 
were observed). 

—Horizontal peak displacements were consis- 
tent with those expected from a tamped, 
burried, spherical explosive of the same 
yield as the DIHEST array. 

—Significant upward late-time motions asso- 
ciated with free surface effects would be 
produced by the DIHEST. 

d. The DIHEST technique simulates, in a variety 
of geological formations, direct-induced ground mo- 
tions expected from a nuclear device effecting 
buried strategic systems or a part of such systems. 
DIHEST coupled with HEST provides a capability 
of producing airblast as well as direct-induced 
ground motions on a variety of test items. Such 
testing is very expensive ( $2 xlO 6 ), but may be a 
good value as the only viable technique under the 



ban on atmospheric nuclear testing. A chronology 
of DIHEST testing is included in table B-10. 

B-4. High-explosive surface burst. 

a. When nuclear experiments in the atmosphere 
were suspended in May 1963 by a presidential or- 
der preceding the nuclear test ban treaty, it was a 
crucial time for the suspension because a series of 
nuclear and high-explosive experiments were under 
way to compare ground motions from a nuclear 
surface burst with ground motions from a high- 
explosive surface burst. 

6. Flat Top was a series of 20-ton HE experi- 
ments using blocks of TNT with the center of the 
spherical mass at the ground surface. Two experi- 
ments were conducted at an alluvium site and the 
third at a limestone rock site. In alluvium, the 
peak radial displacement measured near the sur- 
face was 24 in. at the 500 psi overpressure station. 
In limestone, the displacements were somewhat 
less (data were not reliable for the experiment in 
limestone). 

c. Over the past decade, many HE field tests 
using built-up charges of high explosives have been 
conducted to satisfy the needs of specific programs. 
These tests often have been conducted in remote 
areas over land and water and in an arctic environ- 
ment. Measurements included transient overpres- 
sure, ground shock and ground motion, water shock 
and water waves, as well as crater size and ejecta 
distribution. 

d. Charge sizes up to 500 tons of TNT have been 
used. There are certain limitations involved with 
even these large-sized HE tests. For example, if the 
charge is detonated on or close to a soil surface, 
the crater formed interferes with the placement of 
buried or surface targets. 

e. Two installations having permanent facilities 
to conduct large HE tests, 100 to 500 ton, are the 
Nevada Test Site (NTS) located near Las Vegas, 
Nevada, and operated by the U.S. Atomic Energy 
Commission; and the Defence Research Establish- 
ment, Suffield (DRES) located near Medicine e 
Hat, Alberta, Canada, and operated by the Defence 
Research Board of Canada. Large HE tests con- 
ducted at remote sites other than NTS and Suffield 
are expensive, and every effort is made to include 
in each test as many of the requirements as possi- 
ble of the various defense and armed forces agen- 
cies. 

/. HE tests to investigate various effects of air- 
blast and ground shock are tabulated in table B-12. 
The tabulation begins with the Flat Top series and 
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TABLE B-12. GRABS TESTS: GIANT REUSABLE AIRBLAST SIMULATOR 



DATE 


TEST 


TEST DESCRIPTION 


July 1971 

through 

1973 




GRABS tests have been conducted at Kirtland AFB, using a buried concrete 
cylinder 18 ft dia. x 30 ft deep, flush with the surface. Test objectives 
were to study structure/media interaction using a stiff structure of 
simple geometry within which well-defined soil media could be provided. 
The GRABS facility is a blast pressure chamber using a HEST-type 
environment. i; 




D-1 


Development test: to test the facility and obtain soils data. 




D-2 


Development test: to test the facility and obtain soils data. 




PI-1 


Two models, stacked one on top of the other, in a dry sand medium. 




PI-3 


Single model, 15 ft high, in dry sand medium. 




PI-4 


Same model as PI-3, in wet sand medium 




PII-1 


20-ft articulated structure, resting on floor of dry sand medium. 




PII-2 


20-ft articulated structure, flush with surface, wet sand medium. 




PII-3 


Upper half of PII-2 structure, resting on floor, dry sand medium. 



U.S. Army Corps of Engineers 



NOTE: This series of tests is expected to continue. 



periment to experiment, however, the data scatter 
when scaled can differ by an order of magnitude; 
testing a prototype facility for defined criteria 
could be very difficult. This data scatter in test 
results may show the uncertainties in present-day 
design criteria. 

c. Tamped nuclear experiments have been recom- 
mend for testing Air Fence Launch Control Facility 
and Launch Facility structures, and may produce 
the best simulation of direct-induced motions. 
Since Calculations substantiate test results for peak 
stress, it is reasonable to expect that tests can be 
designed to produce specified peak stress levels at 
specified structure locations. The uncertainty in ex- 
pected ground motion (primarily displacement) 
must be considered in the design of the test. 

B-6. Underground nuclear tunnel 
experiments. 



goes through the Middle Gust series. 

B-5. Underground nuclear tamped 
burst. 

a. Since the test ban treaty, two important un- 
derground nuclear experiments have been 
performed. HARDHAT and PILEDRIVER. These 
experiments were performed in the granite rock in 
Nevada, at depths of about 1500 ft. A series of 
tunnel liners were located at varying ranges from 
ground zero (GZ), all in the high stress region, 0.5 
kb to 4 kb. (Various liner concepts were used, such 
as liners with backpacking and rockbolting, and 
concrete liners poured directly against the rock.) A 
great deal has been learned about the cavity liner 
design for direct-induced motions. 

b. The experiments have also demonstrated the 
capability to produce large ground motions 
expected from megaton surface bursts. From ex- 
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a. Underground tests of nuclear weapon effects 
require the construction of a tunnel complex for 
installation of the weapon and the test specimens 
and apparatus, combining several experiments to 
get maximum use of the site. For moderate in- 
cremental program cost, structural testing experi- 
ments are added to events primarily designed for 
other experiments. The MIGHTY EPIC/DIABLO 
HAWK pair of events is such an instance. 

6. The MIGHTY EPIC event at NTS was de- 
signed primarily to evalulate the radiation physics 
of a nuclear device using a line-of -sight (LOS) pipe. 
An LOS pipe is a conical tube aligned with its 
projected apex at the center of the nuclear device 
and extending radially outward for as much as a 
thousand feed or more. An LOS pipe is an expen- 
sive device to fabricate; in addition to requiring 
safety closures and test stations in various loca- 
tions, it must also be capable of sustaining a hard 
vacuum on the inside (required to prevent attenu- 
ation of radiation along the interior of the pipe). 
MIGHTY EPIC is a landmark event to designers of 
hardened structures because of other programs 
that were conducted at the same time. 

c. For the structures program, several types of 
cylindrical and spherical structures were tested. 
Three separate drifts were constructed at ground 
ranges from the nuclear device that would produce 
the design stress level, one-half design level, and 
twice design level. One of the novel features of this 
experiment was the placement of nine copies-as 
identical as could be constructed— of a spherical 
reinforced concrete structure. They were installed 
in sets of three at each stress level. It was expected 
that the close-in models would fail, the mid-range 
models would sustain moderate damage, and the 
distant units would be unscathed. In general, the 
actual test results followed theory, and the replica- 
tion of test structures at each stress range gave a 
data base for evaluating the random variations of 
design parameters. Additional benefits will accrue 
from the DIABLO HAWK event (see e below). 

A Another objective of the MIGHTY EPIC event 
was to test theory that structures sited in a hard 
medium overlaid by an extensive layer of more 
porous medium would experience a more moderate 
environment than if they were sited at an equal 
depth in an all-hard medium. MIGHTY EPIC was 
sited in a quartzite bed overlaid with tuff. A useful 
amount of data was collected with which it is 
possible to evaluate the concept. The serious re- 
searcher is encouraged to seek out the data in the 
classified MIGHTY EPIC progress reports. 



e. DIABLO HAWK, a follow-on to MIGHTY 
EPIC planned for the same test bed, is designed 
and will be located to provide a repeat loading of 
the structure tested in the MIGHTY EPIC event. 
This second loading of the structures will provide a 
stress wave propagating approximately 90 deg to 
the original MIGHTY EPIC stress wave direction. 
Such cyclic loading is rare in dynamic testing, out- 
side of small-scale models in laboratories. Access to 
the structures was gained by retunneling into the 
test bed and opening the access ports. Extensive 
photographic coverage and dimensional mapping 
were performed after the MIGHTY EPIC event and 
have become "pretest measurements" for the DIA- 
BLO HAWK event, scheduled for some time in 
1978. 

B-7. Giant reusable airblast simulator 
(GRABS). 

a. The GRABS is a cylindrical silo, 18-ft dia. by 
48-ft deep and lined with reinforced concrete, con- 
structed in a massive limestone formation at Kirt- 
land AFB. New Mexico. It was developed to provide 
a two-dimensional zxisymmetric test facility for 
evaluating the computational predictions of the 
vertical motions of missile silos from overpressure. 

6. A preselected amount of soil (up to 30 ft) is 
placed in the lower portion of the facility along 
with instrumentation and model structures, and a 
modified HEST is used to generate overpressure 
within the facility. This technique has the dem- 
onstrated ability to reproduce reasonably the air 
pressure history up to the -1800 psi overpressure 
level from a 1-Mt surface burst. For a predeter- 
mined time, the technique confines within the fa- 
cility a high-explosive detonation in a finite vol- 
ume. This is done by using soil as a surcharge 
above an explosive cavity. The peak pressure, ini- 
tial decay, and late-time decay can be controlled by 
variations in explosive charge density, initial explo- 
sive cavity volume, and the height and density of 
the surcharge material, respectively. The detona- 
tion projects the surcharge into the air, thereby 
eliminating the need for reaction devices. The sur- 
charge support system is designed so that no 
braces extend to the test-bed; thus only the detona- 
tion shock is transmitted to the test-bed. Table 
B-12 is a list of tests conducted in the GRABS 
facility. 

B-8. EMP simulation methods. 

a. The EMP associated with a nuclear explosion 
can be considered an insidious effect, for it can do 
damage to unprotected systems at distances from 
the blast where other effects are not bothersome 
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(high-altitude bursts, for instance). However, its 
effects can be guarded against, and its consider- 
ation should be carried throughout the design and 
construction phase to achieve the best results for 
the least cost. Retrofitting EMP protective devices 
to a finished facility can be very expensive. Testing 
for design and hardness verification should be done 
at each significant step of the way. A summary of 
available facilities and techniques follows. 

b. In the last decade, a large amount of research 
has gone into EMP and its cause, effects, and miti- 
gation. Three documents of major significance to 
designers of hardened facilities summarize the 
state of the art quite well: GE-TEMPO, 1974 a; 
Schlegel, 1972; Bridge-Emberson, 1974. The follow- 
ing sections are extracted from the DNA EMP 
Handbook because it sets forth cogently the idea 
that attitude and prior analysis are so vital to the 
conduct of meaningful tests; less dedicated testing 
becomes mere busy work. 

c. In general, EMP testing is expensive, time 
consuming, and filled with pitfalls that make it 
difficult to obtain credible results. Factors involved 
in deciding whether to engage in a test program 
and how to set it up include: 

— Priority of assessing vulnerability of the sys- 
tem under consideration 

—Simulation requirements based on system defi- 
nition and postulated threat 

—Costs involved 

—Real-time requirements 

Once the decision to test has been made, its success 
or failure depends, to a great degree, on the exper- 
ience and attitude of the individuals directing it. A 
good test is one that 
—Is viewed with a hypothesis at hand, derived 
by analysis; 

—Includes within its concept of accomplishment 
a statistical view of the problem; 

—Is conducted by experimenters who are con- 
stantly prepared to reject their hypothesis 
during the test, based on the evidence at 
hand, and who are willing to develop and 
investigate a new hypothesis if required. 

d. A test should constantly be monitored by 
qualified analysts. For large-scale tests the an- 
alysts should be on-site as well as off -site to assess 
the progress of the test as it proceeds, redirect it if 
necessary, and see that the results are credible. 
The importance of analysis before the test cannot 
be overstressed, because it forms the necessary ba- 
sis for continuous evaluation in terms of success or 
failure as the experiment proceeds. The required 
depth and extent of this analysis are matters of 



judgment, but even rough approximations are use- 
ful. Analysis is extremely difficult if it attempts to 
achieve high accuracy; even a transferfunction pre- 
diction for a real system can prove to be a chal- 
lenging exercise. However, gross estimates must be 
made to determine feasibility in reasonble time and 
cost; these are based on generic groupings of the 
system and parametric analysis of these groupings. 
In this type of analysis, there is no substitute for 
experience; the experienced analyst has at his fin- 
gertips a storehouse of analytical tools, experimen- 
tal data, and knowledge of what has and has not 
succeeded in the past. This perspective is a vital 
factor in achieving reasonably analytical results in 
a short time. 

e. In EMP tests, nothing is precise. Neither the 
simulation, nor the instrumentation, nor the data 
reduction, nor the system behavior itself is ab- 
solute. All vary, and many introduce errors. It is 
essential that this be recognized early in the pro- 
gram. The fact that precision is impossible should 
not be interpreted as an excuse for carelessness. To 
the contrary, every reasonable effort should be 
made to minimize errors everywhere so that at the 
end of the program an assessment can be made 
that is as free from uncertainty as possible. 

/. The problem must also be viewed statistically, 
because the system under test may vary from unit 
to unit. Early recognition of this system variation 
will permit a more orderly pursuit of the end goal. 
If errors and system variability are not bounded, 
the assessment of vulnerability can be totally in 
error. If errors and uncertainties are not mini- 
mized, the system may be penalized by excessive 
hardening requirements to cover not only the EMP 
effect but the large errors and uncertainties asso- 
ciated with the assessment of that effect. 

g. It is important to think carefully at the begin- 
ning of test planning about what is desired as an 
output of the test and exactly how that result will 
be used in the assessment. Voltage and current 
waveforms might be collected (time domain) at low 
level with the idea of using linear extrapolations in 
time domain are invalid if the system is nonlinear. 
To cover this uncertainty, a high-level simulation 
using a reasonable pulse waveform is mandatory. 
At the same time, an analytical model of the sys- 
tem is highly desirable. Alternatively, transfer 
functions from external fields to voltages (for cur- 
rents) in critical locations might be obtained 
(frequency-domain) with the idea of using these in 
analytical predictions. However, accurate transfer 
function determnation may require a simulation 
experiment different from that used in voltage or 
current determinations, and again the validity of a 
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transfer function depends on system linearity. 
Thus, obtaining both representative timedomain re- 
sponses and system transfer functions with one 
simulation may not be compatible goals, and it is 
important to know what is needed and to have an 
alternative if the system proves to be nonlinear. 

h. If the system behavior is truly nonlinear at 
the higher levels, then probably the only way to 
predict this behavior is by developing an analytical 
system model, including models of the characteris- 
tics of the important nonlinear devices or system 
segments. From this model, valid time-domain pre- 
dictions can be made on a computer. 

i In general, the problem of assessing a system's 
vulnerability to EMP is not unlike a massive 
trouble-shooting problem in which the system may 
have several faults that must be located and cor- 
rected. The process is inherently an iterative one in 
which the experimenter learns, little by little, what 
the causes and effects are in the system and how 
detrimental effects can be prevented. Results to- 
ward this end will be achieved sooner if the data 
are collected carefully right from the beginning, 
since poor data can lead to erroneous conclusions, 
confusion, and delays. 

j. The two approaches that can be used in at- 
tempting to establish criteria for what should be 
measured are to: 

Work directly on circuits believed to be most 
critical and measure the input voltages of 
these circuits (a hazardous approach) 

Work into the system from the outside, measur- 
ing first what the EMP induces on the 
exterior of the system, then systematically 
tracing the current flow from there to cab- 
ies, and along these cables to subsystem 
black boxes 

k. The direct circuit approach is plagued by 
many traps. Voltage at the input of a particular 
circuit ignores the fact that the EMP may be com- 
ing in on a different circuit-interface wire such as 
ground or power or some other lead. To evaluate 
the "input" of EMP to that circuit, it would be 
necessary to measure the current or voltage on 
every wire leading to it. This may so load down the 
circuit (probes have capacitances that may not be 
negligible in their effects at these frequencies) that 
the data are useless. Connecting several probes into 
a circuit without bringing in some EMP currents is 
also a difficult matter. One of the greatest 
weaknesses of this approach, however, is that it 
may provide the semblance of data assessment 
without any real understanding of why the signal 



is at the circuit, why it looks as it does, and how it 
can best be eliminated. This type of premature 
assessment of systems is hazardous. 

/. When a system is worked into from the out- 
side, an overall understanding of system behavior 
is obtained sooner. At the same time, enough data 
can be obtained at subsystem or black-box inter- 
faces to permit more detailed investigations to be- 
gin in the laboratory. Though there are some ex- 
ceptions, this procedure is usually best. 

m. What should be instrumented are those 
things that will answer these key points: 

— How much current is induced on the exterior 

of the system. 

— How much of this external current gets inside 

the system. 

— How this energy gets inside. 

— Where this energy goes and how it is distrib- 

uted. 

— What this energy does to each subsystem. 

— How any detrimental effects can be elimi- 

nated. 

n. These questions should guide the selection of 
points to be measured. Having selected all of the 
desired measurement points, it is probable that 
there will be too many to instrument at one time 
because of internal space limitations or other con- 
siderations. In this case, the measurements are 
best performed a few at a time, because instrumen- 
tation does load the system and change current 
amplitudes and distributions to a certain extent, 
and thus should be kept to a minimum. 

o. Electromagnetic scale modeling is an impor- 
tant alternative to full-scale testing under the fol- 
lowing conditions: 

— Test facilities or available equipment are at a 

premium. 

— The system to be tested is very large. 

— The system dedication cost for full-scale test- 

ing is high. In 

addition to the advantages of modeling under these 
conditions, benefits can be derived as follows: 
—Sensors can perhaps be better placed during 
full-scale testing as a result of model ex- 
periments. 

—Design modifications or cable reroutes can be 
made prior to fullscale testing. 

—Electromagnetic angles-of -arrival can be deter- 
mined for worst- and best-case conditions. 

—Effects of change in the conductivity of sur- 
rounding media can be explored to an ex- 
tent not possible in full scale. 

—Estimates can be made of some of the re- 
sponses of a complex system prior to full- 
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scale testing. 
—Quantitative data can be obtained to validate 
analysis. 

p. It should be pointed out that because of the 
difficulty in introducing minute openings or poor 
bonds into models, and since these often control 
inferior fields, the usefulness of modeling is or- 
dinarily limited to the measurement of limited 
value, and are generally appropriate only in con- 
firming previous analysis. However, once the ex- 
terior fields, voltages, and currents are known for a 
complex structure, perhaps having cable runs, ana- 
lysis can often yield internal field quantities of 
interest. 

q. In actually setting up a scale model test, the 
following should be kept in mind. 

(1) Broadband pulse response determination in- 
volves much more than does a steady-state, single- 
frequency response test. 

(2) Special electromagnetic illumination sour- 
ces are required that are coherent, have uniform 
time delay, and use antennas with constant effectiv 
height. 

(3) Special modeling techniques are required 
for studying exposed conductors passing over or 
within a lossy dielectric, such as earth. A pulse- 
type waveform can theoretically be replaced by a 
continuous wave (CW) source with a sensing sys- 
tem that references the sensed CW signal to a 
reference phase from the source. Complex Fourier 
transfer functions can be developed by computer 
processing the recorded data. However, long sweep 
times are required to ensure that all narrow band 
responses are adequately explored, and the actual 
physical implementation of such an approach in 
the microwave band poses additional difficulties. 
On the other hand, the use of scaled real-time 
waveforms allows quick development of actual re- 
sponses, from which complex Fourier transfer func- 
tions can also be developed with the aid of comput- 
ers. 

r. Scaling relationships are derived from the the- 
ory of electrodynamic similitude. We define the 
modeling factor, M, by: 



where 

D a = A dimension of the actual system to be 
modeled 

D,=The same dimension of the scaled system 
B-22 



For example, when a 300-ft long structure is scaled 
down to a 30-ft long model, this 1/10 scale model 
has a modeling factor M 10. 

s. The relations between scaled and actual quan- 
tities are: 

Model Size D.= d./M (B-2) 

Frequency a>,= Ma> a (B-3) 

Conductivity o-, = Mo; (B-4) 

Dielectric Constant (B-5) 

Permeability M t = M» (B-6) 

t A plane wave propagating in the positive z 
direction in an imperfect dielectric can be char- 
acterized by H(z,t)=H m e- Mi(wt - x) (B-7) 

where H m is the wave amplitude at z 0, a describes 
the wave attenuation in the dielectric, and B is the 
phase constant, a and B are real and given by: 



a = w|m*/2^ 

0= wJ^€/2^: 



1 + aVwV +1 



1/2 



1/2 



(B-8) 



(B-9) 



The wavelength, X , in the dielectric is given by 



X=2tt//9 



(B-10) 



and from equations B-8 through B-10, it is seen that 
o , 0 , and X scale as follows: 



a,=Ma. 
0.=M0. 
X.= X./M 



(B-ll) 
(B-12) 
(B-13) 



u. A summary of EMP simulation test facilities 
is shown in table B-9. The listing is not all- 
inclusive by any means but does give an indication 
of the types of fixed facilities (relative) that are 
available. In the extensive literature available, 
methods are discussed for establishing test arrays 
at the facility to te tested. 

v. Note that in the tables the simulator type 
includes stationary or fixed simulators, and porta- 
ble ones. The only simulator listed that can be 
moved about and is not fixed is RES-I, which is 
helicopter-transportable. Pulse variability figures 
represent uncertainties in amplitude and decay 
times, and jitter represents the uncertainty in the 
firing time. Cycle time is the time between pulses. 
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B-9. Blast simulation techniques for 
testing air-entrainment systems and 
blast closures. 

a. Comments are limited to methods available 
for testing overpressure on the performance of har- 
dened air entrainment systems and blast closures 
under conditions simulating a nuclear attack. Tests 
specifically intended for the assessment of debris 
effects do not presently exist; therefore, comments 
in this area are directed to desirable tests and 
possible test techniques. 

6. Shock-tube, implosion driver, and underground 
field tests have been instrumented to test overpres- 
sure effects on blast valves and air-entrainment 
systems (tables B-l through B-4). The blast simula- 
tor at the WES is capable of producing environ- 
ments with an upper pressure limit of 1000 psi; the 
physical size of the test chamber requires subscale 
models. The current pressure rise time is much too 
low for meaningful air-entrainment system tests, 
and it is not considered a useful facility for testing 
present or future air-entrainment systems. 

c. Surface debris entering air-entrainment sys- 
tems during a nuclear attack will come from four 
sources: 

— Particulate material entrained in the air 

shock flow 

— Surface material from cratering moved past 

the entrance 

— Crater ejecta descending from aerial trajec- 

tories 

— Particulate fallout from the nuclear cloud 

d. Surface debris from overpressure flow will be 
relatively small particles carried into the air en- 
trainment systems by the flow velocity behind the 
air shock. This material will be distributed 
throughout the air entrainment system, although 
very little material should reach the vicinity of the 
blast valve. Test simulation of this debris will be 
very difficult because a method of entraining the 
proper particle size and velocity distributions is not 
available. Some of this debris will be removed from 
the facility during the negative overpressure phase, 
and it is not expected that the residue will con- 
stitute a threat to the facility. Future improve- 
ments in the analysis of two-phase shock flows 



may lead to tests that might accurately assess the 
distribution of this material in air entrainment 
systems (table B-10). 

e. Cratering causes most of the debris impinging 
on facilities. For sufficiently large weapons and 
close detonations, the air-entrainment system sur- 
face entrance may be engulfed in debris material 
from cratering. A substantial amount of this debris 
may collect in the debris pit of the air entrainment 
systems, but very little should enter the blast-valve 
trigger or delay lines. 

/. Simulation of this debris flow is complicated 
by lack of test data from nuclear blasts, and lack 
of adequate scaling principles from HE to nuclear 
blast. Tests of some value could be performed by 
using large chemical explosives in close proximity 
to full or scale models of air-entrainment systems. 
Similarity of test site media to that of proposed 
facility locations would be important, since scaling 
of explosive effects between various types of soils is 
not well understood. Tests of this type would be of 
use in determining the likely distribution of var- 
ious sizes of debris particles, the nature of bloc- 
kages of the surface entrances, and methods of 
eliminating such blockages. To reduce the expense 
of these tests, combine them with other large HE 
programs, such as those conducted at DRES, Al- 
berta. 

g. Debris from the nuclear cloud will settle into 
the facility for some time after an attack, but will 
be concentrated primarily in the debris pit. Fallout 
distribution was studied in some detail at the time 
of surface nuclear tests, and can now be predicted 
with some assurance. A test that might be ap- 
plicable for assessing the distribution of fallout 
material in an air-entrainment system could be 
performed by allowing equivalent size particles to 
settle into the entrance while the ventilation sys- 
tem was operating at a level appropriate to postat- 
tack conditions. 

h. The distribution and importance of initial at- 
tack debris could be significantly altered by subse- 
quent blasts. If multiple attacks are part of the 
design threat, place debris equivalent to that from 
an initial attack in an air entrainment system 
prior to performing a blast simulation 
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GLOSSARY 



Acceptance 
Region: 



Air-Entrainment 
System: 



Bias Errors: 



Binomial 
Distribution: 



Blast Attenuator: 



The range of a statistic (in this 
text, the t-statist ic) for which 
the hypothesis is accepted as 
true. 



Accomplishes a continous or pe- 
riodic transfer of air (gas) be- 
tween the atmosphere and the 
facility. 



Systematic errors resulting 
from imcomplete or inaccurate 
modeling, measurement errors, 
etc. 

A distribution where the values 
r/n for r 0, 1, N have the 
frequencies 



n! 



r!(n-r)! 



P r (l-P) n " r 



where p is the probabilty occur- 
rence. 

A device for reducing in an air 
entrainment or exhaust duct. 
Wall friction, flow restriction, 
and expansion chambers are of- 
ten used 



Confidence The range within which an es- 

Interval: timate will fall for a specified 

level of confidence, e.g., 90% of 
the time the difference between 
the mean of a normal popula- 
tion and the mean of an n- 
number sample will fall 
between ± 1.96 o-/y/n. 



Continuous A steady state excitation signal, 

Wave: e.g., cos cot, as opposed to a 

transient or impulsive 
excitation. Frequency sweep or 
frequency stepping may be used 
to cover the broad frequency 
range of interest. 



Correlation The ratio of the covariance of 

Coefficient: two random variables to the 

square root of the product of 
their variances, i.e., 



E [(X-Mx)(y-My)] 



Blast Valve: 



Burst Conditions: 



Chi Square (X 2 ) 

Goodness-of-Fit 

Test: 



Confidence 
Coefficient: 



A valve that prevents entry of 
overpressure into hardened fa- 
cilities. 

A description of the location of 
point of burst relative to the 
ground surface and to the tar- 
get. 



A means of measuring the dis- 
crepancy between the Probabil- 
ity Density Function (PDF) ex- 
hibited by a data set and that 
for Gaussian, i.e., normal, or 
other hypothetical distribution. 



The confidence level associated 
with the confidence interval. 



Delay Line: That portion of an air- 

entrainment system between a 
sensor and blast valve that al- 
lows valve closure before over- 
pressure arrival. 

Deterministic: having an assigned value; with- 
out uncertainty. 

DIHEST: Direct-Induced High-Explosive 

Simulation Technique. A 
method using a buried array of 
high-explosives to producje a 
specified particle velocity time- 
history in the soil a given dis- 
tance from the array. Coupled 
with HEST it can produce air 
blast and direct-induced ground 
motions on various targets. 
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Direct Injection: 



DRES: 



EMP: 



Expansion 
Chamber: 



Failure Modes: 



Failure-Oriented 
Analysis: 



Glossary- 2 



Excitation of a system or sys- 
tem element by direct applica- 
tion of electromagnetic signal, 
force, or displacement rather 
than by the application of such 
loads that pass through and are 
modified by intermediate 
elements such as protective 
structures and shock-isolation 
system or shielding. 

Defense Research 
Establishment, Suffield. A 
semi-remote research facility 
near Medicine Hat, Alberta, 
Canada, where experiments re- 
quiring extensive land area can 
be conducted. Site of previous 
100 ton and 500 ton high explo- 
sive tests. 

Electromagnetic Pulse. 
Associated primarily with the 
high intensity radiation and 
conduction fields induced by nu- 
clear explosions. Can produce 
extremely high currents in con- 
ducting element, disrupting or 
distroying electronic 
components. 

A type of blast attenuator for 
air entrainment/ exhaust 
systems which relies on the 
pressure-averaging effect of 
chamber fill time to mask or 
attenuate the maximum over- 
pressure associated with sharp 
peaks 



Identifiable mechanisms of sys- 
tem or system element failureks 
e. g., communication system 
failure modes may include blast 
damage to antenna, EMP Bur- 
nout of antenna lead, ground- 
shock damage to antenna lead, 
EMP burnout of transmitter, 
shock and vibration damage to 
transmitter, power outage, etc. 

An analysis that examines 
occurrence/ nonoccurrence of 
failures capable of preventing 
mission function. 



Fault Trees: 



Functional Block 
Diagram: 



GRABS: 



HARDHAT: 



Hardness 
Compliance: 

Hardness 
Verification: 



HEST: 



A graphical presentation of the 
relationship between loads and 
failure modes and of their rela- 
tion to system composition from 
the element through overall 
system levels. 

A pictorial presentation of the 
operating relationship between 
system elements. 



Giant Reusable Air-Blast Sim- 
ulator. A large Cylindrical silo 
in massive limestone at Kirt- 
land Air Force Base, New 
Mexico. Designed to provide a 
2-dimensional axisymmetrical 
test bed for air blast overpres- 
sure stimulus. 

A code name for an 
underground tamped nuclear 
event to evaluate hardened 
structural designs. 

Satisfaction of hardness 
requirement. 

The process of determining that 
a system or system element has 
at least the resistance 
(hardness) claimed for it. 



High-Explosive Simulation 
Technique. A method for simu- 
lating nuclear air-blast-induced 
ground motion. 



Hydraulic Surge: Water hammer. 



Ignorance 
Factor: 



Impendance: 



A measure of that portion of 
the uncertainty resulting from 
incomplete knowledge of the 
phenomena or from bias or sys- 
tematic error in measuring the 
phenomena, the uncertainty not 
related to the truly random 
natur of the phenomena. 

Mechanical impendance, the ra- 
tio of the acceleration of a me- 
chanical system to the 
sinusoidal force exciting it 
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Intake Structure: 



Level of 
Significance: 



That portion of an air- 
entrainment system where air 
enters the system. 

a, the chance of rejecting a true 
hypothesis; the complement of 
the level of confidence. 



Network Logic 
Diagrams: 



A graphical means of present- 
ing the functional relationships 
of system elements, loads, and 
failure modes; similar to flow 
diagrams used in computer pro- 
gramming. 



Lognormal 
Distribution: 



A statistical distribution in 
which the logarithm x of the 
variable y is normally distrib- 
uted; its probability density 
functions: 



Nonparametric: 



*(y)- 



Low-Level 
Transient: 



where o-^ ny is the standard de- 
viation of n y and my is the 
median value of y. 

Excitation by a transient signal 
simulating the threat time his- 
tory but at amplitudes below 
those specified for the threat 
level. 



Normal 
Distribution: 



*<x). 



Used here to denote those situ- 
ations where either (1) the dis- 
tribution is not a member of a 
known class (normal, lognormal, 
binomial, etc.) or (2) there are 
not enough data to allow iden- 
tification of the distribution. 
Under these conditions 
estimates of parameters such as 
the mean and variance cannot 
be obtained. However, 
confidence intervals for quantit- 
ies can be obtained. 

Distributed in a Gaussian man- 
ner; i.e., having the probability 
density function 



-L exp -(x-m)V(2ct 2 ) 



Minimum Sample 
Size: 



Mission Critical 
Functions: 



Monte Carlo 
Method: 



Mounts and 
Fasteners: 



The smallest sample from 
which a valid statistical infer- 
ence can be made. 



Those functions that are neces- 
sary to the execution of a sys- 
tem primary mission. 



A technique that obtains pro- 
babilistic approximations to 
problems by executing a large 
number of simulation problems 
with parameters defined by sta- 
tistical sampling. 



The mechanical components 
used to connect equipment to 
protective structures, platforms, 
racks, shock-isolation systems, 
etc. 



NTS: 



One-Sided 
Procedure: 



Overburden: 



Penetration: 



PILEDRIVER: 



where m =mean and <r 2 
riance. 



=va- 



Nevada Test Site, under the 
control of DOE. 

A statistical test wherin the hy- 
pothesis will be accepted if the 
statistic satisfies a one-sided 
condition of the form t x. 



Overlying soil or rock. . 

An opening that pierces the 
protective shell of a hardened . 
facility, such a a conduit for 
communication or power cables. 

Code name for an underground 
tamped nuclear event to evalu- 
ate concepts of hardened struc- 
ture design. 
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Population: 



Probabilistic 
Assessment: 



Probability 
Density Function: 



Probability 
Distribution: 



Probability of 
Success: 



Protective 
Facilities: 



Pulse-Train 
Simulation Test: 



Radiation 
Shielding: 



The set of objects or measur- 
able effects having some com- 
mon observable properties; for 
example, all Minuteman Launch 
facilities, or all thermo-nuclear 
explosions. 

methodology for addressing the 
influence of uncertainties in 
loads and resistances when as- 
sessing the hardness and 
survivability /vulnerability of 
hardened facilities. 



A mathematical statement of 
the frequency of occurrence of 
possible values of a random 
variable (see, for example, Nor- 
mal Distribution). 

A function that assigns to each 
possible value of a random vari- 
able the probability of its occur- 
rence. For continuous rando 
variables, the integral of the 
probability density function, 
also known as the cummulative 
distribution function. 



In this text, probability that a 
resistance exceeds its 
corresponding load, i.e., the 
probability that failure will not 
occur. 



Facilities whose function is to 
protect material, equipment, 
personnel, and mission capabil- 
ity from the harmful effects of 
nuclear weapons. 

A test in which system response 
is excited by a train of pulses 
designed to simulate the threat 
driving forces. 



Material that prevents penetra- 
tion or reradiation of nuclear 
radiation environments— a, ff , 
and T-radiations from nuclear 
weapons. 



Random Errors: 



Random 
Uncertainties: 



Repetitive Pulse: 



Resistance: 

Resistance 
Design Goal: 



Shock-Isolation 
System: 



Stationary Field: 



Survivability: 



Survival 
Probability: 



System: 



System 
Engineering: 



Errors distributed according to 
chance, as opposed to bias or 
systematic errors. 

Uncertainties in the value of an 
attribute resulting from its ran- 
dom nature. 

Excitation by a pulse train 
which covers the frequency 
range of the threat environment 
but may not simulate the time 
history of the threat pulse. 

Ability to withstand nuclear 
weapon-effect loads. 

The level of resistance to be 
achieved to satisfy survivability 
requirements. Often stated in 
terms of local free-field nuclear 
weapon-effect environment am- 
plitudes. 

A system— mechanical, hydrau- 
lic, pneumatic or hybrid— that 
attenuates the shock and vibra- 
tion environment transmitted 
throug its elements. 

An EMP simulation for small 
components not sensitive to the 
transient nature of threat EMP 
pulses. 

The probability that a 
facility/Subsystem/component 
failure-mode will functionally 
survive a nuclear-weapon attack 
and retain its physical integrity 
during the specified endurance 
period. 

The probability that system or 
system element resistances ex- 
ceed their corresponding loads. 

A combination of elements, sub- 
systems, or pieces of equipment 
integrated to perform a specific 
function. 

A science dealing with the de- 
sign and performance of inter- 
connected components, sub- 
systems, and systems. 
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System Hardness 
Level: 



Systematic Error: 



Systematic 
Uncertainties: 



Tap: 



Thermal 
Shielding: 



TTCP: 



The maximum load level for 
which a system retains 
functional capability at a pre- 
scribed level of confidence. 



An error resulting from bias in 
measurement, and not from 
chance. In this text, ignorance 
or imcompleteness of knowledge 
is treated as a systematic 
(nonrandom) error. 

Uncertainties due to unknown 
systematic errors or ignorance 
(incomplete knowledge). 



1 tap — 1 dyne-sec/cm 2 



1 ubar-sec 



Two-Sided 
Procedure: 



14 Xl0 _6 psi-sec 



Material that provides protec- 
tion from the thermal 
environment— radiation and 
fireball immersion— of nuclear 
weapon explosions. 



Tripartite Technical Coopera- 
tion Programme. Member 
nations are Canada, the United 
Kingdom, and the United 
States, s. 



Type I Error: 
Type II Error: 
Uncertainty: 

Variance: 



Weapon Range: 



A statistical test wherein the 
hypothesis is accepted if the 
statistic falls between limits, 
i.e., ti <t <t2, in contrast to a 
one-sided procedure where the 
statistic must satisfy an 
inequality of the form t <ti (or 

t >t2). 

Statistical probability of reject- 
ing a hypothesis when it is true, 
i.e., an a error. 

Acceptance of a hypothesis as 
true on statistical grounds when 
it is false, i.e., a B error. 

The amount (estimated) by 
which the predicted value may 
vary from the observed or true 
value. 

The square of the standard de- 
viation; for a finite sample the 
variance (S 2 ) is defined by 



_ ^ (Xj-X) 2 

^ N-l 



S 2 



i=n 

The horizontal distance from 
the burst point to the target 
point, also called the "offset." 
The distance between burst 
point and target point is the 
slant range. 



Threat Scenario: 



A description of the expected 
nuclear attack, including num- 
ber of weapons, their yields, 
burst conditions, and timing of 
their detonation. 



Transfer 
Function: 



A function relating free-field 
nuclear weapon -effect 
environments to local loads on a 
system or system element. 



t Statistic: 



A statistic for comparing the 
sample and population means 
when the standard deviation is 
unknown. 
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QUESTIONNAIRE 



1. If you want to receive future changes to this manual or to correct your address, return this page to: 
Commander 

US Army Corps of Engineers 
ATTN: DAEN-ECE-T 
Washington, D.C. 20314 



2. Do you have any other questions, comments, or suggestion? (Department of Army respondents: Use DA 
Form 2028) 



3. Your Address; (Attach marked mailing label with corrections indicated) 



4. Your Telephone Number: 
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(DAEN-ECE-T), WASH DC 20314 
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